2 原点からの距離が球面調和関数の絶対値の二乗を表すプロット

lengthplotSphericalHarmonic[l_, m_] := ParametricPlot3D[
    {Abs[ SphericalHarmonicY[l, m, theta, phi] ]^2 x[theta, phi],
      Abs[ SphericalHarmonicY[l, m, theta, phi] ]^2 y[theta, phi],
      Abs[ SphericalHarmonicY[l, m, theta, phi] ]^2 z[theta, phi],
      EdgeForm[]
      }
    ,
    {theta, 0, Pi},
    {phi, 0, 2 Pi}
    ,
    PlotPoints -> 100,
    Boxed -> False,
    Axes -> None,
    PlotLabel -> SequenceForm["l=", l, ",m=", m]
    ]