function [] = prog1main()
ss=[1,10,30,60];
for i=1:4
 s=ss(i);    
 x(1,1)=0;    
 x(1,2)=24; 
  while x(1,2)<37
    x(1,1)=x(1,1)-s; 
    x(1,2)=x(1,2)-s*f(x(1,2));
  end  
 disp(x)  
end
  
 function Fy = f(y) 
  Fy=-0.0076355*(y-20);
