def find_hessian(point, vars=None, model=None):
 model = modelcontext(model)
 H = model.fastd2logp(vars)
 return H(Point(point, model=model))

def find_hessian ( point , vars = None , model = None ) : model = modelcontext ( model ) H = model . fastd2logp ( vars ) return H ( Point ( point , model = model ) )
returns hessian of logp at the point passed .