The following are transcripts of sessions identical to Problems 2.1 a and b, but the initial interior points are set to the average of the boundary.
For mesh spacing 0.1:
kestrel1:~/CPS615.dir/hw4> f90 -o laplace laplace.f90 kestrel1:~/CPS615.dir/hw4> laplace MESH RESOLUTION INFO: dx is 0.100000000000000 dy is 0.100000000000000 PROGRAM COMPLETE! total number of Jacobi iteration steps: 204 infinity norm between u and uold for last Jacobi step is 2.900397369831609E-006 infinity norm between u and the analytical solution is 1.618677728844276E-003 kestrel1:~/CPS615.dir/hw4>
For mesh spacing 0.05:
kestrel1:~/CPS615.dir/hw4> f90 -o laplace laplace.f90 kestrel1:~/CPS615.dir/hw4> laplace MESH RESOLUTION INFO: dx is 5.000000000000000E-002 dy is 5.000000000000000E-002 PROGRAM COMPLETE! total number of Jacobi iteration steps: 683 infinity norm between u and uold for last Jacobi step is 2.981035881255423E-006 infinity norm between u and the analytical solution is 4.456010982192282E-004 kestrel1:~/CPS615.dir/hw4>