function [yvec,tvec] = generate_fun(fun,a,b,h)
% Task : plots the function fun


t = a;
y = feval(fun,t);
tvec = [a];
yvec = [y];

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


return;