6f6479e7-5946-49a3-a287-e9da95fb44b5 10 4 5 false false 0 decimal

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

PRIMER 3.8

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 :

A 0 0.018 00000 0.125 - 0.045 - 0.045 - 000 0.797 - 0.797 - 0000045 mat :

B 014988274274007153530003943940000045 mat :

z 0.35 0.1 0.1 0.1 0.35 51 mat :

P 101.33 :

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

Tr Ψ 0.0364 Tr φ * Tr log10 - :

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

T x γRK Nc x length : C A B T / + : nul Nc el 0 : C C nul transpose stack : C C C transpose + : j 1 Nc range s C j k el x k el * 1 x j el - (* (k 1 Nc sum : m 1 Nc 1 - range n m 1 + Nc range y j n ≡ (j m ≡ (| 0 C m n el x m el * x n el * if : s s y - : 21 line for for γ j el s exp : 31 line for γ 71 line :

x P T kγ γ T x γRK : k γ Tc Pc α T Po * P / vectorize : k 31 line :

T 350 :

Variables - initial guess

Φ 0.5 :

Initial guess

k Tc Pc α T Po vectorize :

x z k 1 - (Φ * 1 + / vectorize :

y k x * vectorize :

k 688.851 101.8051 90.1438 91.5645 14.2433 51 mat

Initial guess

x 0.001 0.0019 0.0022 0.0022 0.0459 51 mat

y 0.699 0.1981 0.1978 0.1978 0.6541 51 mat

Equations - using initial guess

f z Φ y * 1 Φ - (x * + (- : 4.8925 1016 -^* - 1.3769 1017 -^* 1.7517 1017 -^* 1.6984 1017 -^* 051 mat

g y x / vectorize x P T kγ - : 681.7827 100.7947 89.2397 90.6594 14.1025 51 mat

h x sum (y sum (- : 1.8935 -

all f,g,h values should be close to zero in the solution for Φ,x,y

HOW TO SOLVE for Φ,x,y ???

Here is the solution from other software

x 0.049 0.0973 0.0953 0.1047 0.6538 51 mat :

y 0.6044 0.1023 0.104 0.096 0.0933 51 mat :

Φ 0.542 :

and the functions quite close to zero as required

f z Φ y * 1 Φ - (x * + (- : 2.68 105 -^* - 1105 -^* - 1.54 105 -^* - 1.54 105 -^* 9106 -^* - 51 mat

g y x / vectorize x P T kγ - : 0.0054 - 0.0003 0.0002 - 0.0007 - 2.9043 105 -^* - 51 mat

h x sum (y sum (- : 0.0001