Utilities ...Vect(var,n)

a8938ce0-1b9f-402b-b955-97fee0e2785d 86 4 5 decimal

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

Label nodes on graph ... User project

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

Vect(size,n) ... Vect(color,n) var n Vect n 1 ≡ var ct 1 : j 1 n range a ct 1 el var : ct ct 1 + : 21 line for v a 1 el : k 2 a rows range v v a k el stack : for v 51 line if :

Vect(sizechar,n)... Vect(colorchar,n) size n Vect n 1 ≡ size ct 1 : j 1 n range a ct 1 el size : ct ct 1 + : 21 line for v a 1 el : k 2 a rows range v v a k el stack : for v 51 line if :

x h λ CONVOL simple convolution j 2 x rows range c j 1 - el x j el λ h j 1 - el (* * : for x x 1 el c stack augment 31 line :

x 012345678910111 mat :

h 1.23456789 2.3456789 3.456789 4.56789 5.6789 6.789 7.89 8.9 101112111 mat :

Cost of polluants . . . user style

1. polluants [x] are integers .

2. cost [h] is user per application .

]]>

3 appVersion 0.98.6179 ≡

cost factor λ 1 :

Ω Clear 1

Graph Label system polluants x h λ CONVOL : nodes x : NodesLabel nodes var2str vectorize : X x : Y polluants 2 col : 12 mat n polluants rows : 51 line

Convolution infinitesimal sequence Graph X Y augment : x f X Y x ainterp : x c 10 x ≤ (x 11 ≤ (& cases : data Graph . 14 black plot : 41 sys

plot symbol converted2Numbers

size

symbol color

nodes_label Graph NodesLabel 10 n Vect red n Vect augment :

User model cumulative cost of polluants

x f x c * data 1 25 1 25 - 2 2 mat nodes_label 4 1 sys

Graph 001 1.23456789 2 4.6913578 3 10.370367 4 18.27156 5 28.3945 6 40.734 7 55.23 8 71.2 99010110112 mat

plot arguments x f x c * data 125125 - 22 mat nodes_label 41 sys

Y sum 430.1

Inventory cost

j 12 range :

Y j el sum 1.23

j 36 range :

Y j el sum 61.7

data x 010 0.1 range : j 1 x rows range y j el x j el f eval : for x y augment 31 line :

View Graph . 14 black plot :

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 :

Σ data Σ :

j 71 :

x f x 010 int 373.42

X j : Cost Σ j row : 12 mat

track 02565256 Cost 11 el Cost 12 el 32 mat :

data Σ track nodes_label View 5 1 sys

cumulative cost @ day 'j' Cost 7 130.07 12 mat

x y Clear 1

'maple'... apetizer

x y f0 11 x + (3^+ 3 y * 1 x + (* - y 3^+ abs ln : x y f1 11 x + (3^+ 3 y * 1 x + (* - y 3^+ abs ln 0.5 + : x y axes x y / - : 31 line

f0 f1 axes 3 1 sys

y F x y f0 0 ≡ x solve maple :

Domain of the triplet u 1.587401051968190 :

x f3 x F 3 el Re 0 x ≤ (x u ≤ (& ignore imag cases :

x f2 x F 2 el Re 0 x ≤ (x u ≤ (& ignore imag cases :

φ0(u)-φ2(u)... 2.06083896700751*10^{-7} 1.587401051968186 1.587401051968195 21 sys

'u'is like the Goat radius constant i.e: can NOT be solved exactly !

x φ0 x F 1 el Re : x φ1 x f2 : x φ2 x f3 : x φ3 x F 4 el : 41 sys

u 2 / 0.7937

u φ0 0.2599

plot x u 2 / + φ0 1 + x u 2 / + φ1 1 + x u 2 / + φ2 1 + x u 2 / + φ3 1 + u 2 / u φ0 1 + + 20 red 15 mat u 2 / - 0 + 100 black 15 mat 61 sys :

Smath QuickPlot complete looping ... ignore Maple sol [5, 6]

plot

`

Goat radius constant

... a more mathematical alternative solution R L Goat 212 / R 2^* 2 L 2 R * / asin * 2 L 2 R * / asin * sin - (* (* L 2^L 2 R * / acos * + :

R 1 :

L 1.158728473 :

the red  XY coordinates. system 0 L L π R 2^* 2 / π R 2^* 2 / 023 mat transpose :

x c 10 x ≤ (cases :

R x Goat x c * system 2 R * π R 2 ^ * + 40 black 1 5 mat 3 1 sys

elegant mathematical alternate well posed solution
... enables plot the goat function. R L Goat π R 2^* 2 / - 0 ≡ L 1 1.2 solve 1.158728473

π R 2^* 2 /

L

ASHREA 'approxim Чебышев' distribution

'N' even ... increase for smaller centroid-ring area D 1 : N 24 : 12 mat

x_D i 1 N 2 / range su i 1 + el D 2 / D N 2 i * - 4 N * / sqrt * - : for j N 2 / 1 + (N range sv j N 2 / - el D 2 / D 2 j * N - 4 N * / sqrt * + : for su sv stack 31 line :

data n Decimate collect every 'n' pairs u 0 : v 0 : 12 mat range 1 data rows n range : k 1 : i range u k el data i 1 el : v k el data i 2 el : k k 1 + : 31 line for u v augment 31 line 41 line :

x_D x_D eval : j 1 x_D rows range : XY j x_D augment : sub XY 2 XY rows range 12 range el : 41 sys

sanity ASHREA source sanity sub 3 Decimate : 2 0.021 4 0.067 6 0.118 8 0.177 10 0.25 12 0.356 14 0.644 16 0.75 18 0.823 20 0.882 22 0.933 24 0.979 122 mat

t U D ashrea U 1 el U 2 el 21 mat * :

ashrea 022 - 022 mat :

'N' governs
the sub 'Ø'
blue circle N 200 : span 1 N π * / : α π : (overlap ≡ 31 sys

t U D 2 U 2 el * 2 U 1 el * - 21 mat

p 0.5 0.644 0.75 0.823 0.882 0.933 0.979 71 mat :

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

Nrings Ω 00013 mat : p p transpose W p span 21 mat : ODE W 0 α N t U D dn_GearsBDF : Ω Ω ODE stack : 31 line for out Ω 3 col Ω 2 col augment : 31 line :

Φ Nrings Nrings rows N - (Nrings rows range 12 range el :

jCheby 1 p rows range : Cheby p jCheby 0 * augment : C Cheby . 14 black plot : Φ Φ . 8 blue plot : 41 sys

Disposition of the ASHREA of 'Unary-rings' ≡ Area '0.979' is a bit tangent2duct wall ... check Annubar

Φ Nrings C 1 0 + 200 black 1 5 mat 4 1 sys

Φ Nrings C 10 + 200 black 15 mat 41 sys

Delicious X_Y spline plot [2D Isometric.sm]

image

Stem(data)... plot(data,char,size,clr)

data char size clr plot k 1 data rows range r3 k el char : r4 k el size : r5 k el clr : 13 mat for data r3 r4 r5 augment 21 line :

Spectrum stem plot for QuickPlot unicolor data Stem stemX 1 data rows range : i 1 stemX rows range stem i el data i 1 el data i 2 el data i 1 el 022 mat : for i 1 stem rows 1 - range stem i el stem i el stem i 1 + (1 el 1 el 012 mat stack : for V stem 1 el : j 2 stem rows 1 - range V V stem j el stack : for V V data data rows row stack : 71 line :

Animate a system ...

1. Set canvas before installing the plot argument 2. Tansit/isolate via file ...

timing 3 min ≡

Laplace linear diffusion x t f t diff x t f x 2 diff ≡

F α 0.25 : β 20 : 13 mat f => spatial initial conditions i 141 range f 1 i el 0 : f 1 β el 1 : 21 line for recursive finite differences 'x' sol's t 150 range x 240 range f t 1 + x el f t x el (α f t x 1 - el 2 f t x el * - f t x 1 + el + (* + : for for f 61 line :

isolate/accelerate F F eval :

CWD C:\SmathFile CurrentDirectory :

F CWD Diffusion wfile 1

D CWD Diffusion rfile :

step explorer j 15 :

t Diffusion j 141 range : j D t row transpose augment 21 line :

Highlight a specific diffusion step s j Diffusion . 12 black plot :

τ 1102 range :

ω t Diffusion . 12 black plot :

1 Diffusion Stem 2 Diffusion Stem j Diffusion Stem s 4 1 sys

Ω D 20 col :

Ω 1 0.5 0.375 0.312 0.273 0.246 0.226 0.209 0.196 0.185 101 mat

t 1 ≡ t diffusion t Diffusion Stem : t diffusion if 1 Diffusion Stem t 1 ≡ t diffusion t Diffusion Stem : t diffusion ω 2 1 line if 16 0.185 24 0.185 2 2 mat 4 1 sys

Ω Ω rows row 0.0796 11 mat

Laplace Diffusion Anim system t 1 ≡ t diffusion t Diffusion Stem : t diffusion if 1 Diffusion Stem t 1 ≡ t diffusion t Diffusion Stem : t diffusion ω 21 line if 16 0.185 24 0.185 22 mat 41 sys

M Lo Hi Scale Hi Lo - (M M Min - M Max M Min - / (* Lo + :

image D 01024 Scale :

Inverse 512 image - :

iVBORw0KGgoAAAANSUhEUgAAACkAAAAzCAIAAACFT7ZcAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsMAAA7DAcdvqGQAAAWsSURBVFhH7de5T5ZbEAZwQURxRRY3FBUQg+ICkbihsUDcwR2NKyKoSIwxGBckITHRRBMtbKSAChpjZewMrQmthTWdrX8Axf0l58ubk6u5Xj5ebnVPMZk887znYb4zM+cw6//1r9eBAwfYgYGBxP+PVmFh4fXr13/+/EmV5UMysRldeXl569atu3///sTExIoVK1g+BJ5hzNDKycmZP3/+li1bHj9+/O3bt8nJSZYPgYtmeDOxZs+e7efduXNnX1/f169ff/z4wfIhcNEMbybWnDlzli1btn///v7+/i9fvnz//p3lQ+CiGV7qy086b9681atXNzU10fv06dP4+DjLh8BFZ+pn95MuWLBg/fr1hw8ffvbs2YcPH8bGxlg+BC46Iz+7hFTykiVLNmzYcOzYMWc8MjLy+fNnlg+Bi+Kkn7od8/Pzi4qKampqWlpanjx5MjQ09PHjR5YPgYvipK+dm5vrOEtLS2tra1tbWx89evT+/fvR0VGWD4GL4mBmvkllScVB6uDly5dv27bt1KlTDx8+fPfu3fDwMMuHwEVxMNNM3V4OUimtWrWqvr6eUm9v75s3bwYHB1k+BC6Kk/KR20vvLlq0SCPt2LHj9OnTDx48ePXqlaRZPgQuioOZprYjVETKuLy8vKGhgZIZ/uLFi7dv37J8CFwUBzO1I5eEI5w7d+7SpUvdGcYnpXv37j1//vz169csHwIXxcFM7ciDtgLWQhUVFbt27aLU09Pj8n758iXLh8BFcTDT1FY+CrikpKSqqmr37t2U7t69a5RKmuVD4KI4mKmVm8NTPgpY+xpee/bsoXTnzp2nT59KmuVD4KI4mPjpHHnQXrhwoZuqurp67969lG7fvm2kSJrlQ+CiOJjpaIfDVrqaxytl48aNQburq8tIkTTLD9qiOJj4KRx50Fa6ixcvXrlypaHd2NhoknR2dmprSbN8CFwUBzOdUg/aSlfjGluJdkdHh7aWNMtPtHEw0yl1x6Zo7eVJVFZWtmnTpn379p08efLGjRvaWtIsHwIXxcHE99W0jtwfHrQLCgoMDSNz8+bNQdu7WFtLmuUHbVEcTPygnX3qQVvRallDY82aNS5KGm7Ma9eudXd3S5rlQ+CiOJj4odSnpe3YFK2WLS4uNq69hYP21atXtZakWX7QFsXBxJ9uqf9Ne+3atXb3JPVKuXLlyq1btyTN8iFwUZw0tTWMcWFYuiq2bt0atC9fvqy1TDSWH7RFcTDxp9tmodDsYlwYlon2iRMnLl26dPPmTUmzfEiijYnvq+xL3Z8ctDWMvQxLr+BE++LFi1pL0iw/0cbBxE/aLJvUY22jyo6uSI8yGsePH79w4UJ7e7ukWT4ELoqDiZ+CtlbRrEaVd2Cs3dbWprUkzfJjbUx8X2XfZr5RLIm2S6KysnL79u16yX8C58+f11qSZvkQuCgOZqKdZbkl2gaFMRlrHz169Ny5c1pL0iwfEmvjh/EyLW1tGrRdUJ4lsbbWkjQba+NgBu0sW9wHjipoGxRGdKx95MiRs2fPai1Js3xIrI2fjJcpH3nQVqja1C5GtMsx1j5z5ozWkjTLj7Ux8X2VtHj22oZU0PYcq6uro+E/Xg8VrSVplg+Bi+IE7TDaUtA2ol3MsbbHgtaSNMuPtTHxU9AOQy3W9j45dOiQC1trSZrlQ+CxdvajDVuZaJKg7XrwKPAMjbWVt6TZWBsHEz9oZ9NmibYR8at2c3OzC1t5S5rlQ37Vzma8oMbahrO9PEhibZcmVUmz/FgbE99XsbaV2f2fF55D8o0GDdquRTt6fgftgwcPuri0FmGWDwnaOJj4QXvKLR5rG0+Gc6ztf4CgrbVos0EbHmv7KhltU9ZWoom2a9FDLNZ2cWktSbP8WBsTP9GeWqn/UbupqcnFpbwlzfIhKWsbDkaji+G32spb0uxvtX3157E6a9Zfczw2Si9fRQkAAAAASUVORK5CYII= image

J 0 Ω rows 1 - (range :

plot J Ω augment J - Ω augment 21 sys :

plot iVBORw0KGgoAAAANSUhEUgAAACkAAAAzCAIAAACFT7ZcAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsMAAA7DAcdvqGQAAAYgSURBVFhH7de3b1VLEAZwF/wBtHS0NJQUlIgOUSEaqEwBDQVIWAIJCYtgIRACRM5R5BwNbshR5CAwmJxztAET3k/scrQPX865xtevelOMRt9+Z7872pnZvVX/W9l24MABvra2Nov/I3vz5s3y5cu7du1KlRdD4lqn2tevX+/cuTNjxozu3bs/efKEF0PgkdFJ9uPHj+bm5kuXLk2ePLlnz55dunThxRC41ciruNn627dvr1+/Pnny5MSJE3v37t2tWzdeDIFb7RR5m7IvX748ffr04MGD48eP79u3b48ePXgxBG410OI3FbGw4/fv31taWu7fv9/Q0ECvf//+vXr14sUQuFWcQI5fdtDCXkw1ffjw4datW3v37p0wYcLAgQP79OnDiyFwqziR3XH5uM3PpFtbW/VSY2Pjrl27nPHgwYP79evHiyFwqzhZ6izu8hcWN/hpdvz8+fPLly+vXr26ffv2urq6IUOGDBgwgBdD4FZxUm0W92qXxU9/mRp2nM+ePdNO27Zt01fDhg0bNGgQL4bAreKEak8t7limxY9+mVQc5MePHw2T8+fPb9myZcqUKcOHD6+urubFELhVHMzfUmdx30KL9MTs5SCV0sOHD8+cOUNp6tSpI0eOHDp0KC+GwK3i/HbkmcXdcywS/2320rvv3r3TSKdPn968efO0adNqamokzYshcKs4mCW1WdQoaZHSxhyhIlLGd+/ePXXqFKXp06ePGTNmxIgRvBgCt4qD2fbIM4tKv1lcbGOScISfPn169erV7du3T5w4QWnmzJljx44dNWoUL4bAreJgljzyzKJeZhEuZUFbAWuhpqam48ePU5o1a9a4ceNGjx7NiyFwqziY+dosqgaLWCmzi/JRwM+fP79x48axY8cozZ492yiVNC+GwK3iYP6p3FKLwiwCpczhKR8FrH0Nr6NHj1KaO3fupEmTJM2LIXCrOJj4OUeeWbna79+/d1Ndv379yJEjlObNm2ekSJoXQ+BWcTArox0OW+lqnsePH1+7di1oL1iwwEiRNC8O2lZxMPELj5yVpa103759++jRI0P78OHDJsnChQu1taR5MQRuFQezsNSDlaWtdDWusZVpL168WFtLmhdn2jiY5ZQ6K9B2bIrWXp5EDx48uHLlyqFDh7Zu3bpkyRJtLWleDIFbxcHE91Xhkedp++FB2zvQ0DAyL1++HLSXLVumrSXNi4O2VRxM/KCdn3qxtqLVsobGvXv3XJQ03Jje5HPmzJE0L4bAreJg4odS75C2Y1O0WvbFixfG9cWLF4P2ihUrtJakeXHQtoqDiV9OqbdD2x8Au3uSeqWsXLly/vz5kubFELhVnEpqaxjjwrB0VVy4cCFor1q1SmuZaLw4aFvFwcQvp83ytEOh2cW4MCwz7R07dqxevXrRokWS5sWQTBsT31eFpf5HbT85aGsYexmWXsGZ9po1a7SWpHlxpo2DiZ+1WU7qZWkbVXZ0RXqU0di5c+fatWuXLl0qaV4MgVvFwcSvgLZW0axGlXdgqr1u3TqtJWlenGpj4vuqsM3ytBVLpu2SuHnz5rlz5/SSfwLr16/XWpLmxRC4VRzMTDu/3Iq1DQpjMtXevXv3hg0btJakeTEk1cYP46VD2to0aLugPEtSba0laT7VxsEM2oUtXlrbB44qaBsURnSqvWfPno0bN2otSfNiSKqNn42XnCPP01ao2tQuRrTLMdXetGmT1pI0L061MfF9lbX432sbUkHbc+zs2bM0/OP1UNFakubFELhVnKAdRlsFtI1oF3Oq7bGgtSTNi1NtTPwKaIehlmp7n9TX17uwtZakeTEEnmqXM9r+qK1MNEnQdj14FHiGptrKW9J8qo2DiR+089usQNuIaKu9b98+F7byljQvhrTVLhwvJbRRU23D2V4eJKm2S5OqpHlxqo2J76tUm8XdEyut7ZB8o0GDtmvRjp7fQXv//v0uLq1FmBdDgjYOJn7Qzm/xAm3jyXBOtf0HCNpaizYftOGptq+y0dZubSWaabsWPcRSbReX1pI0L061MfEz7ZxS/xvthoYGF5fyljQvhlRY23AwGl0MJbWVt6T5ktq+KhyrVVVV/wAfFH5vZNvUkAAAAABJRU5ErkJggg== Inverse

'though roots' on running canvas .. Inst_Titration SOLVE.sm

1. concentration strong base [mol/L]
2. concentration strong acid [mol/L]
.....
3. ion product of water [mol/L]
4. Volume [mL] of acid solution being titrated cNaOH 0.05 : cHCl 0.1 : Kw 1014 -^: vHCl 50 : 41 sys

roots solve

Titration curve before Equivalence point x ≤ [100 mL] pH1 x 11002 range : i 1 x rows range sol i el H Kw H / - vHCl cHCl * x i el vHCl + / - x i el cNaOH * x i el vHCl + / + 0 ≡ H 1 roots : for sol sol rows el 1107 -^* : x sol log10 vectorize - augment 41 line :

Titration curve beyond Equivalence point  x ≥ [100 mL] pH2 u 100.001 200101 range : i 1 u rows range sol i el OH Kw OH / - u i el cNaOH * u i el vHCl + / - vHCl cHCl * u i el vHCl + / + 0 ≡ OH 100 roots : for u 14 sol log10 vectorize + augment 31 line :

Titration curves ... from 300 roots canvas

pH1 pH2 2 1 sys

1007 + 10 black 15 mat

x pH1 H Kw H / - vHCl cHCl * x vHCl + / - x cNaOH * x vHCl + / + 0 ≡ H 1 roots : u pH2 OH Kw OH / - u cNaOH * u vHCl + / - vHCl cHCl * u vHCl + / + 0 ≡ OH 100 roots : 21 sys

1 x ≤ (x 100 ≤ (&

100 u ≤ (u 200 ≤ (&

<= pH 48 pH1 log10 - 14128 pH2 log10 + 21 mat 1.58 11.90 21 mat

Convivial roots on running canvas ...

γ 1.4 : (Air isentropic ≡ k 0.88128485 : (Pitot constant ≡ p static pressure ≡ qc dynamic pressure ≡ x Pitot dynamic ratio qc/P ≡ x qc p / ≡ (supersonic Pitot ratio ≡ 61 line

Observe the canvas solutions:

as 'x' in Supersonic(x) runs on the plot canvas,

the 'roots' solves and plots each canvas pixel result ...

each result on the basis of 96 ppi [pixel per inch].

]]>

Plot limit x c 10 x ≤ (x 5 ≤ (& cases :

t0 1 time :

x Subsonic 2 γ 1 - / (x 1 + (γ 1 - γ /^1 - (* sqrt 0 x ≤ (x k ≤ (& Supersonic cases : x Supersonic U x 1 + (117 U 2^* / - (2.5^* sqrt / k - 0 ≡ U 1 roots k x ≤ Subsonic cases : 21 line

x Subsonic x Supersonic x c * k 1 + 20 black 1 5 mat 3 1 sys

1 time t0 - 69 s *

supersonic

5 Supersonic 2.07

subsonic

2D surface spline matrix ...

Spline(M,x,y)

M x y Spline A 0 : B 0 : C 0 : 13 mat A M 112 M cols submatrix : B M 2 M rows 11 submatrix : C M 2 M rows 2 M cols submatrix : j 1 A length 2 range CC j el B C 1 C rows j j submatrix y linterp : for A transpose CC x linterp 61 line :

data x/y 0.05 0.1 0.2 0.5 1245102030451 1.8 1.7 1.65 1.55 1.45 1.4 1.35 1.35 01.3 1.3 1.3 1.3 0.9 1.85 1.75 1.7 1.6 1.55 1.5 1.45 1.45 1.4 1.4 1.4 1.4 0.7 1.9 1.85 1.75 1.7 1.65 1.6 1.6 1.55 1.55 1.55 1.55 1.55 0.5 1.95 1.85 1.8 1.75 1.75 1.7 1.7 1.7 1.7 1.7 1.7 1.75 0.3 1.95 1.9 1.85 1.8 1.8 1.8 1.8 1.8 1.8 1.85 1.85 1.85 0.2 1.95 1.9 1.85 1.8 1.8 1.8 1.8 1.8 1.85 1.9 1.9 1.95 0.16 1.95 1.9 1.85 1.85 1.8 1.8 1.85 1.85 1.85 1.9 1.9 1.95 813 mat :

sanity
sanity
user data 1 0.5 Spline data 20 0.5 Spline data 0.75 0.45 Spline 31 mat 1.75 1.7 1.76 31 mat

mesh the ► 'x' direction at will mesh the ▼ 'y' direction at will

1. Construct the 3D plot

Clear the system u v Clear 1

<= Introduce the 3D surface plot u v q u v data u v Spline 31 mat :

CreateMesh paramaters L 1 - : H 1 : N 31 : 31 sys

N mesh within the bounds [L, H] in each directions M u v q L H L H N N CreateMesh :

Euler system ∠(x)

Roll ∠ Ω1 1000 ∠ cos ∠ sin - 0 ∠ sin ∠ cos 33 mat :

∠

Pitch ∠ Ω2 ∠ cos 0 ∠ sin 010 ∠ sin - 0 ∠ cos 33 mat :

Yaw ∠ Ω3 ∠ cos ∠ sin - 0 ∠ sin ∠ cos 000133 mat :

user angular display α π 8 / : β 0 π * : φ 0 π * : 13 mat

<= Rotate around Y θ 0.75 :

Euler combined γ α Ω1 β Ω2 * φ Ω3 * eval :

2. 2D interpolated surface plot

M f B M θ Ω3 * γ * eval : B 1 col B 2 col augment 21 line :

x y axes x y / - :

plot the 2D surface plot M f :

plot axes 2 1 sys

°C vx 12345616 mat :

мм vy 1234551 mat :

Вт/м else 49172535454101928395051120314254512233547606142638516656 mat :

vx vy else 3.5 4.5 InterpBilinear vx vy else 53 InterpBilinear 21 mat 30.5 4221 mat

Prepare
SplineMatrix U vy else augment (: M vx augment U stack : 21 sys

The spline matrix M 1234561491725354524101928395035112031425445122335476056142638516667 mat