The following is a transcript from the session where I solved the Laplace equation using a mesh spacing of 0.10 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 0.100000000000000 dy is 0.100000000000000 PROGRAM COMPLETE! total number of Jacobi iteration steps: 360 infinity norm between u and uold for last Jacobi step is 2.986653990788568E-006 infinity norm between u and the analytical solution is 1.636604399569030E-003 kestrel1:~/CPS615.dir/hw4>