N = 1000
for i = 0 to N - 1
n$ = str$(fibonacci(i))
j = val(left$(n$,1))
actual(j) = actual(j) +1
next
print
html "
| Digit | Actual | Expected |
"
for i = 1 to 9
html "
| ";i;" | ";using("##.###",actual(i)/10);" | ";using("##.###", frequency(i)*100);" |
"
next
html "
"
end
function frequency(n)
frequency = log10(n+1) - log10(n)
end function
function log10(n)
log10 = log(n) / log(10)
end function
function fibonacci(n)
b = 1
for i = 1 to n
temp = fibonacci + b
fibonacci = b
b = temp
next i
end function