<p> Utilities </p>

a1703741-2523-4731-94f4-58be6b8cddaa 12 4 5 decimal

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

An easy/accurate isocontour example ... Auto calc

<p> Utilities </p>

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 :

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 :

Mesh U, V  within the bounds L, H L H N xd N intervals U 0 : dx H L - N / : i 1 N 1 + range U i el L dx i * dx - (+ : for U 51 line :

The contour extractor , pair wise in X direction M ε cntr CONTR count 1 : i 1 M rows range j 1 M cols range u count el i : v count el j : 21 line M i j el cntr - abs ε ≤ count count 1 + : cases for for u v augment 31 line :

data Strip2 nest striped in 'R' u data 1 col : v data 2 col : k 0 : i 1 data rows range u i el v i el * 0 ≠ R k k 1 + : el u i el v i el 12 mat : 11 line Undefined if 11 line for 21 line Un-nest R v R 1 el : j 2 R rows range v v R j el stack : for v 81 line :

data n Decimate collect every 'n' pairs 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 31 line :

Generate the symmetric square matrix x y f U V Sym r 1 : i 1 U rows range c 1 : j 1 V rows range M r c el U i el V j el f eval : c c 1 + : 21 line for r r 1 + : 31 line for M 31 line :

M Low High Scale High Low - (M M Min - M Max M Min - 106 -^+ / (* Low + :

r 1 :

L r - :

H r :

high 'N' for clean isocontour. N 96 :

U L H N xd :

V U :

x y f 4 x 2^* 7 y 5^* - 1 x cos * - :

M x y f U V Sym eval :

M M 0255 Scale :

iVBORw0KGgoAAAANSUhEUgAAAGEAAABhCAIAAABJZFj0AAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAAJcEhZcwAADsMAAA7DAcdvqGQAAA+RSURBVHhe7Z2Fj13V2oc37m7F3SkMFAYtMsXdXYu7u7u7c5XiVrwUGLwUh+KaIsEthISQEELuM+dZfVmzzxnJDfv7cpP9/gVrP31/sk7P2VP88ccfv/zyy/fff//5559/9NFH77zzzmuvvfbiiy8+++yzTz75ZGdn58iRIx988MH77rvv7rvvvuuuu+64447bbrvtlltuufnmm2+88cYbbrhhWGOuv/76fzfmX9n8s/r5R/f5e2P+1pjrrrvu2muvveaaa66++uqrrrrqyiuvvPzyyy+99NKLL774wgsvPO+8884+++zTTz/9lFNOOeGEE44++ujDDjvsgAMO2GuvvXbZZZdtt9120003XWeddQYPHlz8/vvvP//88zfffPPpp59+8MEHb7755iuvvPLcc889/fTTjz/++MMPPzxixIgAdPvtt996663Quemmm6QjmsCRHzQfD13pQMSBi2hyOldcccVll112ySWXAOiCCy4Q0BlnnHHqqaeeeOKJxx577BFHHHHwwQfvt99+u++++0477bTVVltttNFGa6655sorr1z89ttvP/3001dfffXxxx+/++67r7/+eizRo48++tBDD91///333HNPDoj1kY5ogkvprIzHrXRkEQMUuQQadgc6rA+ALrroIgGdc845Z5555mmnnXbyyScff/zxLtGBBx64995777bbbttvv/0WW2yxwQYbdHR0rLDCCsWvv/76ww8/fPHFFwjt7bfffvXVV59//nmW6LHHHmOJHnjggXvvvTcAuT7ujnREEzjSScedtTQevaKBiAMX0eR01FczIFR2zDHHsESHHHIIS7TnnnuG0NZdd93VVlutvb29wIy+++67zz77TKG9/PLLo0ePdolU2fDhw/GgHBDrE3QCTY4gHbYxnrjSEYcDFEc00nF9zj///HPPPReJAUgbQmXHHXfcUUcd5RLts88+Q4cO3XHHHRHaxhtvvPbaa2NGgwYNKjSjTz755L333hszZgxCe+aZZ3QilgiV3XnnnZg0EhOQ6+PuBJocRzpsdtzqBhClCS6icXegw/qUAJ100kmqTCfaf//9cetdd91VoW244Yaa0dJLL11gRl9++eXYsWMVGm791FNP5UuEynJArI/KCjpyKR3as7YcH6CKAYojmpwO+jrrrLMEhE8DCJVh1UceeeShhx7KEu2777649c4774zQNttss/XWW2+NNdbAjJZaaqnixx9/JPU//PBDhYZbP/HEEy4RWabKMOkcEOtTohNQSmeN4dBVj1AYuIgm6Lg+pBgelANCZYcffngsEW6t0DbZZBNSXzMaOHBgQTPCjN5//30TDaHh1sQZS4RVozJsiBTDg3JAJTrBpXRcxhNXPRJx4BJomFgfACkxAGFDqiyWaI899sCtt9tuO4RG6q+11lqY0bLLLrvYYosV3377LWZE6lMdSTSFRmnEiVgiuiIqI8XwoBxQ0Mm5pPO2OnGlA44YuTCiYXdifXJAZBkqw6pdIiIft6YWbbPNNght/fXXHzJkCGa0zDLLLLLIIsXXX3+tGVEdSTSFRifCiVwiVYZJ54DcnRxN4EiHbUw6b8UDjhi4BBp3RzrEfA5IlZH3LhGRj1vvsMMOCo3Ux4xWXHHFtra2hRZaqMCwaUaY0UsvvTRq1CiFRicizlwiVYZJ54BcnxKank5c9UjEgYtoSnSI+RyQKjvooINcIt06hEbqY0bLL7/8kksuucACCxQYNs0IM3rhhReojgoNtybOXCJVRorlgFwfAOVo0pFbHbrSkYgDFyfouD7EfA5IlXE7c4l0a4VGvSb1MaPllltuiSWWmG+++QoMm2aEGZH6VEeFhlsTZy6RKiPFckCuj4CkkxNJ521MOnLFA5EY0MTuuD4AIsVyQKqMYu0S4dbUIoVG6mtGtEcMe5555ikwbO76mBGpT3VUaLg1nSiWCJWRYjmgEp3g0nxix3NXOq4ME2hifQBEiglIG1JlFGuXCLemFiE06jWprxnRHjHsueeeu8Cw33rrLZoRqd/Z2anQcGs6EXGGE7FEqIwUw6RzQDmd4JLOm524eXyGKia4iEY6rE8OCBtSZfkS4dYKjXodZkR7xLDnnHPOQsOmGWFGjzzyCNUxhEac5UukSQegnE7OJZ133IkrHXHEyIUBTU6HmFdiAlJlsUS6NULjHkvqa0a0Rwx7wQUXnH322QsaNoZNM9KMqI4kmkLLlwiVAQgPKgEKOjmUdOTGpFNXORJx4OJIR0CsTw4IlWHVFGuWiMjHrUNopH5HRwdmRHvEsOeff/5ZZ521oGFr2JoR1ZFEQ2i6tUsUKssBuT5Bp0QkHbn7pONXM0BxAo3rkwPChkJlFGuXCLfeeuutN998c4SmGa200kq0Rwx73nnnHTBgQEGocZXFsGlGmlEIjU5EnMUSoTJMGg8KQCU6zSDSqascicQEmpwOJp0DQmX5EoXQSH2uIJgRV1kMe9FFFyXUZp555oJQ07BpRpoR1TGEphOxRKqsJSDp5FDSebPx6JWOXBjQSCcHhMQEpMpiiXRrhUbqa0Ya9sILLzzXXHPNNNNMBbcQGraGrRlRHUk0hUYniiVSZXqQgEp0mnGkU1c5EokJNPn6CEgbQmVYtUtE5LNECo16rRmtssoqtEcM21CbYYYZCkKNhs1VduTIkTQjzSiEZpyxRKqsJSDplKCkI3cfH6CikQsDGuk0A1Jl+RLh1ltuuSVCM/UxozBsbiGE2vTTT1+88cYbhBpXWQ2bZoQZUR0VGm4dS6TKMOkcUOxOziWdujHp4FWOUBzRMNIRkBILQFzxXSLcWqFRr0l9zYj2iGEvvvjihNpss8027bTTFmPGjDHURowYwV0/zIgLmkIzzmKJ9CAAxfoEnRKUdOrKJsHIBi5O0GF9ckCqzCUKt1ZomtGqq66qYRNq3NQI/mmmmaYg+A01DZtmhBlRHU00hRZLhMoE5AYFoECTjt/qAWLSc/zVIxcGNNLJAWlD+RLp1iE0zUjDbmtr4xZCqM0yyyxTTz11QfCPGjWKUOMqq2HTjDAjqmMIzSUKlbUElKNJp84mHb/KEQ0TdHoCFEuEWys0Ul8z0rANNW5qBP9UU01VcJsl+A01rrIatmak0HBr4ixfIlTWDCink07dmHTwikcuTk5HQNqQKuN25hLh1gqN1NeMVl99ddojhj1w4EBCzeCfYoopCsoRwc8txFDTsDWjXGg6USxRAMrpJCrduaSDVzlCcUDjBJ0AlC+Rbo3QuMeGGWnYgwYNItS4qRH8M8444+STT15Qjgx+Q23YsGEYNs0IM6I65kLLlygHFLvTDCUdvOJJVBoTuyMgJSagfIlCaGFGQ4YM0bANNYJ/jjnmoBxNNtlkBTd+brMEv7cQQo32aDPCjEi0EFq+REosX5+gkw7emHTwikcuTtBpBpQvUQhNM+KaFoZtqBH8lqNJJ520oEBSjiL4DTUNWzMKoeVLFB6Ur08C0+ofttKRiyMdJgAxAsqXKISmGYVht7e3E2rc1Ah+ytF00003ySSTFBRIyhG3WYPfUNOwNaMQWr5EobLYoJxOOntj0tkrngSmMdDJAYXK8iUKoWlGGjYNG0aEGoy48VOOKJATTzxxQYG0HBn8hloYtmYUQoslagaU00kHzyYdv8oRDQOdEiBV5hLp1iE0zUjDJtS8qVmOBgwYQIHsYjR69OhgRPAbarlhY0YhtHyJWgJKSPoBJT3QXzrSYaQTgGKJcqHlZoRhy4jgh5HlyAI50UQTFZTszs5OCiTlyOCPUNOww4wUmk7kEgWgfH0Sg2oolCYhySbWpwSIJcqFlpuRjAg1g99yFIwmnHDCLkaUbAqk5cjgN9Q07DCjXGilJSoBSk/QmHT2ikc0Tg6o5RIhtDAjDdtQM/gtRxRIS/YEE0xQcBGxZFuODP481GCEGZloCq15iZrppLO3mvQoFYx0mBKgWKJghNBywybUDP5gRIGkZE855ZR/MqJkw8hyFMFvqGnYuRnByCUqbVAznXT2KidRGTfQaQlItw6hhWEHI4LfcmSBDEbjjz9+wWXNi4gF0nIUwW97zM0ohFZaohKg9ATdJz1HNSMdJgCVGLFEJTPCsA01gt9yFIwo2VxEuKz9yYiLSDOjPNRglJuRQiuprJlOOn7Fk9g0pidAudAwo2ZGliMKpCU7GI033niFF1oZUbItkJSjPPjDsFsKLV+ixKYaOglDzyMgRkAlRixRbtgyiuAvMeIiwmWNC21ixIWWy5oXkZyRwR+hlptRs9DyJUrPNG7SE1Q5iU1j8iUKRrnQwrAj+KNAehHJGRVMzsiLiCW7P4xaLlEC05j0BFVOAjNuAlBpifrJyIsIl7VujPxgJBjFRSQKZM6opRn1tETpIaqZhKT7CIhpyQih5aEW5cgC2cyIS/+fjPxgxAutjPKSDaMI/l4Y9R9QeqAKpgSoT0YEvwVSRlxEghGXfhh1fXgUjIYPH96SkSX7L2GUnqOyERDTCyMM2+CPchSMuIj8N4ziIiKjPPibDTtnlKhkk56jykmEemCUh1rvjLj0/z8wSg9R5SQ8jQlGAPqfYcSkR6lyEqH/4z2q/agbozrX+si1uh/JqMd+VPfsvnt2fV/r132tvvf3ce+vPz/q+/Oj+nPIvj+HrD/P7vvz7Pr/Rfr+f5H6/9eskb39/1r9/7RGm9eR1v9PW/9/f16RuI60+P/++nsjfX9vpP7+kbbd2/eP6u+xadu9fY+t/j6klmS0tf4+ZP29Wi3JaGv9vdr6+9lakrZt/Je/n11/z19L0raN//L3/Ovfi2hJYdt+0tbt9yL1747CksK2iTYrUvrdUf37NeQWlqRtG23Ef/r9Wv07SOWmJWnbRhvxT0Xq+h1k/XvakJuWRNvGtok24t+KVP8uux+/y65/36/caABakrZttKXf99fvidC5kRsNQEvStrm1pfdE1O8bwbmVG5bExQ1Lokli20Rbet9I/d4aVgnnVm5aUtg20db13pr6/UesUsiNBqAltefvPxpbv0erz/do1e9jY5UoASE3LYkmiW1zI+HWVr/Xrx/v9avfD+kq4dzKjQagJdEksW1ubfV7RvvxntH6fbWuEiUA51ZuNAAtqc331dbvPXaVhvby3uP6/dmuErUb5w65cSnBkmiS2Hb9HvZ+vIe9fp9/rJLOTVFCbjQALInLLU2y/rsQXYqLVcK5KUrIjcJNA8CSaJL13xfpx98Xqf9ODRnHKtEDqN2UAJwbuXG/pQFoSfXfO0o/9nCVKAE4t3KjcHf4947qv5uF4jBvXYlVwrkpSsiNBqAl1X9/rcuYUJyrxA0O5w650QAGDx78H2FuftRlRgRJAAAAAElFTkSuQmCC

rotate 90°

ε 5 :

level 128 :

set 'N' low Experiment 1: ε 4 ... level 64 Experiment 2: ε 4 ... level 128 Experiment 2: ε 5 ... level 128

Q M ε level CONTR :

XY Q Strip2 : Export the project ≡

Shift/center the system

XY rows 97

~ true 9 D zero 0.65613675303 - :

v i 1 XY rows range v i el XY 2 col i el : for v 2 v * v Min - v Max - v Max v Min - / : 21 line :

i 1 XY rows range u i el XY 1 col i el : for u 2 u * u Min - u Max - u Max u Min - / : 21 line

Export u v augment eval :

pts Export . 9 red plot :

c 00 + 100 cyan 15 mat :

Export Export - zero 1.25 zero 1.25 - 2 2 mat 3 1 sys c pts 2 1 sys

x spline u v x cinterp :

0.65613675303 - spline 0.000000000