[Bug-gsl] ode-initval/bsimp is always high order

[Tom Quinn]

The bsimp_alloc() routine uses a very small EPSILON to set the order of 
the integrator, and there seems to be no interface to change this.

The result is that even for small steps and large error tolerances, a 
minimum of 140 function evaluations are taken per timestep.  For my 
application, this is more than an order of magnitude more than are 
necessary.

Tom Quinn       Astronomy, University of Washington
Internet:       address@hidden
Phone:          206-685-9009