四阶龙格库塔法解微分方程组 Mathematica

[x,y]=dsolve('Dx=x-y-x*(x^2+y^2)','Dy=x+y-y*(x^2+y^2)','x(0)=2','y(0)=1')得到的结果是解析解没有找到.用数值解.在Matlab下

在Matlab下输入：edit,然后将下面两行百分号之间的内容,复制进去,保存%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

告诉我邮箱,发你代码

注意,12.19t这种写法不对,应该写成12.19*t.结果如下>>symsS1S2t>>[S2,S1]=dsolve('DS2=12.19*t-24209-0.3*S2-0.01*S2','DS1=

把方程和初值条件都列出来这一大堆东西几个人看得懂?再问：公式（1）（2）就是方程啊，我写的那两个就是对应的方程。倒数第二张图片有实验数据。还需要什么参数吗？非常抱歉，我学的专业不太扎实，问题可能有点弱

常微分方程的解有标准格式,需要根据议程形式假设其通解,比如假设y=c1x+c2x'+c3,代入方程解通解,再加边界条件x(0)=y(0)=1解特解.

怎么说呢...可能是因为所求系数(整型)太大了,也可能是你初始量设置的有问题.具体这会儿我没参考书,要不然就告诉你怎么改进.那句"Error,(incombine/radical)integertoo

functionhhX0=[17,195,1500,4,62,200];TSPAN=[04];[t,X]=ode45(@myfun,TSPAN,X0)plot(t,X)legend('x','y','

dsolve('Dx1=-8/3*x1+x2*x3','Dx2=-10*x2+10*x3','Dx3=-x2*x1+28*x2-x3')ans=[x3(t)=0,x2(t)=0,x1(t)=C1*ex

解方程很简单,调用ode45即可关键是求运动时间,假设z方向坐标变为0即为落至地面运动结束,可求得运动时间的近似值及落地点坐标代码如下：clearallclcc=0.07;u0=25;a=pi/6;f

dsolve('m*D2y-a-b*f*Dx-sqrt(Dx*Dx+Dy*Dy)*Dy','m*D2x-b*f*Dy-c*sqrt(Dx*Dx+Dy*Dy)*Dx')正常情况使用上面的命令就可以解决问

lz需要给出y_0的初值~这个方程满足李普希兹条件,因此,解存在唯一并且可以唯一延拓到边界,应用lax等价定理可以知道,向前欧拉法具有二阶的收敛速度……所以可以尝试用向前欧拉法编写：我刚学matlab

先编一个M函数文件 function x=rr(t,u) %下面的参数值自己设置合适的,否则不一定有解,画不出图来. b1=2; d1=-3;&nbs

clearall;clc;load15.mat;t=celiang(:,1);u1=celiang(:,2);symsu2i5iixu5%定义u1,L1,L2,R1,C1,C3,ii等vU1=celi

好像没有.你可以转写成y',y''型的再用NSolve.当然,也可以用函数D表示导数（包括偏导数）,我想dy,dx形式的转写成用函数D表示要比转写成y',y''型方便.函数D可以在帮助查到.

你之前错将x[0]==x0写成过x[0]=x0吧.这导致了变量赋值.先执行一下Clear[x]再执行你的代码.

求系数矩阵的特征值,特征向量；特征向量求出后,构造基解矩阵,就ok.看书,照猫画虎做一遍,就会了.再问：能给出具体过程吗？再问：这个我真不会

把第三个式子也微分.>>[x,y,z]=dsolve('Dx=(-2)*A*x+u*y','Dy=2*A*x+2*u*z-(A+u)*y','Dx+Dy+Dz=0','x(0)=1,y(0)=0,z(

楼主是想得到含字母形式的解吧?如果我理解的不错,其实有MN个这样的微分方程组,因为j=1,2,…,M,i=1,2,…,N,而且y是关于t的函数,对否?当然了,因为是符号解,其实只要求出一组这样的方程组

你的常数让人不知所措,好歹给组数据吧数据自己编的a0=1;a1=1;a2=1;a3=1;a4=1;a5=1;n=1;key=dsolve('Dx=a0*a1*s*x-a2*x','n*Ds=-n*a0