function [tt,yy,steps,res] = Driver(IC,n,h,TMax,tolerance)
 
  y = IC;
  tt=0;
  yy=[];
  steps=[];
  res = [];
  
   while tt <= TMax
      [YNew,counter,residual] = ODESOL(y,h,tolerance);
      yy = [yy;YNew];
      tt = tt + h;
      y = YNew;
      steps = [steps;counter];
      res = [res;residual];     
  end
  
  
 
      

      
 
 


    