球谐函数

Question

为了制作我的课堂笔记，我需要在 Asymptote 中可视化球谐函数。使用来自https://mathworld.wolfram.com/SphericalHarmonic.html

triple fs(pair p){
real x=sin(p.x)*cos(p.y)*(sqrt(3/(4*pi))*sin(p.x)*cos(p.y))^2;
real y=sin(p.x)*sin(p.y)*(sqrt(3/(4*pi))*sin(p.x)*cos(p.y))^2;
real z=cos(p.x)*(sqrt(3/(4*pi))*sin(p.x)*cos(p.y))^2;
return (x,y,z);
}
surface s=surface(fs,(0,pi),(0,2pi),100);
draw(s);