code: def Thurber(b, x, y=0):
      b = read_params(b)
      model = (b[0] + b[1] * x + b[2] * x**2 + b[3] * x**3) / (
          1 + b[4] * x + b[5] * x**2 + b[6] * x**3
      )
      return model - y
  
imports: import os
  from numpy import arctan, array, cos, exp, log, sin
  from lmfit import Parameters
label_desc: Use image node_id 23 for calling a global function with example usage: Thurber(b, x, y) and returns: unknown