Added do.gradient switch which is currently off by default as the
gradients are fragile and affecting hunt for minimum
Using `monotone' instead of linearize, no longer retain full basis
outside xmin/xmax, but setting monotone removes two columns from full
basis second from left and right) that asymptote at the edge leading
to non-monotonic functions
Cleaning up verbose output
Added verbose option, passing of method to optim, compare analytical
and numeric derivatives with compareDerivative() function in maxLik.
Set fnscale=-1 in optim to maximize function, changed gradient,
sum.log etc. to match
Removed code for lbound.pd/ubound.nd,deriv being fed to optimizer etc.
Gradients properly implemented thanks to compareDeriv in maxLik I was
able to categorically confirm! Wow, the nail in the coffin...
Removed par.lb and lower from code