OPEN "tvp.txt" FOR INPUT AS #1
OPEN "press.txt" FOR OUTPUT AS #2
OPEN "veloc.txt" FOR OUTPUT AS #3
OPEN "displ.txt" FOR OUTPUT AS #4
CLS
DIM tin(20), pin(20), t(150), p(150), v(150), d(150), tim(50), pre(50), vel(50)

pi = 3.14159265#

REM for 6PPC with a 68 gr bullet and a 22" barrel.

dia = .243
g = 386
mass = 68 / (7000 * g)
area = pi * dia ^ 2 / 4
fudge = .78
INPUT #1, n
FOR i = 1 TO n
INPUT #1, tin(i), pin(i)
tin(i) = (tin(i) - 1.29) * 1.5 / 5880
pin(i) = (pin(i) - 5.88) * 50000 / 2.03
NEXT i
FOR i = 1 TO n
PRINT USING "   ##  #.####### ########.#"; i; tin(i); pin(i)
NEXT i

REM Make a table of 146 points.
PRINT #2, USING "####"; 145
PRINT #3, USING "####"; 145
PRINT #4, USING "####"; 145


tmax = .00145
dt = tmax / 145
t = -dt
s = 2
v = 0
n = 1
FOR i = 1 TO 145
t = t + dt
IF t > tin(n + 1) THEN n = n + 1
m = (pin(n + 1) - pin(n)) / (tin(n + 1) - tin(n))
b = pin(n) - m * tin(n)
p = m * t + b
p = fudge * p
v = v + p * area * dt / mass
s = s + v * dt
PRINT USING " #### #.#######, ######. ######.### ####.###"; n; t; p; v / 12; s
PRINT #2, USING " #.#######, #######.#"; t; p / fudge
PRINT #3, USING " #.#######, #####.###"; t; v / 12
PRINT #4, USING " #.#######, #####.###"; t; s
t(i) = t: p(i) = p: v(i) = v / 12: d(i) = s



NEXT i

nt = 0
dis = 0
FOR dis = 1 TO 30
FOR i = 1 TO 145
IF d(i) > dis THEN 10 ELSE 20
10 m = (t(i) - t(i - 1)) / (d(i) - d(i - 1))
b = t(i) - m * d(i)
tim(dis) = m * dis + b

m = (p(i) - p(i - 1)) / (d(i) - d(i - 1))
b = p(i) - m * d(i)
pre(dis) = m * dis + b

m = (v(i) - v(i - 1)) / (d(i) - d(i - 1))
b = v(i) - m * d(i)
vel(dis) = m * dis + b


GOTO 30
20 NEXT i
30 PRINT USING "#.######   ###.#   #####   ######"; tim(dis); dis; vel(dis); pre(dis) / fudge
NEXT dis