用MATLAB解微分方程y^3-y 2t y=0,y(0)=1,0

%也不知道是你抄错了,还是这本书太烂,如果是后者,赶紧换一本.%把下面的保存到当前目录,文件名为cir1.mfunctiondy=cir1(t,y)Vi=20;R=5;C=70e-6;dy=zeros

%给你参考一下我当年做的一个类似的程序%cleardt=0.004;t=0:dt:5;a=-pi:0.2:pi;l=100;r=8;d=10;%导弹尺度参数：长度、半径、火焰尺度L=l;D=l/3;c

function[x,y]=Classical_RK4(odefun,xspan,y0,h,varargin)x=xspan(1):h:xspan(2);y(1)=y0;fork=1:length(x

一维导热问题,应该采用有限差分法求解.具体方法随便找个数值计算方法书都有

设缉私艇航速为v1,走私船航速为v2,航行时间为t,以缉私艇出发点为原点,正东方向为x轴,正北方向为y轴,缉私艇经t时间行驶到（x,y）点.经过的航迹为y=f（x）；y'=(v2*t-y）/(c-x)

请查看dsolve命令,输入helpdsolve回车即可!

这次看看,程序通了.functionhhh[t,x]=ode45(@xprim2,[0,20],[30;20]);plot(t,x);xlabel('timet0=0,tt=20');ylabel('

恩你没有给初值不要解答呀记住加分哦%bydynamic%2009.2.16%%我这里假设初值全部为1%下面是解析解root=dsolve('D2x=2*x*y','D2y=y^2-x^2-20^2',

前两天刚回答了一个类似的问题：再问：非常感谢你，那个问题也是我发的，但我处理二阶就不会了，失败了。听同学说没有解析解，想问下怎么写二阶形式的求解，只要ODE45的方法就行了，别的不麻烦您了再答：>>o

dy(1)是对y(1)求导数的意思后面一样了

A=[-53;100-301];[VD]=eig(A);lambda=diag(D);C=V\[50;100];t=linspace(0,1,1000)';y1=C(1)*V(1,1)*exp

这只是求偏导数而已,你别搞笑哈.给你代码symsxyzz=(x^2+y^2)*exp((x^2+y^2)/x/y);z_x=diff(z,x)z_xx=diff(z,x,2)z_xy=diff(dif

参考代码：a=0.5; b=0.6;dIdt=@(t,I)a*(1-I)*I-b*I;ode45(dIdt,[0 10],0)存在的问题是,如果初值为0,可以计算出导数dI/dt=

这样可以吗?t=1:6m=[556569727885]qi=[151418192120]ppi=[212434]v=[202225273133]plot(t,[m;v],'o-')vsp=spline

    这个方程是解不出解析解的.    我用ODE45给您算了一下.    先

先定义地字符,再求eg：symsxy;dsolve(Dy=x)simiplify(y)或直接dsolve('Dy=x');simiplify(y)

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

functionf=ffun(t,x)f=zeros(3,1);f(1)=-x(1)+0.5*(x(3)+x(2))*(x(3)+x(2)+sin(x(3)));f(2)=-0.5*x(1)*(x(3

如图2,改变gain=100 为 gain=4.改变gain= 1000为 gain= - 3

functionhahaha[t,r]=ode45(@myfun,[01],0)functiondrdt=myfun(r,t)drdt=-i*(50.24-235.5*t)*r+3.14*(1-r^2