We start with preprocessing using the bounding box of the RPP surface
patch coupled with some level of subdivision. Then we solve a system
,
using Newton's method and verify the
results with the third equation
. Projection of
onto a planar approximation of the surface (faceting) may
provide a good initial approximation.