250817c8-45fe-4148-9aa1-45419088b1ba 6 4 6 decimal

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

Two definitions of the same piece-wise function

A 3.9083 103 -^* :

B 5.775 107 -^* - :

C 4.18301 1012 -^* - :

R.0 100 :

t R1.Pt100 R.0 1 A t * + B t 2^* + t 100 - (C * t 3^* + (* t 0 < R.0 1 A t * + B t 2^* + (* cases :

t R2.Pt100 t 0 < R.0 1 A t * + B t 2^* + t 100 - (C * t 3^* + (* R.0 1 A t * + B t 2^* + (* if :

Both produce same results

0 R1.Pt100 100

195.75 - R1.Pt100 20.3542

100 R1.Pt100 138.5055

0 R2.Pt100 100

195.75 - R2.Pt100 20.3542

100 R2.Pt100 138.5055

Newton-Raphson method (yes, I know there are plugins for this)

1 f g v e NR x1 x2 m x2 f x1 f - x2 x1 - / : v x1 x2 xNext x1 v x1 f - x1 x2 m / (+ : nMax 100 : n 0 : x1 g : x2 g e + : n nMax ≤ (x1 f v - abs e > (& ((n n 1 + : x1 v x1 x2 xNext : x2 x1 e + : 31 line while x1 81 line :

Solving for "t"

r T1.Pt100 t R1.Pt100 0 r 109 -^NR 11 line :

r T2.Pt100 t R2.Pt100 0 r 109 -^NR 11 line :

100 T1.Pt100 0

20 T1.Pt100 196.572 -

138 T1.Pt100 98.6675

100 T2.Pt100 #

20 T2.Pt100 #

138 T2.Pt100 #

Function with "if" doesn't work

Another implementation

t R.n R.0 1 A t * + B t 2^* + t 100 - (C * t 3^* + (* :

t R.p R.0 1 A t * + B t 2^* + (* :

t R3.Pt100 t 0 < t R.n t R.p if :

r T3.Pt100 r 100 < t R.n 0 r 109 -^NR t R.p 0 r 109 -^NR if 11 line :

r T4.Pt100 r 100 < t R.n 0 r 109 -^NR t R.p 0 r 109 -^NR if :

0 R3.Pt100 100

195.75 - R3.Pt100 20.3542

100 R3.Pt100 138.5055

100 T3.Pt100 0

20 T3.Pt100 196.572 -

138 T3.Pt100 98.6675

195.75 - R3.Pt100 T3.Pt100 195.75 -

Function without line() doesn't work

100 T4.Pt100 0

20 T4.Pt100 0

138 T4.Pt100 0