You can pick a uniform random orientation without trig functions by first generating a random point in the unit disk via rejection sampling, then projecting it onto the boundary [0].

Of course, using rejection sampling for disk points will give you an estimate for π more directly.

[0] https://mathworld.wolfram.com/CirclePointPicking.html