clear	
	g = '2*t*y';
	function out = f(y,t)
	 out = evstr(g)
	endfunction


	t0 = 0.1;
        y0 = 1;
	h = 0.05;
	n = 20;
Yr =[];
Tr =[];	
        y = y0;
	t = t0;

	for j = 1:n
	k = f(y,t);
	y = y + h*k;
	t = t+h;
	Yr(j) = y;
	Tr(j) = t;
	end
        
	R = [Tr Yr]