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 ordinatesh=(b-a)/M;T=zeros(1,M+1);Y=zeros(1,M+1);T=a:h:b;Y(1)=ya;for j=1:MY(j+1)=Y(j)+h*feval(f,T(j),Y(j));endE=[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)