Así pues, se quiere parametrizar el plano de forma que un punto del mismo esté definido por dos parámetros $$ \vec{r} = \mathbf{ S}(u,v) $$
y luego establecer límites a $u$ y $v$ para definir sus límites. En el caso de un plano, se necesita un sistema de coordenadas local con dos direcciones mutuamente ortogonales $\mathbf{e}_1$ y $\mathbf{e}_2$ que también son ortogonales a la normal del plano $ \mathbf{n} \cdot \mathbf{e}_i = 0$ . El centro de origen está en $\vec{r}_0$ .
$$ \vec{r} = \vec{r}_0+ \mathbf{e}_1 u + \mathbf{e}_2 v $$
Ahora la pregunta es: cómo establecer un sistema de coordenadas local en un plano definido por la ecuación $ax +by + cz+d =0 $
La normal del plano es $$\mathbf{n} = \frac{1}{\sqrt{a^2+b^2+c^2}} \pmatrix{a\\b\\c}$$ y el punto del plano más cercano al origen es
$$ \vec{r}_0 = \pmatrix{ \frac{-a d}{a^2+b^2+c^2} \\ \frac{-b d}{a^2+b^2+c^2} \\ \frac{-c d}{a^2+b^2+c^2} } $$
La elección de las direcciones del vector unitario local es algo arbitraria, pero tiene sentido utilizar algo relacionado con los límites que se quieren definir. El único requisito para $\mathbf{e}_1 = \pmatrix{e_x & e_y & e_z}$ es que tiene $e_x a + e_y b + e_z c = 0$
Una forma de conseguirlo es encontrar la dirección más cercana al x -eje
$$\mathbf{e}_1 = \pmatrix{ \frac{b^2+c^2}{a^2+b^2+c^2} \\ \frac{-a b}{a^2+b^2+c^2} \\ \frac{-a c}{a^2+b^2+c^2} }$$
o el y -eje
$$\mathbf{e}_1 = \pmatrix{ \frac{-a b}{a^2+b^2+c^2} \\ \frac{a^2+c^2}{a^2+b^2+c^2} \\ \frac{-b c}{a^2+b^2+c^2} }$$
o el z -eje
$$\mathbf{e}_1 = \pmatrix{ \frac{-a c}{a^2+b^2+c^2} \\ \frac{-b c}{a^2+b^2+c^2} \\ \frac{a^2+b^2}{a^2+b^2+c^2} }$$
y por último, definir la segunda dirección a partir del producto cruzado de $\mathbf{n}$ y la primera dirección
$$ \mathbf{e}_2 = \text{unitvector}( \mathbf{n} \times \mathbf{e}_1 ) $$
Resumen
Con el esquema anterior se puede pasar de $u$ y $v$ a un punto del plano $\vec{r}=\mathbf{S}(u,v)$ . Sus límites se establecen en términos de $(u,v)$ como coordenadas planas. Por ejemplo $u=1$ producirá una línea en el plano que está desplazada del origen del plano en una unidad a lo largo de la línea local v dirección.