code: def approx_fhess_p(x0, p, fprime, epsilon, *args):
      f2 = fprime(*((x0 + epsilon * p,) + args))
      f1 = fprime(*((x0,) + args))
      return (f2 - f1) / epsilon
  
imports: from __future__ import division, print_function, absolute_import
  import warnings
  import sys
  import numpy
  from scipy._lib.six import callable, xrange
  from numpy import atleast_1d, eye, mgrid, argmin, zeros, shape, squeeze, vectorize, asarray, sqrt, Inf, asfarray, isinf
  import numpy
  from .linesearch import line_search_wolfe1, line_search_wolfe2, line_search_wolfe2, LineSearchWarning
  from scipy._lib._util import getargspec_no_self
label_desc: Use image node_id 14 for calling a global function with example usage: approx_fhess_p(x0, p, fprime, epsilon) and returns: unknown