data

22316b2a-513e-4706-94aa-3a6fa1905e6f 4 6 12 decimal

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

data

data 0 2113.3 1 2111.7 2 2109.2 3 2105.2 4 2099.4 5 2092.1 6 2084.7 7 2078.6 8 2075.1 9 2073.9 10 2074.1 11 2074.4 12207413 2072.5 14 2070.1 15 2067.2 16 2064.4 17206218 2060.5 19 2059.6 202 mat :

Thiele(X,Y,Δ) ... algo

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 :

ListCollect(data,list) ... plot(data,char,size,clr)

data list ListCollect collects index list select list : 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 :

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 :

list 13579111315171920111 mat transpose : xy data list ListCollect : raw data o 5 black plot : support xy . 15 blue plot : 41 line

data rows 20

Address the system

support raw 2 1 sys

X xy 1 col : Y xy 2 col : XY X Y augment : 31 sys

'j' governs"
Cfr(K,X,x) j XY rows :

j 11

Cfr(K,X,x) ..... Export Optimiz NONE

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 :

Expand symbolic NONE.
Export the project
c/w X,K fill-in K X x Cfr K 1 el x 0 - K 2 el x 2 - K 3 el x 4 - K 4 el x 6 - K 5 el x 8 - K 6 el x 10 - K 7 el x 12 - K 8 el x 14 - K 9 el x 16 - K 10 el x 18 - K 11 el / + / + / + / + / + / + / + / + / + / +

recast the Cfr(K,X,x) algo style dependent upon 'K' x f K X x Cfr :

@ this stage, f(x) is algo style, awaiting for next 'K' calculated below

Calculate/assign 'K' from ... Thiele(X,Y,0)

Thiele K's coefficients
@ support points X X Y 0 Thiele 0 2113.3 2 0.4878 - 4 9.9982 6 0.4556 8 18.271 - 10 2.4424 - 12 1.5877 14 2.8399 16 7.5897 - 18 2.2942 - 19 12.293 - 112 mat

K X Y 0 Thiele 2 col : 2113.3 0.488 - 9.998 0.456 18.271 - 2.442 - 1.588 2.84 7.59 - 2.294 - 12.293 - 111 mat

x c 11 - x ≤ (x 20 ≤ (& cases :

'data' are the ksmoothed corn thermies ... Thiele fitted

raw x f x c * support 3 1 sys Export Corn Thermies fitted Thiele Continued fraction

X data 1 col :

x Corn K 1 el x X 1 el - K 2 el x X 3 el - K 3 el x X 5 el - K 4 el x X 7 el - K 5 el x X 9 el - K 6 el x X 11 el - K 7 el x X 13 el - K 8 el x X 15 el - K 9 el x X 17 el - K 10 el x X 19 el - K 11 el / + / + / + / + / + / + / + / + / + / + :

'data' are the ksmoothed corn thermies ... Thiele fitted

raw x Corn x c * support 3 1 sys

x Corn x 019 int 39536

rational model  ... α from Mathcad GenfitMatrix.
No way found to initialize al_nleqsolve x α Φ α 1 el α 2 el x * + α 3 el x 2^* + α 4 el x 3^* + α 5 el x 4^* + α 6 el x 5^* + 1 α 7 el x * + α 8 el x 2^* + α 9 el x 3^* + α 10 el x 4^* + α 11 el x 5^* + / :

α 2113.3 7698867.75 - 2325181.6 309845.84 - 18772.336 461.17 - 3643.464 - 1098.816 146.84 - 8.933 0.22066 - 111 mat :