9461744b-9cfb-4cee-a7b8-7755cd01e3b9 48 6 9 decimal

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

vx vy a b c plotG a char ≡ b size ≡ c color ≡ 13 mat n vx length : plot vx 1 el vy 1 el a b c augment : i 2 n range plot plot vx i el vy i el a b c augment stack : for plot 51 line :

Tr φ 36 Tr / 42 Tr ln * + 35 - Tr 6^- :

Tr Ψ 0.0364 Tr φ * Tr log10 - :

Tc Pc α T Pvap 10 Pc log10 0.1183 T Tc / φ * - 7 T Tc / log10 * + α 7 - (T Tc / Ψ * -^:

Tc 433.8 556.4 553562 619.6 51 mat :

Pc 31.6 4540 48.6 3851 mat 101.33 * :

α 6.771 6.742 6.844 6.848 7.330 51 mat :

z 0.35 0.1 0.1 0.1 0.35 51 mat :

P 101.33 :

T 350 :

T P Tc Pc α k Tc Pc α T Pvap P / :

i 1 Tc length range k i el T P Tc i el Pc i el α i el k : for

k 6.798095 1.004688 0.889606 0.903627 0.140563 51 mat

x unknown 1.8 - :

Φ f z i el k i el 1 - (* k i el 1 - (Φ * 1 + / (i 1 z length sum :

δ 1.1636 :

0 x ≤ x 1.5 ≤ & x f rejected if δ 200 δ 200 - 2 2 mat x unknown 3 1 sys

δ f 7445.470396

Plot spec's 0 x ≤ (x 1.5 ≤ (& x f rejected if δ 200 δ 200 - 22 mat x unknown 31 sys

x φ x Ψ 2 1 sys

0 f 1 f 21 mat 1.708323 1.864072 - 21 mat

Given f(x) => solve for 'x'

f L H N Inverse dx H L - N / : i 1 N 1 + range ω i el L dx - (dx i * + : for j 1 ω rows range sol j el x f ω j el - 0 ≡ x 0 1.1636 solve eval : for ω sol augment 51 line :

Sol f 1.7 1.8 - 35 Inverse :

Working range of f(Φ)

0 x ≤ ( x 1 ≤ ( & x f rejected if 1.1636 3 1.1636 3 - 2 2 mat 2 1 sys

Sanity of "Sol" wrt support points Φ f 1.7 ≡ Φ 0 0.1 solve Φ f 1 ≡ Φ 0 0.1 solve Φ f 0 ≡ Φ 0.4 0.6 solve Φ f 1 - ≡ Φ 0.8 1 solve Φ f 1.8 - ≡ Φ 0.95 1 solve 51 mat 0.000695 0.087618 0.481969 0.895438 0.995056 51 mat

Sol 1.7 0.000695 1.6 0.009489 1.5 0.019195 1.4 0.029956 1.3 0.041944 1.2 0.05537 1.1 0.07049 1 0.087618 0.9 0.107138 0.8 0.129522 0.7 0.155338 0.6 0.185263 0.5 0.220059 0.4 0.260505 0.3 0.307234 0.2 0.360414 0.1 0.419315 0 0.481969 0.1 - 0.545337 0.2 - 0.606151 0.3 - 0.661926 0.4 - 0.711424 0.5 - 0.754488 0.6 - 0.791601 0.7 - 0.823506 0.8 - 0.850987 0.9 - 0.874759 1 - 0.895438 1.1 - 0.913536 1.2 - 0.929471 1.3 - 0.943586 1.4 - 0.95616 1.5 - 0.967422 1.6 - 0.97756 1.7 - 0.986728 1.8 - 0.995056 362 mat

vx Sol 1 col : vy Sol 2 col : sol vx vy . 9 black plotG : 31 line

Global solution => think modelling

sol 1.8 - 2 1.8 - 0.5 - 2 2 mat 1.7 0.5 1.7 0.5 - 2 2 mat 3 1 sys

User on-line "solve/collect" Φo OnLine 1.8 - Φo ≤ 1.7 ≤ Φ Φ f Φo ≡ Φ 0 1.1636 solve eval : i 1 z length range x i el z i el k i el 1 - (Φ * 1 + / : y i el k i el x i el * : 21 line for header Φo Φsol x.....y 13 mat : header Φo Φ x y augment 13 mat stack 51 line :

Φo 1.25 - :

Φo OnLine Φo Φsol x.....y 1.25 - 0.9367 0.0544 0.37 0.0996 0.1 0.1115 0.0992 0.1099 0.0993 1.7955 0.2524 52 mat 23 mat