      SUBROUTINE INITAL
C...  
C..,  BOUNDARY VALUE ODE INTEGRATED AS AN INITIAL VALUE ODE
C...
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      COMMON/T/S,NFIN,NORUN/Y/Y1,Y2/F/DY1DS,DY2DS/PARM/A
C...
C...  PARAMETER A
      A=11.417542D0
C...
C...  INITIAL CONDITIONS
      Y1= 0.0D0
      Y2=-0.873849D0
      RETURN
      END
      SUBROUTINE DERV
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      COMMON/T/S,NFIN,NORUN/Y/Y1,Y2/F/DY1DS,DY2DS/PARM/A
C...
C...  ODES
      DY2DS=-A*((1.0D0-Y2**2)**0.5D0)*Y1
      DY1DS=Y2
      RETURN
      END
      SUBROUTINE PRINT(NI,NO)
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      COMMON/T/S,NFIN,NORUN/Y/Y1,Y2/F/DY1DS,DY2DS/PARM/A
C...
C...  TEST FOR THE END OF THE RUN
      IF(S.GT.0.55D0)RETURN
C...
C...  MONITOR THE RUN
      WRITE(*,*)NORUN,A,Y1,Y2
C...
C...  AT THE END OF THE RUN, PRINT Y1 TO DETERMINE HOW CLOSELY THE
C...  BOUNDARY CONDITION IS SATISFIED
      WRITE(NO,1)NORUN,S,A,Y1,Y2
1     FORMAT(10X,' NORUN = ',I5,/,
     +       10X,'     S = ',F10.2,/,
     +       10X,'     A = ',F10.6,/,
     +       10X,'    Y1 = ',F10.6,/,
     +       10X,'    Y2 = ',F10.6,/)
      RETURN
      END

