from sympy import symbols, sin, pi, Ynm, re, lambdify
from spb import *
theta, phi = symbols('theta phi')
r = re(Ynm(3, 3, theta, phi).expand(func=True).rewrite(sin).expand())
graphics(
    surface_spherical(
        abs(r), (theta, 0, pi), (phi, 0, 2 * pi), "real",
        use_cm=True, n2=200,
        color_func=lambdify([theta, phi], r)),
    backend=PB)