Add FE convergence

The RHS is:

3 * (cos(x-1) * cos(y-1) * cos(z-1)) / pi
Plot FE convergence
Also calculate the real solution