return column vector (t(i+1) & y(i+1)) instead of row vector
1 次查看(过去 30 天)
显示 更早的评论
function [y,t] = euler_timestep(f,t0,tf,alpha,N)
h=(tf-t0)/N;
t(1)=t0;
y(1)=alpha;
for i=1:N
t(i+1)=t(i)+h;
y(i+1)=y(i)+h*f(t(i),y(i));
end
0 个评论
采纳的回答
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!