FEM to solve for 2D Laplace's Equation

FEM solution to a steady-state PDE (example: Laplace's) applied to an electric potential

%% Laplace's Equation: nabla^2 u = 0 (version 2: acquire matrix results)
% 2 space dimensions: uxx + uyy = 0, where u = V (electric potential)
% by Prof. Roche C. de Guzman
clear; clc; close('all');

%% Given
% source geometry (circle)
r = 0.5; % radius [m]
C = [0 0]; % center (h,k)
% domain geometry (rectangle)
f = 1.5; % factor relative to the circle
% electric potential
Vg = 0; % ground
Vs = 110; % source [V]