Thanks for you post. I found it very interesting.

However, I think it is very important to obtain estimates of standard errors of parameters using the inverse of hessian matrix.

I think you could extend this example using the CALL NLPFDD