Finite Element Method

  The two-dimensional Laplace equation around a cylinder is given. 1- Implement the fully implicit Finite 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.