home *** CD-ROM | disk | FTP | other *** search
- FUNCTION COT(X)
- REAL*8 X,COT
- COT=1.D+0/TAN(X)
- RETURN
- END
- SUBROUTINE SEARCH(MAXNOY,INP,J,KJ,M)
- INTEGER NY1(3000)
- COMMON /I2BYT/NY1
- NX=M*J+KJ
- NY=M*KJ+J
- INP=0
- DO 100 JJ=1,MAXNOY
- IF (NY1(JJ).EQ.NX.OR.NY1(JJ).EQ.NY)INP=JJ
- 100 CONTINUE
- IF(INP.EQ.0)MAXNOY=MAXNOY+1
- IF(INP.EQ.0)INP=MAXNOY
- RETURN
- END