But something is very weird with the function values.
Without constraints:
Quasi-Newton BFGS minimisation
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Line search: Backtracking line search.
k: 0 xk: [ 11.8464633, 5.383340723, 10.37739135] fk: 22.1840340996
Parameter values: [11.744468102051812, 5.3843668638102162, 9.8277882468793187]
Function value: 19.089946986005284
Iterations: 8
Function calls: 57
Gradient calls: 9
Hessian calls: 0
Warning: None
With constraints:
Logarithmic barrier function
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
k: 0 xk: [ 11.8464633, 5.383340723, 10.37739135] fk: 73094199769.8
Entering sub-algorithm.
Quasi-Newton BFGS minimisation
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Line search: Backtracking line search.
k: 0 xk: [ 11.8464633, 5.383340723, 10.37739135] fk: 73094199769.8
Parameter values: [11.846463301422499, 5.3833407234324602, 10.377391345839921]
Function value: 73094199769.763443
Iterations: 1
Function calls: 91
Gradient calls: 2
Hessian calls: 0
Warning: None
task #7824(https://gna.org/task/index.php?7824): Model parameter ERROR estimation from Jacobian and Co-variance matrix of dispersion models.