<p> ListCollect ... Jacobian </p>

c58b03b8-4a17-4929-9ce9-1de653331837 16 15 12 decimal

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

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 :

sol(init,Tol) ... AlgLib Solver for a system of nonlinear equations
adapted from Viacheslav original init Tol sol k 1 α f rows range u k el α k el : for α Jac u f u Jacobian : StepMax 0 : ε 1015 -^: 12 mat u init StepMax Tol α f α Jac al_nleqsolve : 41 line :

Bar plot unicolor style, bar width user. data dt BarSize v data 2 col : t v bar t dt - 0 t dt - v t dt + v t dt + 042 mat : j 1 data rows range j 1 ≡ B data j 1 el data j 2 el bar : B B data j 1 el data j 2 el bar stack : if 11 line for B 41 line :

al_nleqsolve suffers ln0, can't rescue ε 1012 -^:

XY 3 ε 4 0.1653 - 5 0.1313 6 0.219 7 0.8097 8 1.2045 9 1.302 10 2.1009 11 2.6004 12 3.7002 13 4.3001 112 mat :

1. selct the support points
2. collect for al_nleqsolve
3. collect variates [X, Y] list 24681115 mat : xy XY list ListCollect : X XY 1 col : Y XY 2 col : 12 mat 31 sys

data XY o 5 black plot : support xy . 14 blue plot : 21 line

data support 2 1 sys

The support points x xy 1 col : y xy 2 col : 12 mat

xy 4 0.1653 - 6 0.219 8 1.2045 10 2.1009 13 4.3001 52 mat

Automated al_nleqsolve ... NUMERIC/SYMBOLIC x α Φ α 1 el x * 1 x α 5 el - α 3 el / exp + ln (α 2 el^(* α 4 el + : init 0.01 11 0.1 - 115 mat transpose : α f ct 1 : j 1 x rows range Q ct el x j el α Φ y j el - : ct ct 1 + : 21 line for Q 31 line : 31 line

Model borrowed from Tesla project

t0 1 time :

Tol 106 -^:

α init Tol sol : 0.00717 0.80738 0.06782 0.16339 - 4.87525 51 mat

1 time t0 - s 25

Absolute residuals Residuals i 1 XY rows range Δ i el Y i el X i el α Φ - eval : for X Δ augment 21 line :

FIT U 015 0.01 range : i 1 U rows range vy i el U i el α Φ eval : for U vy augment 31 line :

data FIT 0.1653 + support 3 1 sys

Residuals 0.1 BarSize

ConjugateGradient diverges, hyperdifficult initialized ... NOT given, L-M failed

NO more exploration needed . . .

]]>

data FIT 0.1653 + support 31 sys