n = 15108097
k = n Mod 10000
T = n Mod 100
pi = 4 * Atn(1)
g = 10
x = 2 - 1 / k
y = 0.5 + 1 / n
v = 6
a = g * x ^ 2 / 2 / v ^ 2
b = -x
c = g * x ^ 2 / 2 / v ^ 2 + y
c = a + y
N1 = (-b + Sqr(b ^ 2 - 4 * a * c)) / (2 * a)
'N1 = Atn(N1)
aN1 = Atn(N1)
DaN1 = 180 * aN1 / pi
MsgBox DaN1

N2 = (-b - Sqr(b ^ 2 - 4 * a * c)) / (2 * a)
'N2 = Atn(N2)
aN2 = Atn(N2)
DaN2 = 180 * aN2 / pi
MsgBox DaN2

'
y = -1
a = g * (1 + N1 ^ 2) / 2 / v ^ 2
b = -N1
c = y
X1 = (-b + Sqr(b ^ 2 - 4 * a * c)) / (2 * a)
MsgBox X1

a = g * (1 + N2 ^ 2) / 2 / v ^ 2
b = -N2
c = y
X2 = (-b + Sqr(b ^ 2 - 4 * a * c)) / (2 * a)
MsgBox X2

'N1 = Atn(N1)

t1 = X1 / v / Cos(aN1)
MsgBox t1

'N2 = Atn(N2)
t2 = X2 / v / Cos(aN2)
MsgBox t2


Dim f(99)
a = 0
b = X1
'b = X2
n = 88
h = (b - a) / n
'
For kk = 0 To n
x = a + kk * (b - a) / n
f(kk) = Sqr(1 + (N1 - g * x * (1 + N1 ^ 2) / v ^ 2) ^ 2)
'f(kk) = Sqr(1 + (N2 - g * x * (1 + N2 ^ 2) / v ^ 2) ^ 2)
Next kk
'
i2j = 0
For j = 1 To n / 2 - 1
i2j = i2j + f(2 * j)
Next j
'
i2jm1 = 0
For j = 1 To n / 2
i2jm1 = i2jm1 + f(2 * j - 1)
Next j
'
L1 = (f(0) + 2 * i2j + 4 * i2jm1 + f(n)) * h / 3
'
MsgBox L1

'b = X1
b = X2
n = 88
h = (b - a) / n
'
For kk = 0 To n
x = a + kk * (b - a) / n
'f(kk) = Sqr(1 + (N1 - g * x * (1 + N1 ^ 2) / v ^ 2) ^ 2)
f(kk) = Sqr(1 + (N2 - g * x * (1 + N2 ^ 2) / v ^ 2) ^ 2)
Next kk
'
i2j = 0
For j = 1 To n / 2 - 1
i2j = i2j + f(2 * j)
Next j
'
i2jm1 = 0
For j = 1 To n / 2
i2jm1 = i2jm1 + f(2 * j - 1)
Next j
'
L2 = (f(0) + 2 * i2j + 4 * i2jm1 + f(n)) * h / 3
'
MsgBox L2

aV1 = L1 / t1
MsgBox aV1

aV2 = L2 / t2
MsgBox aV2