3616.sagews

def spherical_to_cartesian(r,azimuth,inclination):
    return (r*cos(azimuth)*sin(inclination), r*sin(azimuth)*sin(inclination),r*cos(inclination))

data=[(1,0,0),(1,pi/2,pi/2),(1,pi,pi/3),(2,pi,pi/5),(1.5,3*pi/2,2*pi/3)]

points([spherical_to_cartesian(*p) for p in data],size=30)