4 球面調和関数の実部の値を値を球面上の濃淡でプロット

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