function [yvec,tvec] = Euler(fun,a,b,h,y0)
% Task : Euler's method to solve an ODE


t = a;
y = y0;
tvec = [a];
yvec = [y0];

while t < b
	yprime = feval(fun, [t; y]);
	t = t + h;
	y = y + h*yprime;
	tvec = [tvec t];
	yvec = [yvec y];
end


return;