code: def f(x, y, pts):
      z = np.zeros_like(x)
      for p in pts:
          z = z + 1 / (np.sqrt((x - p[0]) ** 2 + (y - p[1]) ** 2))
      return 1 / z
  
imports: import time
  import matplotlib.pyplot
  import numpy
label_desc: Use image node_id 2 for calling a global function with example usage: f(x, y, pts) and returns: unknown