Added maxfun keyword to Optimizers for scipy.optimize routines
to be able to make larger number of evaluationg for the
negative log likelihood function.