pian jindianjiao.f90

计算卫星位置的fortran程序

SUBROUTINE P_JDJ(M,E,EK)
IMPLICIT NONE
INTEGER M1
REAL*8 M,E,E0,EK
   M1=0
    E0=0
	EK=M+E0*SIN(E0)
!    WRITE(*,*)E0,E,E-E0
    DO WHILE(ABS((EK-E0))>=10E-10)
	  M1=M1+1
	  E0=EK
      EK=M+E*SIN(E0)
	END DO
!	WRITE(*,*)"M1=",M1,"E=",E
	RETURN
END