d64f37d1-2e32-459c-8d7d-32d7e98d91fe 51 4 5 decimal

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

vx vy char size color plotG n vx length : plot vx 1 el vy 1 el char size color augment : i 2 n range plot plot vx i el vy i el char size color augment stack : for plot 41 line :

http://www.chemicalaid.com/tools/reactionstoichiometry.php?equation=S+%2B+O2+%3D+SO3

z n RepEnd m z rows : i 1 n range w i el z 1 el : for i n 1 + (m n + (range w i el z i n - el : for i m n + 1 + (m 2 n * + (range w i el z m el : for w 51 line :

This Smath document reproduces Mathcad 11 De Boor spline. Takes much longer to compute. De Boor reproduces De Casteljau more compact algorithm. Bézier-Bernstein were "Alma Mater" companions. Bézier and De Casteljau did not know each other until the founder of ADOBE met Bézier. The idea of De Boor was to code De Casteljau algorithm to be executable on the available computing machinery of his time. A nice lagacy.

X 1.5 35745 0.5 - 71 mat :

Y 010424671 mat :

r 1 :

n 2 :

X X r RepEnd :

Y Y r RepEnd :

m 1 X rows n 2 + (+ range :

# # range

XY 1.5 0 1.5 03150744254 0.5 - 6 0.5 - 692 mat :

X XY 1 col :

Y XY 2 col :

Cntrl X Y . 14 black plotG :

Colibri data for De Boor

XY Cntrl 2 1 sys

u i 1 m rows range u i el m i el : for u 21 line :

m u :

v m N0 two errors code u m el v ≤ u m 1 + el < 10 if 21 line :

Smath does not support an index in the calling function [LHS], neither it supports the calling function being indexed.

1.75 1 N0 #

v j N0 μ j 11 line : u μ el v ≤ u μ 1 + el < 10 if 21 line :

The index must be locally assigned in the EXECUTABLE procedure

μ 2 :

u μ el 2

u μ 1 + el 3

5.5 3 N0 1

Not yet instructed: Smath is rigid wrt the bolean. The above Bolean is underdetermined, Ah !

v j N0 μ j 11 line : u μ el v ≤ (v u μ 1 + el < (& 10 if 21 line :

Now, the instruction is complete and executable for De Boor spline.

2.5 1 N0 0

11 N0 1

1.5 1 N0 1

21 N0 0

v j N1 μ j 11 line : v u μ el - u μ 1 + el u μ el - / v μ N0 * u μ 2 + el v - u μ 2 + el u μ 1 + el - / v μ 1 + N0 * + 21 line :

v j N2 μ j 11 line : v u μ el - u μ 2 + el u μ el - / v μ N1 * u μ 3 + el v - u μ 3 + el u μ 1 + el - / v μ 1 + N1 * + 21 line :

μ j 11 line : v μ N1 21 line μ j 11 line : v μ N2 21 line μ j 11 line : v μ N3 21 line 31 mat

v j N3 μ j 11 line : v u μ el - u μ 3 + el u μ el - / v μ N2 * u μ 4 + el v - u μ 4 + el u μ 1 + el - / v μ 1 + N2 * + 21 line :

Not needed for n = 2 then not collected

v j N μ j 11 line : v μ N1 21 line μ j 11 line : v μ N2 21 line 21 mat :

Collect the respective De boor

n 2

<= Spline order set from above

v Sx v k N n el X k el * (k 1 X rows sum :

v Sy v k N n el Y k el * (k 1 X rows sum :

vx 310 3.05 range :

<= populate for mor points at will

i 1 vx rows range SX i el vx i el Sx : for

i 1 vx rows range SY i el vx i el Sy : for

DeBoor SX SY augment :

De Boor parametric order 2 wrt Smath ... kind of Convex Hull

XY DeBoor Cntrl 3 1 sys

De Boor components

vx SX augment vx SY augment 2 1 sys

jmGiraud: 2015/09/27

DeBoor Σ

DeBoor Deriv

data Σ Finite differences Cumul Integration vx data 1 col : vy data 2 col : Y 1 el 0 : i 2 vx rows range Δx i el vx i el vx i 1 - el - : Y i el Y i 1 - el 12 / vy i 1 - el vy i el + (Δx i el * (* + : 21 line for vx Y augment 61 line :

Cumulative Integral

DeBoor Σ

data Deriv data Derivative BackSubstitution X data 1 col : Y data 2 col : n X rows : deriv 1 el 0 : deriv n el 0 : i 2 n 1 - range deriv i el 0.5 Y i 1 - el Y i el - X i 1 - el X i el - / Y i el Y i 1 - el - X i el X i 1 - el - / + (* : for deriv 1 el deriv 2 el : deriv n el deriv deriv rows el : X deriv augment 101 line :

Finite differences derivative

DeBoor Deriv