Correct the handling of bounds
Add PSwarm to MaxLik.m
Add the option Vectorized to numjac.m to speed-up calculations
Prevent the calculation of covariance and hessian when parameters are at their bounds