Different finite element formulations for Poisson problems are implemented and compared.

The methods are applied to two model problems.

The first problem has smooth solution and the second one presents a square root singularity.

For the problem with a smooth solution, the results confirm convergence orders.

For the singular problem, hp adapted meshes tuned to the singularity are created.