home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 10 / 10.iso / l / l460 / 2.ddi / SRC.DI$ / YPRIME.F < prev    next >
Encoding:
Text File  |  1993-03-26  |  489 b   |  26 lines
  1. C The actual YPIME subroutine in FORTRAN
  2.  
  3.  
  4.       SUBROUTINE YPRIME(YP, T, Y)
  5.       REAL*8 YP(4), T, Y(4)
  6.  
  7.       REAL*8 MU, MUS, R1, R2
  8.       
  9.       MU = 1.0/82.45
  10.       MUS = 1.0 - MU
  11.  
  12.       R1 = SQRT((Y(1)+MU)**2 + Y(3)**2)
  13.       R2 = SQRT((Y(1)-MUS)**2 + Y(3)**2)
  14.  
  15.       YP(1) = Y(2)
  16.       YP(2) = 2*Y(4) + Y(1) - MUS*(Y(1)+MU)/(R1**3) - 
  17.      & MU*(Y(1)-MUS)/(R2**3)
  18.  
  19.       YP(3) = Y(4)
  20.       YP(4) = -2*Y(2) + Y(3) - MUS*Y(3)/(R1**3) - 
  21.      & MU*Y(3)/(R2**3)
  22.  
  23.       RETURN
  24.       END
  25.  
  26.