Quote (TempoONE @ 29 Jan 2013 03:49)
I compiled it and proceeded to see if it works, but it's flawed. Where's your declaration statements? and sub routines?
those above are the SUBROUTINES, the main program follows (putting "." in the first 6 columns instead of space because of space suppresion):
......IMPLICIT INTEGER (A-Z)
......EXTERNAL FUNCTION DELTA,EPSILON
C
C only using integers in this program
C
......DO 1050 K=1,3
......DO 1040 I=1,3
......DO 1030 J=1,3
......DO 1020 P=1,3
......DO 1010 Q=1,3
......RIGHT=DELTA(I,P)*DELTA(J,Q)-DELTA(I,Q)*DELTA(J,P)
......LEFT=EPSILON(I,J,K)*EPSILON(P,Q,K)
......IF(LEFT.EQ.RIGHT) GOTO 1010
C
C we found a combination the formula does not work
C
......WRITE(6,1000) I,J,P,Q,K
.1000.FORMAT(' the identity is not true for I=',I1,' J=',I1,' P=',I1,' Q=',I1,' and K=',I1)
......GOTO 9000
C
C continue to check the other values
C
.1010.CONTINUE
.1020.CONTINUE
.1030.CONTINUE
.1040.CONTINUE
.1050.CONTINUE
C
C we tested all combinations and they worked
C
......WRITE(6,2000)
.2000.FORMAT(' the identity xxx is correct')
C
.9000.STOP
......END
remember that this is FORTRAN V and you have to sort out the compilation/binding depending on your operating system
PS: not column 6 but column 7 the statements start, column 6 is for statement continuation - sorry about saying 6 in previous postThis post was edited by brmv on Jan 28 2013 10:17pm