from sympy import *
from spb import *
phi, theta = symbols("phi, theta")
r = 0.75
params = {
    phi: (-pi/2, -pi, pi),
    theta: (2*pi/3, -pi, pi),
}
graphics(
    arrow_3d((0, 0, 0), (1, 0, 0), rendering_kw={"color": "k"},
        show_in_legend=False),
    arrow_3d((0, 0, 0), (0, 1, 0), rendering_kw={"color": "k"},
        show_in_legend=False),
    arrow_3d((0, 0, 0), (0, 0, 1), rendering_kw={"color": "k"},
        show_in_legend=False),
    arrow_3d(
        (0, 0, 0),
        (r*sin(theta)*cos(phi), r*sin(theta)*sin(phi), r*cos(theta)),
        params=params),
    xlabel="x", ylabel="y", zlabel="z",
    xlim=(-1.5, 1.5), ylim=(-1.5, 1.5), zlim=(-1.5, 1.5), aspect="equal"
)