Write a double do loop for evaluating the Laplacian using the finite-difference approximation
Start from the skeleton laplacian.F90 where the input array is initialized to u(x,y)=x2+y2.
The analytic solution for that function is
for checking the correctness of your implementation.