Question

求用MATLAB解下列微分方程组初值

请把详细的程序写出来，跪谢大神！

Answer 1

改进欧拉
1.函数
function[x,y]=eulerpro(fun,x0,xfinal,y0,n)
if nargin<5
n=50;
end
h=(xfinal-x0)/n;%步
x(1)=x0;y(1)=y0;
for i=1:n
x(i+1)=x(i)+h;
y1=y(i)+h*feval(fun,x(i),y(i));
y2=y(i)+h*feval(fun,x(i+1),y1);
y(i+1)=(y1+y2)/2;
end
end
2.函数
function f=doty(x,y)
f=cos(x*y);
end
3.主函数调用
[x,y]=eulerpro('doty',0,1,1,10)

追问

大神，你这个太高端了，我初学，你能帮我看看图片里哪里不对吗?

就按图片上的先定义函数，然后再调用，可老是出各种错，能帮忙按照图片上的思路写一下吗？