3 球面調和関数の絶対値を球面上の濃淡でプロット

plotSphericalHarmonic[l_, m_] := Block[{saidai},
    saidai = findmaxSphericalharmonic[l, m];
    ParametricPlot3D[
      {x[theta, phi], 
        y[theta, phi],
        z[theta, phi],
        {
          EdgeForm[],
          GrayLevel[     
            (Abs[SphericalHarmonicY[l, m, theta, phi]]/saidai)^2
            ]
          } 
        },
      {phi, 0, 2Pi}, {theta, 0, Pi},
      Lighting -> False,
      Boxed -> False,
      Axes -> None,
      PlotPoints -> 100,
      PlotLabel -> SequenceForm["l=", l, ",m=", m]
      ]]