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
	k1 = f(y,t);
	k2 = f(y+ h*k1, t+h);
        y = y + (h/2)*(k1+k2);
	t = t+h;
	Yr(j) = y;
	Tr(j) = t;
	end
	
	R = [Tr Yr]