x = r * sin(theta) * cos(phi); y = r * sin(theta) * sin(phi); z = r * cos(theta);
d2p(x, y, z, &r, &theta, &phi);
( は狭義単調減少で、 , であるから、 区間 にただ1つの実数解を持つことが分かる。)