欧拉法的matlab程序

欧拉法的matlab程序

f unction oulafa(f,a,b,x0,y0,n,g)

f=inline(f);

g=inline(g);

h=(b-a)/n;

y(1)=y0;

x=a:h:b;

disp('k x(k) y(k) y(x) e ') for i=1:n

y(i+1)=y(i)+h*f(x(i),y(i));

gg(i+1)=g(x(i+1));

e(i+1)=abs(g(i+1)-y(i+1));

fprintf('%d %f %f %f %f\n',i,x(i+1),y(i+1),gg(i+1),e(i+1))

end