Matlab实现——Euler's Method

欢迎前往个人博客 驽马点滴 和视频空间 哔哩哔哩-《挨踢日志》

euler.m

%Program9.1 (Euler's Method)

function E=euler(f,a,b,ya,M)

%Input - f is the function entered as a string ??f??

%         - a and b are the left and right endpoints

%         - ya is the initial condition y(a)

%         - M is the number of steps

%Output - E=[T?? Y??] where T is the vector of abscissas 

%              and Y is the vector of ordinates

h=(b-a)/M;

T=zeros(1,M+1);

Y=zeros(1,M+1);

T=a:h:b;

Y(1)=ya;

for j=1:M

Y(j+1)=Y(j)+h*feval(f,T(j),Y(j));

end

E=[T' Y'];

 

fun.m

 

function f=fun(t,y)

f=3*y+3*t;

 

 

 

 

untitled.m

 

a=0;
b=2;
ya=1;
M=20;
f='fun';
E=euler(f,a,b,ya,M)

欢迎前往个人博客 驽马点滴 和视频空间 哔哩哔哩-《挨踢日志》