oulaqiancha.m

来自「1．编写欧拉前差、后差、梯形公式。 2．编写二阶、三阶龙格库塔法通用程序。 」· M 代码 · 共 15 行

function [yt,th]=oulaqiancha(f,y0,t0,td,h)
y=y0;
t=t0;
yt(:,1)=double(y0);
th(1)=double(t0);
for i=2:fix((td-t0)/h)+1
fn=double(eval(f));
yt(:,i)=yt(:,i-1)+fn'.*h;
th(i)=th(i-1)+h;
y=yt(:,i);t=th(i);
end
%n=length(fn);
%C=zeros(1,n);C(1)=1;
%out=C*yt;