I've edited line 1524 of 3dmain.f to read
     DO 1 I=ICOL,NUMEP1-1
It used to read 
     DO 1 I=ICOL,NUMEP1
which causes the program to run off of the end of the array XMATRIX.

I've changes line 914 to of 3dmain.f print out I5 instead of I3 so that
matricies larger with dimensions larger than 1000 will have their dimensions
printed.
