Here's the electric potential and field overlaid for a uniform cylindrical charge distribution created using FEniCS, an opensource finite element solver. Although this is a simple analytical problem, this technique can be used for much more complicated geometries.

Next up: linear dielectrics, I've tested in both 2D and 3D and I'm ready to write it up =] surprisingly easy!