向後歐拉的MATLAB算法

y(k+1)=y(k)-30*h*y(k+1)

變形求得:

y(k+1)=y(k)/(30*h+1)

故MATLAB程序有:

h=0.05;

x=[0:h:1];

y(1)=1;

for k=1:length(x)-1

y(k+1)=y(k)/(30*h+1);

end

plot(x,y,'r.-');

title('向後歐拉');

grid on