What about generating 4 random numbers, each uniformly distributed. Make it a vector and normalize it. Now use it as a quaternion rotation, apply it to to vector (1,0,0)