The following is a transcript from the session where I solved the Laplace equation using a mesh spacing of 0.05 and an initial internal value of 0.0. I used the kestrels with the DEC F90 compiler:
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: 1267 infinity norm between u and uold for last Jacobi step is 2.998694651701328E-006 infinity norm between u and the analytical solution is 5.103057795947663E-004 kestrel1:~/CPS615.dir/hw4>