If you just need to compute the hessian matrix of a general function (often called a black box: you pass stuff into the black box and stuff comes out, based on some unknown algorithm.)
As long as your function is reasonably well behaved, you can use my derivest toolbox, found on the file exchange .
I'll show it here, using a function handle.
rosen = @(x) (1-x(1)).^2 + 105*(x(2)-x(1).^2).^2;
Hmat = hessian(rosen,[2 3])
Hmat =
3782 -840
-840 210
If you have an m-file, just do this:
hessian(@fob_8_52, [4 6])
