Thiele's interpolation formula: Difference between revisions
m
→{{header|Sidef}}: updated code
(Completed restoration of formulae made invisible at 18:05, 31 July 2016 (one ρ character, and three π characters restored today)) |
m (→{{header|Sidef}}: updated code) |
||
Line 1,039:
{{trans|Python}}
<lang ruby>func thiele(x, y) {
var ρ =
for i in ^(ρ.end) {
ρ[i][1] = ((x[i] - x[i+1]) / (ρ[i][0] - ρ[i+1][0]))
}
for i
for j
ρ[j][i] = (((x[j]-x[j+i]) / (ρ[j][i-1]-ρ[j+1][i-1])) + ρ[j+1][i-2])
}
}
var ρ0 = ρ[0]
func t(xin) {
var a = 0
for i
a = ((xin - x[i-1]) / (ρ0[i] - ρ0[i-2] + a))
}
Line 1,061:
return t
}
# task 1: build 32 row trig table
var xVal =
var tSin = xVal.map { .sin }
var tCos = xVal.map { .cos }
var tTan = xVal.map { .tan }
# task 2: define inverses
var iSin = thiele(tSin, xVal)
var iCos = thiele(tCos, xVal)
var iTan = thiele(tTan, xVal)
# task 3: demonstrate identities
say 6*iSin(0.5)
Line 1,079:
{{out}}
<pre>
3.14159265358979323846438729976818601771260734312
3.14159265358979323846157620314930763214337987744
3.14159265358979323846264318595256260456200366896
</pre>
| |||