code: def lnprob(p):
      resid = residual(p, just_generative=True)
      return -0.5 * np.sum(
          ((resid - y) / dy) ** 2 + np.log(2 * np.pi * dy**2)
      )
  
imports: import matplotlib.pyplot
  import numpy
  import lmfit
label_desc: Use image node_id 4 for calling a global function with example usage: lnprob(p) and returns: unknown