f = inline('2 * y * t','t','y')
   
    t0 = 0.1;
    y0 = 1;
    t1 = 1;
	h = 0.05;
	n = 0;
    
    if (n <=0)
        n = abs(t1 - t0)/h;
    end
    
    Yr =[];
    Tr =[];
    
     y = y0;
	 t = t0;
     
     for j = 1:n
	 k = feval(f,y,t);;
	 y = y + h*k;
	 t = t+h;
	 Yr(j) = y;
	 Tr(j) = t;
	end
        
	R = [Tr ; Yr];
    R'