链接: 参考链接 l = 3;m = 0;t = 0 : pi / 100 : 2 * pi;le = legendre(l, cos(t));R = le(m + 1, :) .^ 2;p = (0 : 50)' / 25 * pi;x = ones(size(p)) * (R .* cos(t));y = R .* sin(t);z = cos(p) * y;y = sin