<p> Thiele modules </p>

1a90eaac-acdb-4557-8a69-f08298c81fbe 4 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

XY 2.6534 0.0001 2.7823 0.0153 2.9112 0.2767 3.04 1.556 3.1689 4.1571 3.2978 6.6836 3.4267 7.4801 3.5556 6.4108 3.6844 4.4876 3.8133 2.6778 3.9422 1.3961 4.0711 0.6393 4.2 0.2509 4.3288 0.077 4.4577 0.0134 152 mat :

X 0.149 0.299 0.448 0.597 0.747 0.896 1.045 1.195 1.344 1.4933 1.6427 1.792 1.9413 2.0907 2.24 151 mat :

Y 1.5 1.5 1.5 1.5 1.5 1.25 1.071 0.938 0.833 0.75 0.6818 0.625 0.5769 0.5357 0.5 151 mat :

Apply . . . Thiele interpolation

In lieu of fitting XY Gaussian i.e: a resident function

Thiele interpolation is offered for concluding the project.

]]> <p> Thiele modules </p>

X Y Δ Thiele ε 1018 -^: n X rows : i 1 n range M i 1 el Y i el : for j 2 n range i j n range dM M i j 1 - el M j 1 - j 1 - el - : dM dM abs ε < ε M i j el X i el X j 1 - el - dM / : if : 21 line for for i 1 n range c i 1 el X i el : c i 2 el M i i el : 21 line for Δ 0 ≡ c M if 61 line :

X XY 1 col : Y XY 2 col : 12 mat

length Thiele j XY rows : 15

Thiele on length of support points 'j' K X x Cfr N j 1 + el K j 1 + el : c j el N j el : a j el K j el : jj j 1 - (0 range c jj el a jj 1 + el : a jj el K jj el x X jj el - c jj el / + : 21 line for a 1 el 51 line :

recast the Cf on the system  [K,X,x] x Cf K X x Cfr :

Thiele K's coefficients @ support points X X Y 0 Thiele 2.6534 0.0001 2.7823 8.480263 2.9112 0.017077 - 3.04 9.05218 - 3.1689 0.577568 3.2978 0.77966 3.4267 1.85357 3.5556 0.193387 - 3.6844 4.998246 - 3.8133 0.439989 3.9422 6.803144 - 4.0711 0.028857 - 4.2 8.079122 4.3288 0.827207 4.4577 1.543274 152 mat

Thiele K's coefficients. K X Y 0 Thiele 2 col :

K K 6 round vectorize eval :

x λ 1 X 1 el x ≤ (x X 15 el ≤ (& cases :

Very slow plot x D x Cf x diff x λ * :

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

support XY . 12 black plot :

support x Cf x λ * 2 1 sys