09e01732-0d21-4023-a247-0377e0418d3d 9 15 12 decimal

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

x Γ x Gamma :

α 1.3536 :

β 153.92 :

x g 1 β α^α Γ * / x α 1 -^e x - β /^* (* :

MCD sanity

x Y D x g :

init 0 :

sol init 020002000 x Y D dn_GearsBDF :

t0 1 time :

Smath 6179 cumulative ODE solver ... 2000 data

sol 2 1 sys

PDF for cumulative CDF

x g

0 g 0

sol sol rows row 2000.000 1.000 12 mat

sol rows 2001

1 time t0 - 1.1 s *

The hyperfast Golden-Ratio search [min/max]

x f x g : a b G 32 / b * 12 / a * - 12 / 5 sqrt * a b - (* + : 21 line

TOL 106 -^:

Code adapted from Valery Ochov Mathcad 11 f a b test TOL Extrema y1 y2 ≥ max in range ≡ y1 y2 ≤ min in range ≡ x1 a b G eval : y1 x1 f eval : x2 b a G eval : y2 x2 f eval : 41 line iter 0 : iter 1000 < a b - abs TOL ≤ break test b x2 : x2 x1 : y2 y1 : x1 a b G eval : y1 x1 f eval : 51 line a x1 : x1 x2 : y1 y2 : x2 b a G eval : y2 x2 f eval : 51 line if iter iter 1 + : 21 line if while a b + 2 / iter 21 mat 61 line :

t0 1 time :

X f 0128 y1 y2 ≥ TOL Extrema 1 el :

1 time t0 - 0.2 s *

x g X X f + 100 red 1 5 mat 2 1 sys

1. All modules/functions in "Symbolic" optimiz. 2. Once in target: shrink the range 3. respect the search y1≥y2≡"max in range" y1≤y2≡"min in range"

X X g 21 mat 54.42607 0.00355 21 mat