function [yvec,tvec] = Adam_Bash(fun,a,b,h,y0)
% Task : 2-step Adams-Bashforth method to solve an ODE


t0 = a;
yprime0 = feval(fun, [t0; y0]);

y1 = y0 + h*yprime0;
t1 = a + h;
yprime1 = feval(fun, [t1,y1]);

tvec = [t0 t1];
yvec = [y0 y1];

t0 = t1;
t1 = t1 + h;

while t1 < b + h/2
	y0 = y1;
	y1 = y0 + (3/2)*h*yprime1 - (1/2)*h*yprime0;
	
	tvec = [tvec t1];
	yvec = [yvec y1];
	
	yprime0 = yprime1;
	yprime1 = feval(fun, [t1; y1]);
	
	t0 = t1;
	t1 = t1 + h;

end


return;