Adding checking for positive/negative derivatives at the left/right
bound (for monotonically decreasing/increasing at the left/right
boundary one would want these to be FALSE which is the default)
Propogated change from norm to er (naming convention) through plot
Fixed glitch in testing for lbound.pd ubound.nd (equal to zero fine),
fixed plot when ylim passed as an arg