Calculation: PHsIoD - 1

 (x-x0)^2+ (y-y0)^2+ (z-z0)^2  = 
 R^2
 (x-x0)^2+ (y-y0)^2+ (z-z0)^2 =  R^2(x-x0)^2+(y-y0)^2+(z-z0)^2 = R^2
 (x-x0)^2 =  R^2- (y-y0)^2- (z-z0)^2(x-x0)^2 = R^2-(y-y0)^2-(z-z0)^2
 (x-x0)^2 =  R^2- y^2+ 2* y* y0- y0^2- z^2+ 2* z* z0- z0^2(x-x0)^2 = R^2-y^2+2*y*y0-y0^2-z^2+2*z*z0-z0^2
x = saknis( R^2- y^2+ 2* y* y0- y0^2- z^2+ 2* z* z0- z0^2)+x0x = saknis(R^2-y^2+2*y*y0-y0^2-z^2+2*z*z0-z0^2)+x0