<content> Utilities </content>

3bb344d3-ebf8-49aa-8ce2-dad7caf0b4b9 445 4 5 false false false 0 3 decimal

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

NaN 0 :

n\032E\ Id n\032E\ identity :

Τ ΤΤ 0 time Τ 1 el - 11 line :

x\032E\ TryNaN x\032E\ 0 try x\032E\ NaN if 11 line :

Τ Τ Τ IsDefined told\0324\ Τ 2 el : Τ 2 el 0 time : 0 time told\0324\ - 31 line Τ 0 time 0 time 12 mat : Starts 21 line if 11 line :

x\032E\ IsMatrix x\032E\ rows 0 ≥ TryNaN NaN ≠ eval 11 line :

Τ Τ Starts

f\032E\ x\032E\ y\032E\ feval try(f\032E\ num2str (x\032E\ num2str, y\032E\ num2str),'NaN) concat str2num 11 line :

a\032E\ b\032E\ n\032E\ pRange a\032E\ b\032E\ a\032E\ - n\032E\ / 0 n\032E\ range * + eval 11 line :

B\032E\ N\032E\ pRange B\032E\ 1 el B\032E\ 2 el N\032E\ 1 el pRange B\032E\ 3 el B\032E\ 4 el N\032E\ 2 el pRange 12 mat eval 11 line :

x\032E\ y\032E\ z\032E\ pAxis 000 x\032E\ 000000 y\032E\ 000000 z\032E\ 00073 mat :

M iFind M M : n 1 M length range : 12 mat M 0 ≡ (vectorize n el n augment eval 11 findrows 2 col eval 0 try 21 line :

M\032E\ y\032E\ pAxis M\032E\ cols 3 ≡ (y\032E\ cols 2 ≡ (& A\0324\ 0 : k\0324\ 13 range : k\0324\ m\0324\ 00013 mat 1.05 M\032E\ 1 col max * 3 identity k\0324\ row * stack y\032E\ * : 13 mat A\0324\ k\0324\ el polyline k\0324\ m\0324\ 11 mat red green blue 31 mat k\0324\ el solid 2 stack : eval 21 line todo if 11 line :

<content> Color utilities </content>

A\032E\ pCycleC B\0324\ A\032E\ sys2mat eval : M\0324\ B\0324\ 1 el 11 mat : 12 mat k 2 : k B\0324\ length ≤ k k 1 + : n 1 pCycleC range M\0324\ M\0324\ rows 1 + el ∞ ∞ 12 mat : for M\0324\ M\0324\ rows 1 + el B\0324\ k el : 21 line for M\0324\ mat2sys.1 eval 31 line :

r\032E\ g\032E\ b\032E\ a\032E\ rgb hex\0324\ 0123456789ABCDEF : d\0324\ 0 : 12 mat c\0324\ H\0324\ n\0324\ 0255 c\0324\ 0 round 12 mat min 12 mat max : m\032E\ n\0324\ 16 mod : 12 mat n\0324\ m\032E\ - 16 mod 0 ≡ hex\0324\ n\0324\ m\032E\ - 16 / 1 + 1 substr hex\0324\ m\032E\ 1 + 1 substr concat hex\0324\ m\032E\ 1 + 1 substr if 21 line : h\0324\ 255 a\032E\ * H\0324\ r\032E\ H\0324\ g\032E\ H\0324\ b\032E\ H\0324\ concat : n\0324\ h\0324\ strlen : 12 mat k\0324\ 1 n\0324\ range d\0324\ d\0324\ hex\0324\ h\0324\ k\0324\ 1 substr findstr 1 el 1 - (16 n\0324\ k\0324\ -^* + : for d\0324\ 232^2 d\0324\ * 232^1 - ≥ (* - eval 51 line :

C\032E\ rgb C\032E\ 1 el C\032E\ 2 el C\032E\ 3 el 1 rgb :

C\032E\ a\032E\ rgb C\032E\ 1 el C\032E\ 2 el C\032E\ 3 el a\032E\ rgb :

T\032E\ n\032E\ a\032E\ pCMap C\0324\ 0 : r\0324\ 1 n\032E\ range : 12 mat z\0324\ r\0324\ T\0324\ f 02550 - n\032E\ 1 - / r\0324\ 1 - (* + (3 identity z\0324\ 1 el row * T\0324\ + a\032E\ rgb eval 11 line : n\032E\ 0 ≡ 0 T\032E\ IsString T\032E\ Jet ≡ m\0324\ 0.25 n\032E\ * 0 round 112 mat max eval : v\0324\ z\0324\ v\0324\ rows 1 matrix : 12 mat x\0324\ 2551 m\0324\ range m\0324\ / * : y\0324\ 255 m\0324\ / m\0324\ 2 / m\0324\ range (* : o 255 m\0324\ 1 matrix + : 13 mat R\0324\ y\0324\ z\0324\ o z\0324\ x\0324\ o y\0324\ reverse stack : G\0324\ y\0324\ z\0324\ x\0324\ o x\0324\ reverse y\0324\ z\0324\ stack : B y\0324\ o x\0324\ reverse o z\0324\ y\0324\ z\0324\ stack : m\0324\ R\0324\ length n\032E\ - : C\0324\ r\0324\ el R\0324\ r\0324\ m\0324\ + el G\0324\ r\0324\ m\0324\ + el B r\0324\ m\0324\ + el a\032E\ rgb eval : 71 line T\032E\ GS ≡ R n\032E\ a\032E\ pCMap G n\032E\ a\032E\ pCMap + B n\032E\ a\032E\ pCMap + rgb\0324\ RGB T\032E\ findstr : (1 - ≡ C\0324\ r\0324\ el T\032E\ : C\0324\ r\0324\ el rgb\0324\ r\0324\ 0 f : if if if C\0324\ r\0324\ el T\032E\ iFind r\0324\ T\032E\ f : if if 31 line :

T\032E\ x\032E\\0324\ pCMap T\032E\ x\032E\\0324\ 1 pCMap 11 line :

T\032E\ pCMap T\032E\ 15 pCMap 11 line :

M\032E\# G\032E\ pClr M\0324\ M\032E\# : C\0324\ 0 : r\0324\ 1 M\0324\ length range : 13 mat G\032E\ 0 ≡ M\0324\ length 1 matrix G\032E\ IsString C\0324\ r\0324\ el G\032E\ : a\0324\ b\0324\ n\0324\ 13 mat M\0324\ min M\0324\ max G\032E\ length 13 mat eval : 11 mat C\0324\ r\0324\ el G\032E\ M\0324\ r\0324\ el a\0324\ - (n\0324\ 1 - b\0324\ a\0324\ - / (* 0 round 1 + el eval : 21 line if if 21 line :

pCMJet Jet pCMap :

pCMRed R pCMap :

pCMGreen G pCMap :

pCMBlue B pCMap :

<content> pView </content>

EulerRot(v) returns the Euler rotation matrix. At least in some convention. v\032E\ pEulerRot 1000 v\032E\ 1 el cos v\032E\ 1 el sin - 0 v\032E\ 1 el sin v\032E\ 1 el cos 33 mat v\032E\ 2 el cos 0 v\032E\ 2 el sin 010 v\032E\ 2 el sin - 0 v\032E\ 2 el cos 33 mat * v\032E\ 3 el cos v\032E\ 3 el sin - 0 v\032E\ 3 el sin v\032E\ 3 el cos 000133 mat * eval 11 line :

az\032E\ el\032E\ pView2 a\032E\ az\032E\ : l\032E\ el\032E\ : 12 mat l\032E\ l\032E\ 180 ° * + 360 ° * mod 360 ° * + 360 ° * mod 180 ° * - : l\032E\ 90 ° * > l\032E\ 180 ° * l\032E\ - : a\032E\ a\032E\ 180 ° * + : 12 mat l\032E\ 90 ° * - < l\032E\ 180 ° * - l\032E\ - : a\032E\ a\032E\ 180 ° * + : 12 mat 0 if if a\032E\ a\032E\ 360 ° * mod 360 ° * + 360 ° * mod : a\032E\ cos l\032E\ sin a\032E\ sin * - a\032E\ sin l\032E\ sin a\032E\ cos * 0 l\032E\ cos 32 mat eval 51 line :

<content> pMesh </content>

P cl ls w cs pLine cs 0 ≡ polyline P 11 mat cl ls w stack polygon P 11 mat cl ls w cs stack if 11 line :

M\032E\# N\032E\# g\032E\ GM\032E\# GS\032E\# wm\032E\ LS\032E\ d\032E\ pShow M\0324\ M\032E\# : N\0324\\032E\ N\032E\# : nu\0324\ N\0324\\032E\ 1 el abs : nv\0324\ N\0324\\032E\ 2 el abs : w\0324\ wm\032E\ GM\032E\# 0 ≠ (* : 15 mat P\0324\ M\0324\ g\032E\ * eval : c\0324\ 1 P\0324\ cols range : M\0324\d M\0324\ d\032E\ col : 13 mat gM\0324\ gS\0324\ 12 mat M\0324\d GM\032E\# pClr gS\0324\ M\0324\d GS\032E\# pClr : 12 mat eval : pShowAxis nL\0324\ 3 : L\0324\ M\0324\ g\032E\ pAxis eval : 12 mat nL\0324\ 0 : L\0324\ 0 : 12 mat if r.np\0324\ 1 pMesh.np range : r.nv\0324\ 1 nv\0324\ 1 - (range : r\0324\ 0 : S\0324\ 0 : ns\0324\ 0 : 15 mat u\0324\ 1 nu\0324\ range v\0324\ 1 nv\0324\ range r\0324\ r\0324\ 1 + : S\0324\ ns\0324\ ns\0324\ 1 + : el P\0324\ r\0324\ r.np\0324\ + c\0324\ el gM\0324\ r\0324\ el LS\032E\ w\0324\ gS\0324\ r\0324\ el pLine eval : r\0324\ r\0324\ pMesh.np + 1 - : 31 line for r\0324\ r\0324\ 1 + u\0324\ nu\0324\ ≡ (- : L\0324\ nL\0324\ nL\0324\ 1 + : el P\0324\ r\0324\ r.nv\0324\ + c\0324\ el gM\0324\ r\0324\ el LS\032E\ w\0324\ 0 pLine eval : r\0324\ r\0324\ nv\0324\ + 2 - : 41 line for nL\0324\ 0 ≡ 01 matrix L\0324\ if N\0324\\032E\ 1 el 0 < S\0324\ reverse S\0324\ if stack eval 11 mat 71 line :

M\032E\ N\032E\ g\032E\ GM\032E\ GS\032E\ pShow M\032E\ N\032E\ g\032E\ GM\032E\ GS\032E\ 1 solid 3 pShow 11 line :

T\032E\ M\032E\# G\032E\# pShow T\032E\ vf ≡ vf\0324\ M\032E\# : G\0324\ G\032E\# : k\0324\ 1 vf\0324\ length range : clr\0324\ 0 : 14 mat dG\0324\ clr\0324\ k\0324\ el vf\0324\ k\0324\ el 2 col vf\0324\ k\0324\ el 1 col - norme : eval G\0324\ pClr eval : C\0324\ k\0324\ el polyline vf\0324\ k\0324\ el G\0324\ dG\0324\ k\0324\ el el solid 115 mat transpose : 11 mat 31 line T\032E\ ip ≡ ip\0324\ M\032E\# sys2mat.1 : G\0324\ G\032E\# : dG\0324\ 1 ip\0324\ length range G\0324\ pClr eval : IP\0324\ 0 : n\0324\ 0 : 15 mat k\0324\ 1 ip\0324\ length range ip\0324\ k\0324\ el 0 ≡ continue j\0324\ 1 ip\0324\ k\0324\ el length range IP\0324\ n\0324\ n\0324\ 1 + : el polyline ip\0324\ k\0324\ el j\0324\ el G\0324\ dG\0324\ k\0324\ el el solid 115 mat transpose : for if for IP\0324\ 11 mat 31 line P\0324\ T\032E\ 1 el : R\0324\ M\032E\# : T\0324\ G\032E\# : 13 mat k\0324\ 1 P\0324\ length range p\0324\ P\0324\ k\0324\ el : p\0324\ 2 el p\0324\ 2 el R\0324\ * : P\0324\ k\0324\ el p\0324\ 2 el p\0324\ 2 el 1 col T\0324\ 1 el - p\0324\ 2 el 2 col T\0324\ 2 el - augment : eval : 31 line for P\0324\ 11 mat eval 31 line if if 11 line :

PMesh returns a maitrx suitable for 3D plots. F\032E\# pMesh F\032E\ F\032E\# : nu\032E\ F\032E\ rows 1 - : nv\032E\ F\032E\ cols 1 - : 13 mat np\032E\ pMesh.np : r.np\032E\ 1 np\032E\ range : r.nv\032E\ 1 nv\032E\ 1 - (range : 13 mat M\032E\ 03 matrix : c\032E\ 13 range : 12 mat u\032E\ 1 nu\032E\ range v 1 nv\032E\ range r\032E\ M\032E\ rows : k\032E\ r\032E\ r.np\032E\ + : 12 mat M\032E\ k\032E\ c\032E\ el F\032E\ 0011015 mat k\032E\ r\032E\ - el u\032E\ + 0110015 mat k\032E\ r\032E\ - el v + el c\032E\ el : 21 line for r\032E\ M\032E\ rows : j\032E\ r\032E\ r.nv\032E\ + : 12 mat M\032E\ j\032E\ c\032E\ el M\032E\ u\032E\ 1 - (np\032E\ nv\032E\ * nv\032E\ + 1 - (* np\032E\ nv\032E\ 1 - (* + j\032E\ r\032E\ - 1 - (np\032E\ * - c\032E\ el : 31 line for M\032E\ 51 line :

PMesth returns a maitrx suitable for 3D plots. F# pMesh F F# : c 13 range : M 0 : RC 0 : nu F rows 1 - : nv F cols 1 - : 16 mat np pMesh.np : r.np 1 np range : r.nv 1 nv 1 - (range : r 0 : 14 mat u 1 nu range v 1 nv range k r r.np + : RC k el 0011015 mat k r - el u + i 0110015 mat k r - el v + (* + : r r np + : 31 line for k r r.nv + : RC k el RC u 1 - (np nv * nv + 1 - (* np nv 1 - (* + k r - 1 - (np * - el : r r nv + 1 - : 41 line for k 1 RC length range : M k c el F RC k el Re RC k el Im el c el : eval 51 line :

PMesh returns a maitrx suitable for 3D plots. f\032E\ X\032E\# Y\032E\# pMesh X\032E\ X\032E\# : Y\032E\ Y\032E\# : nr\032E\ X\032E\ length : nc\032E\ Y\032E\ length : r\032E\ 1 nr\032E\ range : c\032E\ 1 nc\032E\ range : F\032E\ 0 : 17 mat X\032E\ rows X\032E\ cols 12 mat 2212 mat ≡ f\032E\ X\032E\ 1 el X\032E\ 2 el Y\032E\ 1 el pRange X\032E\ 3 el X\032E\ 4 el Y\032E\ 2 el pRange pMesh f\032E\ X\032E\ 1 el Y\032E\ 1 el feval IsMatrix F\032E\ r\032E\ c\032E\ el f\032E\ X\032E\ r\032E\ el Y\032E\ c\032E\ el feval : F\032E\ r\032E\ c\032E\ el X\032E\ r\032E\ el Y\032E\ c\032E\ el f\032E\ X\032E\ r\032E\ el Y\032E\ c\032E\ el feval 31 mat : if pMesh if 21 line :

<content> Initials </content>

pMesh.np 5 :

pCycleC 5 :

TOL 1012 -^:

∞ 105^:

stop 10 / 11 line :

Matlab default view γ.M 37.5 ° * - 30 ° * pView2 :

Friedel default view γ.F 145 ° * 48 ° * pView2 :

2D View γ.O 0 ° * 90 ° * pView2 :

Flag pShowAxis 1 :

Fridel matrix 0.8232 - 0.4194 - 0.3827 0.5677 0.6187 - 0.543 0.009 0.6643 0.7474 33 mat

Mesh plot

<content> implicit </content>

Box B 20 - 7011522 mat :

Grid N 74321 mat * :

Scale sc 0.3 0001000 0.3 33 mat :

V 1500 :

S 1000 :

F π y 2^* x z + (* 3 V * - 3 π * y 2^* / π y * y 2^x 2^+ sqrt y 2^z 2^+ sqrt + (* S - 2 π * y * / - :

z Φ F 11 mat :

x y zo surf 111 mat zo 11 mat Φ al_nleqsolve * 11 line :

γ 30 ° * 60 ° * - pView2 : cm.1 Jet 64 0.7 pCMap : cm.2 01926413 mat 1281 pCMap : x y s.1 x y 100 - surf : x y s.2 x y 100 surf : S.1 s.1 B N pMesh : S.2 s.2 B N pMesh : 71 line

S.1 sc * N γ black cm.1 pShow S.2 sc * N γ cm.2 0 pShow

Alvaro

4 appVersion 1.73.9126.0