向后欧拉的MATLAB算法

比方说
dy/dx=-30y;y(0)=1;
用向后欧拉法,在MATLAB中来求解
令h=0.05,x=[0:h:1];
再用迭代法求y(k+1)时,老是陷入死循环,小妹初试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
温馨提示:答案为网友推荐,仅供参考