code: function totalRayleigh(lambda) {
      return [
          (8.0 * pow(pi, 3.0) * pow(pow(n, 2.0) - 1.0, 2.0) * (6.0 + 3.0 * pn)) / (3.0 * N * pow(lambda[0], 4.0) * (6.0 - 7.0 * pn)),
          (8.0 * pow(pi, 3.0) * pow(pow(n, 2.0) - 1.0, 2.0) * (6.0 + 3.0 * pn)) / (3.0 * N * pow(lambda[1], 4.0) * (6.0 - 7.0 * pn)),
          (8.0 * pow(pi, 3.0) * pow(pow(n, 2.0) - 1.0, 2.0) * (6.0 + 3.0 * pn)) / (3.0 * N * pow(lambda[2], 4.0) * (6.0 - 7.0 * pn))
      ]
  }
docs: returns vec3