PovRay scene: orthocircle.pov
[(x2 + y2 - 1)2 + z2] [(y2 + z2 - 1)2 + x2] [(z2 + x2 - 1)2 + y2] - c12 [1 + c2 (x2 + y2 + z2)] = 0 Where c1 = 0.075 and c2 = 3