How to calculate numerically a volume integral over a complicated 3D shape in Python? [closed]