求助一道改进欧拉法matlab程序的详细解析,每行都有注释的那种有没有大神帮忙
function [t,x1,x2]=ModifiedEuler(a,b,c,d,n) h=(b-a)/n;
t=a:h:b;
x1=[c zeros(1,n)];
x2=[d zeros(1,n)];
for i=1:n y=ODE(x1(i),x2(i));
yn1=x1(i)+h*y(1);
yn2=x2(i)+h*y(2);
dx=ODE(yn1,yn2);
x1(i+1)=x1(i)+(h/2)*(y(1)+dx(1));
x2(i+1)=x2(i)+(h/2)*(y(2)+dx(2));
end