10 DIM T(1002),S(1000),V(1000),A(1000)
20 INPUT "P = ";P
30 INPUT "ECCENTRICITY =";R
40 FOR N = 0 TO 1000
50 T(N) = N/10
60 NEXT N
70 FOR N = 1 TO 1000
80 S(N)=P/(1+R*((COS(T(N))-R)/(1-R*COS(T(N)))))
90 V(N) = (S(N)-S(N-1))/(T(N)-T(N-1))
100 A(N) = (V(N)-V(N-1))/(T(N)-T(N-1))
110 M = M + 1
120 NEXT N
130 FOR N = 1 TO M
140 PRINT "N =";N;TAB(20)"S =";S(N)
150 PRINT TAB(40);"V =";V(N)
160 PRINT TAB(60);"A =";A(N)
170 FOR PAUSE = 1 TO 100
180 NEXT PAUSE
190 NEXT N