Adaptive Plot

f18cb960-2268-4668-b695-e831c8372bcd 120 4 5 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 :

Adaptive Plot

z w v CAngle a w z - : b v w - : 12 mat a Re i a Im * - (b * Re a abs b abs * / acos Re eval 0 try 21 line :

1 f a b APlot ξ f a b APlotAngle APlotDist MaxIters APlotNG / APlotNG APlot :

1 f a b α p Μ Ν APlot δ b a - Μ / : k 1 Μ 1 + (range T k el a k 1 - (δ * + : for Z T f vectorize eval : d p Z max Z min - abs * eval : ng 1 Ν range N T length : k 1 N 2 - (range A k el Z k el Z k 1 + el Z k 2 + el CAngle : for k 1 N 1 - (range D k el Z k 1 + el Z k el - abs : for A max α ≤ (D max d ≤ (& eval break N T length : k 0 : t 0 : κ 0 : τ 0 : 15 mat n 0 N 1 - 2 / 1 - (range t k k 1 + : el T 2 n * 1 + el : θ Z 2 n * 1 + el Z 2 n * 2 + el Z 2 n * 3 + el CAngle : θ α > (Z 2 n * 2 + el Z 2 n * 1 + el - abs d > (| t k k 1 + : el 0.5 T 2 n * 1 + el T 2 n * 2 + el + (* : if t k k 1 + : el T 2 n * 2 + el : θ α > (Z 2 n * 2 + el Z 2 n * 3 + el - abs d > (| t k k 1 + : el 0.5 T 2 n * 2 + el T 2 n * 3 + el + (* : if 51 line for N 2 mod 0 ≡ t k k 1 + : el T N 1 - el : θ α > (Z N el Z N 1 - el - abs d > (| t k k 1 + : el 0.5 T N 1 - el T N el + (* : if 21 line if T t T N el stack eval : Z T f vectorize eval : 51 line if 41 line for Z 61 line :

CSspline

Index in X for the interval containing x, assuming that  
X is in ascending order. x\032E\ X\032E\ CS.L n\0324\ X\032E\ length eval : x\032E\ X\032E\ 1 el X\032E\ n\0324\ el lele x\032E\ 2 X\032E\ 1 el * X\032E\ 2 el - X\032E\ 2 el ltlt 1 x\032E\ X\032E\ n\0324\ 1 - el 2 X\032E\ n\0324\ el * X\032E\ n\0324\ 1 - el - lelt n\0324\ 1 - n\0324\ x\032E\ X\032E\ ≤ vectorize 11 findrows rows - if if CS.L: Value outside of range. error if 21 line :

t\032E\ h\032E\ y\032E\.0 y\032E\.1 y\032E\''.0 y\032E\''.1 CS.1 t\032E\ t\032E\ 1 - 2 / y\032E\''.0 t\032E\ 1 + 3 / y\032E\''.1 y\032E\''.0 - (* + (* h\032E\ * h\032E\ * y\032E\.1 + y\032E\.0 - (* y\032E\.0 + eval :

Cubic Spline interpolation function. M\032E\ x\032E\ CS X\032E\ Y\032E\ V\032E\ T\032E\ 14 mat M\032E\ eval : k\0324\ x\032E\ X\032E\ CS.L : h\0324\ X\032E\ k\0324\ 1 + el X\032E\ k\0324\ el - : x\032E\ X\032E\ k\0324\ el - h\0324\ / h\0324\ Y\032E\ k\0324\ el Y\032E\ k\0324\ 1 + el V\032E\ k\0324\ el V\032E\ k\0324\ 1 + el CS.1 41 line :

Cubic Spline Derivative. M\032E\ x\032E\ CS.D X\032E\ Y\032E\ V\032E\ T\032E\ 14 mat M\032E\ eval : k\0324\ x\032E\ X\032E\ CS.L : h\0324\ X\032E\ k\0324\ 1 + el X\032E\ k\0324\ el - : t\0324\ x\032E\ X\032E\ k\0324\ el - h\0324\ / : 2 t\0324\ * 1 - (V\032E\ k\0324\ el * t\0324\ 2^13 / - (V\032E\ k\0324\ 1 + el V\032E\ k\0324\ el - (* + (h\0324\ 2 / * Y\032E\ k\0324\ 1 + el Y\032E\ k\0324\ el - h\0324\ / + 51 line :

Cubic Spline Second Derivative. M\032E\ x\032E\ CS.D2 X\032E\ Y\032E\ V\032E\ T\032E\ 14 mat M\032E\ eval : k\0324\ x\032E\ X\032E\ CS.L : V\032E\ k\0324\ el V\032E\ k\0324\ 1 + el V\032E\ k\0324\ el - (x\032E\ X\032E\ k\0324\ el - X\032E\ k\0324\ 1 + el X\032E\ k\0324\ el - / * + 31 line :

Cubic Spline Integral. M\032E\ a\032E\ b\032E\ CS.I a\032E\ b\032E\ > M\032E\ b\032E\ a\032E\ CS.I - a\032E\ b\032E\ ≡ 0 X\032E\ Y\032E\ V\032E\ T\032E\ 14 mat M\032E\ eval : i\0324\ a\032E\ X\032E\ CS.L : j\0324\ b\032E\ X\032E\ CS.L : I\0324\ 0 : 13 mat k\0324\ i\0324\ j\0324\ range h\0324\ X\032E\ k\0324\ 1 + el X\032E\ k\0324\ el - : I\0324\ I\0324\ h\0324\ 2 / Y\032E\ k\0324\ el Y\032E\ k\0324\ 1 + el + h\0324\ 2^12 / V\032E\ k\0324\ el V\032E\ k\0324\ 1 + el + (* - (* + : 21 line for a\032E\ X\032E\ i\0324\ el ≠ d\0324\ a\032E\ X\032E\ i\0324\ el - : h\0324\ X\032E\ i\0324\ 1 + el X\032E\ i\0324\ el - : t\0324\ d\0324\ h\0324\ / : 13 mat s\0324\ V\032E\ i\0324\ 1 + el V\032E\ i\0324\ el - (t\0324\ 4 / * V\032E\ i\0324\ el + (t\0324\ * V\032E\ i\0324\ el V\032E\ i\0324\ 1 + el 2 / + (- : I\0324\ I\0324\ s\0324\ d\0324\ h\0324\ * 6 / * Y\032E\ i\0324\ 1 + el Y\032E\ i\0324\ el - 2 / t\0324\ * - Y\032E\ i\0324\ el + (d\0324\ * - : 31 line I\0324\ if b\032E\ X\032E\ j\0324\ 1 + el ≠ d\0324\ X\032E\ j\0324\ 1 + el b\032E\ - : h\0324\ X\032E\ j\0324\ 1 + el X\032E\ j\0324\ el - : t\0324\ d\0324\ h\0324\ / : 13 mat s\0324\ V\032E\ j\0324\ el V\032E\ j\0324\ 1 + el - (t\0324\ 4 / * V\032E\ j\0324\ 1 + el + (t\0324\ * V\032E\ j\0324\ 1 + el V\032E\ j\0324\ el 2 / + (- : I\0324\ I\0324\ s\0324\ d\0324\ h\0324\ * 6 / * Y\032E\ j\0324\ 1 + el Y\032E\ j\0324\ el - 2 / t\0324\ * - Y\032E\ j\0324\ 1 + el + (d\0324\ * - : 31 line I\0324\ if 51 line if if 11 line :

Cubic Spline interpolation vector of coefficients. X\032E\ Y\032E\ CS.V n\0324\ X\032E\ length : c\0324\ 1 el X\032E\ 3 el X\032E\ 1 el - X\032E\ 3 el X\032E\ 2 el - / - : d\0324\ 1 el 0 : b\0324\ X\032E\ 3 el X\032E\ 1 el - 3 / c\0324\ 1 el X\032E\ 2 el X\032E\ 1 el - 6 / * - : c\0324\ 2 el X\032E\ 2 el X\032E\ 1 el - 6 b\0324\ * / X\032E\ 3 el X\032E\ 2 el - X\032E\ 2 el X\032E\ 1 el - / X\032E\ 2 el X\032E\ 1 el - X\032E\ 3 el X\032E\ 2 el - / - (* : d\0324\ 2 el 1 b\0324\ / Y\032E\ 3 el Y\032E\ 2 el - X\032E\ 3 el X\032E\ 2 el - / Y\032E\ 2 el Y\032E\ 1 el - X\032E\ 2 el X\032E\ 1 el - / - (* : n\0324\ 3 > k\0324\ 3 n\0324\ 1 - (range b\0324\ X\032E\ k\0324\ 1 + el X\032E\ k\0324\ 1 - el - 3 / c\0324\ k\0324\ 1 - el X\032E\ k\0324\ el X\032E\ k\0324\ 1 - el - 6 / * - : c\0324\ k\0324\ el X\032E\ k\0324\ 1 + el X\032E\ k\0324\ el - 6 b\0324\ * / : d\0324\ k\0324\ el 1 b\0324\ / Y\032E\ k\0324\ 1 + el Y\032E\ k\0324\ el - X\032E\ k\0324\ 1 + el X\032E\ k\0324\ el - / Y\032E\ k\0324\ el Y\032E\ k\0324\ 1 - el - X\032E\ k\0324\ el X\032E\ k\0324\ 1 - el - / - d\0324\ k\0324\ 1 - el X\032E\ k\0324\ el X\032E\ k\0324\ 1 - el - 6 / * - (* : 31 line for if a\0324\ X\032E\ n\0324\ el X\032E\ n\0324\ 2 - el - (- X\032E\ n\0324\ el X\032E\ n\0324\ 1 - el - (c\0324\ n\0324\ 2 - el * - : b\0324\ X\032E\ n\0324\ 1 - el X\032E\ n\0324\ 2 - el - c\0324\ n\0324\ 1 - el a\0324\ * - : d\0324\ n\0324\ el X\032E\ n\0324\ el X\032E\ n\0324\ 1 - el - (d\0324\ n\0324\ 2 - el * - d\0324\ n\0324\ 1 - el a\0324\ * - b\0324\ / : k\0324\ n\0324\ 1 - (1 range d\0324\ k\0324\ el d\0324\ k\0324\ el c\0324\ k\0324\ el d\0324\ k\0324\ 1 + el * - : for d\0324\ 1 el d\0324\ 1 el X\032E\ 2 el X\032E\ 1 el - X\032E\ 3 el X\032E\ 2 el - / d\0324\ 3 el * - : X\032E\ Y\032E\ d\0324\ 032 stack 14 mat eval 131 line :

Adaptive Cubic Spline interpolation. 1 f a b CS.A x g x i x f * + : Z x g a b APlot : z Z 1 el : k 2 Z length range Z k 1 - el Re Z k el Re ≡ k Z k el 12 mat trace z z Z k el stack eval : if for z Re z Im CS.V eval 51 line :

Setup

MaxIters 100 :

APlotAngle 15 ° * :

APlotDist 10 % * :

APlotNG 4 :

For the cloud:

a b c b a ≤ (a c ≤ (& 11 line : lele

a b c ltlt b a < (a c < (& 11 line :

a b c b a ≤ (a c < (& 11 line : lelt

t0 1 time :

CS Interpolation

ε 1015 -^: X ε 12 0.25 range : x f 1 x ln - sin : M X X f vectorize CS.V : x p M x CS : 51 sys

M X X f vectorize CS.V : x p M x CS : 21 sys

data M 1 el M 2 el augment : plotPts data o 5 black plot : 21 line

collect/specify/plot all 'M' points

plotPts x f 2 1 sys

a b 12 mat 1 9.1 12 mat :

a b 12 mat X min X max 12 mat :

a f 0.8415

a p 0.8415

b f 0.935 -

b p 0.935 -

a f a diff 0.5403 -

M a CS.D 0.545 -

b f b diff 0.039 -

M b CS.D 0.039 -

a f a diff a diff 0.3012 -

M a CS.D2 0.2039 -

b f b diff b diff 0.0156

M b CS.D2 0.0156

x f x a b int 3.3316 -

M a b CS.I 3.3316 -

Alvaro

x p x diff x f x p -

1 time t0 - 26.3 s *

t0 1 time :

$ε 1015 -^: X ε 12 0.25 range : x f 1 x ln - sin : Y X f vectorize : x q X Y x cinterp : 51 sys$ $x p x ε e int #$ x f x q - $1 time t0 - 0.7 s *$ $x q x ε e int 1.3761$ M x CS.D x f 2 1 sys x f x q - x f x p - 2 1 sys M x CS.D2