Poisson's equation

An elliptic PDE on an irregular boundary
u(x,y) has Dirichlet, Neuman or mixed boundary conditions
Arbitrary 2D regions (possibly including holes and cracks)
Leads to very big, but sparse linear system
Compute the sparse system in Mathematica