Function uses the third-order Runge–Kutta scheme
1 次查看(过去 30 天)
显示 更早的评论

function y = odesRK3(f, t, y0)
n=length(t);
y=nan(length(y0),n);
y(:,1)=y0(:);
for k=1:n-1
h=t(k+1)-t(k);
F1=h*f(t(k),y(:,k));
F2=h*f(t(k)+h/2,y(:,k)+(F1/2));
F3=h*f(t(k)+0.75*h,y(:,k)+(0.75*F1));
y(:,k+1)=y(:,k)+(2*F1+3*F2+4*F3)/9;
end
2 个评论
John D'Errico
2022-10-15
This must be like a quiz show. We are given the answer, and now we need to guess the question.
I will guess, "What color is the bluebird?"
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Communications Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!