plotG(vx,vy,char,size,color)

3faf3652-8fed-46f1-81ec-e044979c15c9 26 8 12 decimal

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

..... Pt100 Thiele XFR .....

continuous inverse interpolation

source algorithm

]]> plotG(vx,vy,char,size,color)

vx vy char size color plotG n vx length : plot vx 1 el vy 1 el char size color augment : i 2 n range plot plot vx i el vy i el char size color augment stack : for plot 41 line :

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 :

Thiele(X,Y,Δ)

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 :

eit90 273.15 - 0.00006230541 268.935 - 0.000155487 259.3467 - 0.00119007 248.5939 - 0.00844974 218.7916 - 0.09171804 189.3442 - 0.21585975 38.8344 - 0.84414211 0.01 1 29.7646 1.11813889 156.5985 1.60980185 231.928 1.89279768 419.527 2.5689173 660.323 3.3760086 961.78 4.28642053 1064.43 4.57255961 1084.5 4.62715523 162 mat :

data List collects index list select 34567891011121314112 mat : ct 1 : i select a ct el data i 1 el : b ct el data i 2 el : ct ct 1 + : 31 line for a b augment 51 line :

select the support points XY eit90 List :

length Thiele j XY rows : 12

Cfr(K,X,x)

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 :

The continued fraction ... Optimiz => NONE K X x Cfr K 1 el x X 1 el - K 2 el x X 2 el - K 3 el x X 3 el - K 4 el x X 4 el - K 5 el x X 5 el - K 6 el x X 6 el - K 7 el x X 7 el - K 8 el x X 8 el - K 9 el x X 9 el - K 10 el x X 10 el - K 11 el x X 11 el - K 12 el / + / + / + / + / + / + / + / + / + / + / +

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

data set vx eit90 1 col : vy eit90 2 col : 21 sys

support points X XY 1 col : Y XY 2 col : 21 sys

support X Y . 12 black plotG : data vx vy . 10 black plotG : x c 1 273.15 - x ≤ (x 1200 ≤ (& cases : 31 line

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

x XFR x c * data 259.3467 - 0.00119007 + 40 blue 1 5 mat 961.78 4.28642053 + 40 blue 1 5 mat 5 1 sys

Carlos algorithm: 2018-05-10
Thiele table of inverse differences i 1 Y rows range k i 1 el Y i el : for j 2 Y rows range i j Y rows range k i j el X j 1 - el X i el - k j 1 - j 1 - el k i j 1 - el - / : iter j : 21 line break try for for 21 line

passed to Pt100 CALCULATOR

i 1 iter range Thiele i el k i i el 12 round : for

X Y Δ Thiele

Thiele K 12 mat 0.00119007 1481.16925424985 0.028845074566 - 1311.68667916945 - 2.758880829033 87.121262673755 33.319386171553 27.672759469917 - 11.033135031489 - 329.084599654368 - 2.244251954285 179.627799615139 121 mat 0.00119007 1481.16925424985 0.028845074566 - 1311.68667916944 - 2.758880829033 87.121262673727 33.319386169641 27.672759482964 - 11.03313503424 - 329.084601051645 - 2.244251932054 179.627801033821 121 mat 12 mat

With the Try / on error function, division by zero is avoided when the inverse differences of a column are all the same, and the loop is exited.

k 11 mat Extend Thiele table ≡

kernel internal counter iter 12

k 0.0012 0000000 0.0084 1481.1693 000000 0.0917 447.9842 0.0288 - 00000 0.2159 326.094 0.0513 - 1311.6867 - 0000 0.8441 261.5953 0.172 - 1257.132 - 2.7589 0001 259.6657 0.2035 - 1252.6011 - 3.2047 87.1213 00 1.1181 258.8402 0.2277 - 1249.7616 - 3.5383 88.0143 33.3194 0 1.6098 258.574 0.3314 - 1240.6521 - 4.8701 92.5705 28.736 27.6728 - 1.8928 259.7128 0.3934 - 1236.3532 - 5.5921 95.567 27.4596 34.5003 - 2.5689 264.387 0.5491 - 1226.9619 - 7.1865 103.5242 25.5757 50.3328 - 3.376 272.5094 0.752 - 1215.6591 - 8.8482 114.8177 23.8411 66.5264 - 4.2864 284.9617 1.0118 - 1200.9913 - 10.399 130.9681 21.9348 81.8663 - 128 mat

273.15 - 268.935 - 259.3467 - 248.5939 - 218.7916 - 189.3442 - 38.8344 - 0.01 29.7646 156.5985 231.928 419.527 660.323 961.78 1064.43 1084.5 161 mat_NBP Helium TP Hydrogen TP Neon TP Oxygen TP Argon TP Mercury TP Water MP Gallium FP Indium FP Tin FP Zinc FP Aluminum FP Silver FP Gold FP Copper 161 mat 12 mat