У меня есть точка, заданная в 3D, и я хочу выбрать точки, которые находятся внутри трехмерной фигуры, охватывающей границу, заданную как
ro(phi) = sqrt(cos(2*phi)+sqrt(1.4-(sin(2*phi))**2))
X = ro(phi) * cos(phi)
Y = ro(phi) * sin(theta) * sin(phi)
Z = ro(phi) * cos(theta) * sin(phi)
0 <= theta <= 2pi
0 <= phi <= pi

Есть ли способ сделать это в среде Python?
Буду признателен за любую помощь.