from sympy import symbols, sin, pi
from spb import *
theta, phi = symbols('theta phi')
r1 = 1
r2 = 1.5 + sin(5 * phi) * sin(10 * theta) / 10
graphics(
    surface_spherical(r1, (theta, 0, pi / 2), (phi, 0.35 * pi, 2 * pi),
        label="r1", wireframe=True, wf_n2=25),
    surface_spherical(r2, (theta, 0, pi / 2), (phi, 0.35 * pi, 2 * pi),
        label="r2", wireframe=True, wf_n2=25),
    backend=PB)