code: def fun1(beta, y, x):
      # print(beta.shape, x.shape)
      xb = np.dot(x, beta)
      return (y - xb) ** 2
  
imports: import numpy
  from numpy.testing import assert_allclose, assert_almost_equal
  import statsmodels.api
  from statsmodels.tools import numdiff
  from statsmodels.tools.numdiff import approx_fprime, approx_fprime_cs, approx_hess_cs, _approx_fprime_scalar, _approx_fprime_cs_scalar
label_desc: Use image node_id 3 for calling a global function with example usage: fun1(beta, y, x) and returns: unknown