Numerical Solution of Electrostatics with Dielectrics in 2D