Numerical solution of PDEs with FreeFem++