code: def fun(self, coords):
      x, y = coords
      p1 = (1.5 - x + x * y) ** 2
      p2 = (2.25 - x + x * y**2) ** 2
      p3 = (2.625 - x + x * y**3) ** 2
      return p1 + p2 + p3
  
imports: import time
  import numpy
  from numpy import sin, cos, pi, exp, sqrt, abs
  from scipy.optimize import rosen
label_desc: Use image node_id 1 for calling the Beale obj's underlying member method code with example usage: obj.fun(coords) and returns: unknown