球面調和関数 の可視化を行います。量子力学でいうと l は方位量子数、m は磁気量子数になります。用いる式は以下のとおりです。自明ですが、r が定数の場合は球面となります。 l=3,m=1 SageShell を起動し、以下のとおり入力します。 sage: phi,theta = var('phi,theta') sage: Y31 = spherical_harmonic(3,1,theta,phi) sage: r = (Y31(phi,theta).conjugate()*Y31(phi,theta)).real() sage: parametric_plot3d((r(phi,thet…