Thiele's interpolation formula: Difference between revisions
m
→{{header|Julia}}
Line 718:
idx = (N - 1 - n) * div(N - n, 2) + i
if !haskey(rhocache, idx)
rhocache[idx] = (x[i + 1] - x[i + n + 1]) / (rho(x, y, rhocache, i, n - 1) -
rho(x, y, rhocache, i + 1, n - 1)) + rho(x, y, rhocache, i + 1, n - 2)
end
Line 728:
return 1.0
end
rho(x, y, r, 0, n) - rho(x, y, r, 0, n - 2) + (xin - x[n + 1]) / thiele(x, y, r, xin, n + 1)
end
function thiele_tables()
for i in 1:N
xval_table[i] = (i - 1) * step
tsin_table[i] = sin(xval_table[i])
tcos_table[i] = cos(xval_table[i])
Line 749:
3.1415976324027914
</pre>
=={{header|Kotlin}}==
| |||