' Trapezoidal rule
' for function f(x)= y(x) = x^6
Dim f(999)
s = 16108089
T = s Mod 100
a = 0
b = 1
n = 88
n = 2 * T
h = (b - a) / n
'
For k = 0 To n
x = a + k * h
f(k) = x ^ 6
Next k
'
i = 0
For k = 1 To n - 1
i = i + h * f(k)
Next k
i = i + (f(0) + f(n)) * h / 2
MsgBox i
d1 = Abs(i - 1 / 7)
'MsgBox d1
d2 = Abs(30 / (12 * n ^ 2))
MsgBox d2
dd = Abs(d1 - d2)
'MsgBox dd