1- Implement the fully implicit Finite element solution algorithm and use a direct solver (LU factorization) to solve.
2- Plot Error function versus deltaR and DeltaTheta in a log-log scale. Compute the spatial convergence rate.