老师布置了一个coursework, 是用matlab求这个方程组的数值解:dy/dx=cos(t^2)+x;dx/dt=sin(t^2)+y(初始条件是x(0)=0:y(0)=0). 老师规定要用导数定义,下面是我的代码,运行不成功,请教各位大神帮忙解答下;x(1)=0;y(1)=0;t=0:10^-4:10;for n=2:length(t), y(n)=y(n-1)+(cos(t(n-1).^2)+x(n-1))*(10^-4), x(n)=x(n-1)+(sin(t(n-1).^2)+y(n-1))*(10^-4),endplot(t,x)hold onplot(t,y)
一周热门 更多>