functions in romberg.i - m

             max_doublings= 20  
 
     is the maximum number of times romberg or simpson will split the  
     integration interval in half.  Default is 20.