code: def xyz2rgb(x__, y__, z__):
      
      x2_ = x__ / 100.0
      y2_ = y__ / 100.0
      z2_ = z__ / 100.0
  
      r__ = x2_ *  3.2406 + y2_ * -1.5372 + z2_ * -0.4986
      g__ = x2_ * -0.9689 + y2_ *  1.8758 + z2_ *  0.0415
      b__ = x2_ *  0.0557 + y2_ * -0.2040 + z2_ *  1.0570
  
      def finv(arr):
          """Inverse
          """
          return np.where(arr > 0.0031308,
                          1.055 * (arr ** (1.0 / 2.4)) - 0.055,
                          12.92 * arr)
  
      return finv(r__) * 255, finv(g__) * 255, finv(b__) * 255
docs: XYZ colorspace to RGB