Plot Utilities

7a1cb69c-737c-49e3-929f-b635bad5072d 52 4 5 true true 0 decimal

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

Plot Utilities

vx vy char size clr plotG n vx length : plot vx 1 el vy 1 el char size clr augment : i 2 n range plot plot vx i el vy i el char size clr augment stack : for plot 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 :

ListCollect

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 :

Levenberg Marquardt's second algorithm

x y β Tol βs LMSA Eps Tol : eta1 Eps sqrt : eta2 eta1 : h eta1 : stop 0 : 15 mat m Clear 1 : m x rows : n β rows : 12 mat β fx i 1 m range ff i el x i el ft y i el eval - : for ff 21 line : β derfh i 1 m range j 1 n range β j el β j el h + : jjh1 β fx i el : β j el β j el 2 h * - : jjh2 β fx i el : jjh0 i j el jjh1 jjh2 - 2 h * / : β j el β j el h + : 61 line for 11 line for jjh0 21 line : A n n Zeros : g n Zeros : J m n Zeros : 13 mat Fn 1 Zeros : dL 1 Zeros : dF 1 Zeros : F 1 Zeros : 14 mat f β fx : J β derfh : 12 mat A J transpose J * : g J transpose f * : ng g normi : 13 mat F f transpose f * 2 / : mu eta1 A Diag max * : nu 2 : 13 mat stop ¬ ng eta2 ≤ stop 1 : p A mu n identity * + invert g - (* : break try np β normi : nx eta2 β normi + : np eta2 nx * ≤ stop 2 : stop ¬ xnew β p + : fn xnew fx : Jn xnew derfh : Fn fn transpose fn * 2 / : dL p transpose mu p * g - (* 2 / : dF F Fn - : dL 1 el 0 > (dF 1 el 0 > (& β xnew : F Fn : J Jn : f fn : A J transpose J * : g J transpose f * : ng g normi : mu mu 13 / 12 dF 1 el dL 1 el / * 1 - (3^- Max * : nu 2 : 91 line mu mu nu * : nu 2 nu * : 21 line if 71 line 11 line if if 41 line if 11 line while βs β : 121 line :

x β Clear 1.000

x ft β 1 el β 3 el x β 4 el - / (β 2 el 1 +^* β 3 el x β 4 el - / (β 2 el^- exp * :

x 1234567891011111 mat :

y 0.4897 1.3079 2.84 3.0364 1.7709 0.5491 0.4606 0.3345 0.2073 0.0012 0.0024 111 mat :

XY x y augment : Data XY . 10 black plot : 21 line

ConjugateGradient => β 10573 - 14 mat transpose :

lower 'tol' Tol 0.00001 :

t0 1 time :

x y β Tol βs LMSA 8.512 44.14 55.43 52.06 - 41 mat

8.5418 40.7318 52.0514 48.7382 - 41 mat

1 time t0 - s 4.375

x β f β 1 el β 3 el x β 4 el - / (β 2 el 1 +^* β 3 el x β 4 el - / (β 2 el^- exp * :

x c 10 x ≤ (x 12 ≤ (& cases :

β 8.512 44.1447 55.4279 52.0571 - 41 mat

α 8.5418 40.7318 52.0514 48.7382 - 41 mat :

x α F α 1 el α 3 el x α 4 el - / (α 2 el 1 +^* α 3 el x α 4 el - / (α 2 el^- exp * :

x α F x c * x β f x c * Data 3 1 sys