Kugelkoordinaten passen da irgendwie nicht, dafür müsste sowas wiedastehen. Ich habe mal rumprobiert mit beispielsweise:
und komme dann auf das Integral:
Das Hauptproblem liegt dann halt darin, dM zu bestimmen, hierzu müssen die Differentiale dx, dy und dz als Funktion von r, t und s bestimmt werden, da die Transformation allerdings nicht linear ist, ist es ungemein schwierig, sie umzukehren. Eventuell geht es, wenn die Diffenentiale dr, dt und ds aus der Transformation als Funktion von x, y und z berechnet werden und anschliessend irgendwie umgeformt werden können. Doch bei mehrdimensionalen Differentialen bin ich mir nicht mehr sicher, wie das genau geht, da noch partielle Differentiale hinzukommen könnten.