Мне нужно равномерное распределение точек на 4-мерной сфере. Я знаю, что это не так тривиально, как выбрать 3 угла и использовать полярные координаты.

В 3 измерениях я использую

from random import random

u=random()
costheta = 2*u -1 #for distribution between -1 and 1
theta = acos(costheta)
phi = 2*pi*random

x=costheta
y=sin(theta)*cos(phi)
x=sin(theta)*sin(phi)

Это дает равномерное распределение x, y и z.

Как я могу получить аналогичное распределение для 4 измерений?