just take spherical coordinates(x=alpha, y=beta) and convert them via sin and cos to 3d space (vec_rotate((1, 0, 0), (alpha, beta, 0)))
this won't give you a uniform distribution, but a spherical shape