I have below code

import scipy as scp    
def fun(param, Power, SomLog):
    return (param[0] ** Power[0] + param[1] ** Power[1] + param[2] ** Power[2]) if SomLog else (param[0] ** Power[0] + param[1] ** Power[1])

scp.optimize.minimize(fun, [0.5, 0.5, 0.5], args = ([2, 3, 3], True)).x

However I also want to get the estimate of Hessian matrix at the point of optimized parameters. Could you please point if there is any method to obtain the same.