Replacing calls to E04YCF

This routine is a follow-up to E04FCF, E04FDF and E04GBF to convert the Jacobian of a least-squares minimisation to the covariance matrix of the fitted parameters. See the documentation of E04YCF, of the minimisation routine used previously, and of PDA_LMDIF. Migration hints are not yet available.