from sympy import *
from spb import *
x, y, z = symbols("x, y, z")
r = sqrt(x**2 + y**2)
d = symbols('d')
expr = 10 * cos(r) * exp(-r * d)
graphics(
    surface(
        expr, (x, -10, 10), (y, -10, 10), label="z-range",
        params={d: (0.15, 0, 1)}, n=51, use_cm=True,
        wireframe = True, wf_n1=15, wf_n2=15,
        wf_rendering_kw={"line_color": "#003428", "line_width": 0.75}
    ),
    title = "My Title",
    xlabel = "x axis",
    ylabel = "y axis",
    zlabel = "z axis",
    backend = PB
)