code: def poly_div(N, D):
      dD = degree(D)
      dN = degree(N)
      if dD < 0:
          raise ZeroDivisionError
      if dN >= dD:
          q = [0] * dN
          while dN >= dD:
              d = [0] * (dN - dD) + D
              mult = q[dN - dD] = N[-1] / float(d[-1])
              d = [coeff * mult for coeff in d]
              N = [
                  fabs(coeffN - coeffd) for coeffN, coeffd in zip(N, d)
              ]
              dN = degree(N)
          r = N
      else:
          q = [0]
          r = N
      return q, r
  
label_desc: Use image node_id 2 for calling a global function with example usage: poly_div(N, D) and returns: q,  r