code: def inv_transform(xs, ys, zs, invM):
      """
      Transform the points by the inverse of the projection matrix, *invM*.
      """
      vec = _vec_pad_ones(xs, ys, zs)
      vecr = np.dot(invM, vec)
      if vecr.shape == (4,):
          vecr = vecr.reshape((4, 1))
      for i in range(vecr.shape[1]):
          if vecr[3][i] != 0:
              vecr[:, i] = vecr[:, i] / vecr[3][i]
      return vecr[0], vecr[1], vecr[2]
  
imports: import numpy
  from matplotlib import _api
label_desc: Use image node_id 13 for calling a global function with example usage: inv_transform(xs, ys, zs, invM) and returns: vecr,  vecr,  vecr