DIMENSION A(8), B(8) S=8 V=8 N=11 M=12 AN1=1 AM1=1 AN2=1 AM2=1 CSBEG DO 21 I=1,N AN1=AN1*I 21 CONTINUE DO 22 I=1,M AM1=AM1*I 22 CONTINUE A1=AN1/AM1 CSEND CPBEG DO 6 I=1,N,V Z=N-I+1 IF (Z.GT.V) GO TO 2 GO TO 4 2 Z=V 4 DO 1 J=1,Z A(J)=I+J-1 1 CONTINUE DO 5 J=Z+1,V A(J)=1 5 CONTINUE #CALL EMULVV(A,A,CSHIFT(A,4),4) #CALL EMULVV(A,A,CSHIFT(A,2),2) AN2=AN2*A(1)*A(2) 6 CONTINUE DO 16 I=1,M,V Z=M-I+1 IF (Z.GT.V) GO TO 12 GO TO 14 12 Z=V 14 DO 11 J=1,N A(J)=I+J-1 11 CONTINUE DO 15 J=Z+1,V A(J)=1 15 CONTINUE #CALL EMULVV(A,A,CSHIFT(A,4),4) #CALL EMULVV(A,A,CSHIFT(A,2),2) AM2=AM2*A(1)*A(2) 16 CONTINUE A2=AN2/AM2 CPEND IFLAG=0 DO 31 I=1,N IF (A1.NE.A2) IFLAG=1 31 CONTINUE 32 FORMAT(' A1=',F9.2,' A2=',F9.2,' S=',I2,' IFLAG=',I1) WRITE(6,32) A1,A2,S,IFLAG STOP END |