Thiele's interpolation formula: Difference between revisions

 

=={{header|Julia}}==

{{trans|C}}

const N2 = N * div(N - 1, 2)

const xval_table = zeros(Float64, N)

const tsin_table = zeros(Float64, N)

idx = (N - 1 - n) * div(N - n, 2) + i

if !haskey(rhocache, idx)

rho(x, y, rhocache, i + 1, n - 1)) + rho(x, y, rhocache, i + 1, n - 2)

end

function thiele_tables()

for i in 1:N

tsin_table[i] = sin(xval_table[i])

tcos_table[i] = cos(xval_table[i])

thiele_tables()

</lang>{{output}}<pre>

</pre>