<content> Plot Utilities </content>

fe996056-08f2-4075-be52-f5fee1046fb6 2 Carlos Barradas G 4 5 false false 0 0 decimal

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

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 :

<content> ListCollect </content>

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 :

<content> Levenberg Marquardt's second algorithm </content>

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 Zero : 13 mat Fn 1 Zero : dL 1 Zero : dF 1 Zero : F 1 Zero : 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 ft β 1 el x β 2 el - (2^2 β 3 el 2^* / - exp * :

M x,time,min y,dye concentration, g/gal 30507011014116319132122242925262618278293330360400440182 mat :

x M 51511 submatrix :

y M 51522 submatrix :

x 1114161921242526272933111 mat

y 0131322292618830111 mat

Tol 101 -^:

t0 1 time :

x y β Tol βs LMSA #

1 time t0 - s 12.063