<p> plot </p>

ab837d97-dbcb-4cf5-a9a2-78429657e56c 4 15 12 decimal

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

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 31 line for data r3 r4 r5 augment 21 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 :

Example of Non-linear Fit [Adapted from Smath source]. We will fit a model fnct using the CholeskyMatrixSolver

t0 1 time :

T 236.3 253.4 261.5 270.4 280.6 288.4 299.1 315.2 333.6 353.1 101 mat :

P 1510204060100200400760101 mat :

Install 'Y' non-linear i 1 P length range Y i el P i el log10 : for

Y 0 0.699 1 1.301 1.6021 1.7782 2 2.301 2.6021 2.8808 101 mat

i 1 P length range : NL i Y augment : 21 sys

NL o 5 blue plot

In this non-linear fit, the all thing is about installing 'Y' the non-linear component.

Just to show what equation we are fitting i.e: the model function to fit data points

P log10 b 1 el b 2 el T / + b 3 el T log10 * + b 4 el T 2^* + ≡

[1]

x φ 1000 x x log10 x 2^41 mat :

Memo 100 x x log10 x 2^51 mat

Prepare Cholesky n T length : k x φ length : 21 line

n k 21 mat 10421 mat

Cholesky  solver Num. j 1 k range i 1 n range X i j el T i el φ j el : for for Φ X transpose X * : d X transpose Y * : b Φ 1 -^eval d * : 41 line

Let recast the model function in [1]
to eliminate the log10 x fit 10 b x φ *^eval :

i 1 n range Pc i el T i el fit eval : for

Deviations @ support points Δ i 1 T length range res i el P i el T i el fit - eval : for P res augment 21 line :

b 0.288666886318 - 0.341588684153 - 149.943525856756 0.000241564906 41 mat

2.8867 - 0.3416 - 149.9436 0.0002 41 mat

T 236.3 253.4 261.5 270.4 280.6 288.4 299.1 315.2 333.6 353.1 101 mat

P 1510204060100200400760101 mat

Pc 0.98 5.26 10.25 19.66 37.86 59.07 101.67 204.6 401.83 753.59 101 mat

Δ 1 0.02 5 0.26 - 10 0.25 - 20 0.34 40 2.14 60 0.93 100 1.67 - 200 4.6 - 400 1.83 - 760 6.41 102 mat

Δ 2 BarSize

L H N xd U 0 : dx H L - N / : i 1 N 1 + range U i el L dx - (dx i * + : for U 41 line :

U 225375300 xd :

FIT i 1 U rows range vy i el U i el fit eval : for U vy augment 21 line :

1 time t0 - 1 s *

data T P augment o 5 black plot :

x c 1225 x ≤ (x 375 ≤ (& cases :

Complete plot: data points + the fitted model

data FIT 2 1 sys

Export this formula for the remaining of project.
Transport high decimals [b]... high 'x' in x^2 x q 10 b 2 el b 1 el x φ 1 el * b 2 el / x + (* b 3 el x log10 * + b 4 el x 2^* + (^:

data x q x c * 2 1 sys

300 fit 300 q 21 mat 106.083 106.083 21 mat