Discretización de PDEs, matrices por bloques, simetrizaciones y rangos.

Suponemos $latex Delta u = f$ en $latex 2D$, es decir, $latex frac{partial^2}{partial x^2}u(x,y) + frac{partial^2}{partial y^2}u(x,y) = f(x,y)$. Miraremos como queda la matriz del sistema al discretizar, como simetrizarla y su rango en tres casos: condición Neuman respecto $latex x$ en una frontera, condición Neumann respecto $latex y$ en una frontera y condición Neumann …