Download Programmes Fortran pour distances de rangs, constellations et

Transcript
16
1002
120
600
700
601
701
602
702
800
801
1000
2000
WRITE(IWRITI~~O)VIFICH
PA5DANS
F O R M A T ( ~ H ~ Y ~ ~ ANE
~ , SE
/ TROUVE
@
GOTO 1 0 0 0
WRITE(IWRIT,700)
FORMAT(lH1,'JE
HORS DES L I M I T E S ' )
GOTO 1 0 0 0
WRITE(IWRIT,701)
FORMAT(lH1,'IVM
HORS D E SL I M I T E S ' )
GOTO 1 0 0 0
WRITE(IWRITI~O~)
HORS D E SL I M I T E S ' )
FORMAT(lH1,'IV
GOTO 1 0 0 0
WRITE(I~RIT~BO1)
F O R M A T ( l H l r ' P L U . 5 DE 1 0 CARTESFORMAT')
STOP
END
SUBROUTINELECT(IVM,JEM*UIVD,B)
COMMON/C1/IREAD~IdRIT,IPUIKFOR
I M P L I C I TI N T E G E R ( V )
DIMENSION U(1) 'VD(1)
D I M E N S I O N 7(400)r F M T ( 1 3 0 )
KFOR=KFOR*13
READ(IREAD*2000) (FMT(I),I=lrKFOR)
FORMAT (13A6)
L=O
DO SOO J = l . J E M
READ(IREADtFMT,END=502) (Z(I)rI=l~IVM)
DO 5 0 1 I = l r I V M
501
I F ( V D ( 1 ) mNE.1)GOTO
L=L+1
501
500
502
1
U(L)=Z(I)
CONTINUE
CONTINUE
RETURN
WRITE(IWRIT*l)
F O R M A T ( @ I L MANQUE
RETURN 5
END
DES CARTESDONNEESO
SURROUTINERGEQU(XJ)
COMMON/C4/JE,I,JRE,IV
COMMOPJ X ( 3 0 0 0 ) r X B I S ( 4 2 0 0 0 )
D I M E N S I O NA ( 1 0 0 ) q R ( 1 0 0 )
XRIS(l)=O.
T=O.
DO 1 J = l r J E
H(J)=O.
DO 2 J = l r J E
1
Jl=J-1
2
3
A(J)=X(Jl*IV+I)
MP= J E
DO 13 J=l * J E
IF(R(J))3*3*13
IF(R(J)+1.)4r13*4
4
s=o.
6
E=O e
30
31
Y=A ( J I
IF(Y+1.)30912,30
DO 9 K = l , J E
IF(A(K)+1.)31,9131
If ( A ( K ) - Y ) 7 r 8 r 9
LE F I C H I E R ' r 8 A l )