Measuring the Hessian matrix for lsqnonlin

Hello

I am optimising a multiobjective function using the lsqnonlin function with the trust-region reflective algorithm.

I was wondering how I can get the hessian matrix as an output?

Is it possible to calculate it using the Jacobian matrix at the optimum solution that lsnqnonlin provides?

Thanks