Using cubic interpolation, a new step size is computed for minimizing a function value.
stepit(linemat, ips, dblwrd, MAXSTEP)
A positive step size.
A 3 by 5 matrix containing step sizes, function values and slopes at multiple values.
If 1, previous slope was positive.
Vector of length 2: dblwrd[1] TRUE means step halved dblwrd[2] TRUE means step doubled.
maximum allowed size of a step.