%
% Euler's Method for solving
%   y'=f(t,y), a <= t <= b, y(a)=y0
%
% [tv,yv,n]=eulfun(a,b,y0,h);
% 
% eulfun.m calls M-function fun.m 
% which evaluates f(t,y) for given
% t and y.
%
function [tv,yv,n]=eulfun(a,b,y0,h);
n=floor((b-a)/h);
tv(1)=a;
yv(1)=y0;
for ii=2:n+1;
    fty=fun(tv(ii-1),yv(ii-1));
    yv(ii)=yv(ii-1)+h*fty;
    tv(ii)=tv(ii-1)+h;
end;