<content> FillContour </content>

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

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

2D Contour Plots

n a pCmapJet r 1 n range : m 0.25 n * trunc : M 0 : 13 mat v z v rows 1 matrix : o 255 m 1 matrix + : 12 mat x 2551 m range m / * trunc vectorize : y 255 m / m 2 / m range (* trunc vectorize : 12 mat R y z o z x o y reverse stack : G y z x o x reverse y z stack : B y o x reverse o z y z stack : k R length n - : M r el B r k + el 256 G r k + el * + 2562^R r k + el * + 2563^255 a * trunc * + : M r el M r el 232^2 M r el * 232^1 - ≥ (* - eval : 91 line :

a b n rng a b a - n / 0 n range * + 11 line :

2 f B N pGrid r 1 N 1 el range : c 1 N 2 el range : Z 0 : 13 mat X B 1 el B 2 el N 1 el 1 - rng : Y B 3 el B 4 el N 2 el 1 - rng : Z r c el X r el Y c el f : 13 mat 21 line :

Filled Contour G g pFillContour X Y Z 13 mat G : n.1 X length 1 - : n.2 Y length 1 - : S 0 : 14 mat r 1 n.1 range : c 1 n.2 range : n 1 n.1 n.2 * (range : 13 mat a Z min eval : d g length 1 - Z max a - / eval : 12 mat Z Z a - (d * trunc 1 + vectorize eval : S r c el polygon X r el Y c el X r el Y c 1 + el X r 1 + el Y c 1 + el X r 1 + el Y c el X r el Y c el 52 mat 11 mat g Z r c el el solid 1 g Z r c el el stack : S n el 11 mat 61 line :

<content> Contours </content>

¦AuxIPE 0931061551212515610390116 mat :

¦AuxIPΔ 0110 0.5 1100 0.5 25 mat :

¦AuxIPT 0112211331122111042434344323424000000200003000000000030000400000416 mat :

¦AuxIPB 0110011001100110001111000011110000001111111100000101010110101010416 mat :

Contour Plot G NL pContour X Y Z 13 mat G : h 15 range : CL 0 : 13 mat L G 3 el max G 3 el min NL 2 + rng 2 NL 1 + (range el eval : nl 1 NL range M 0 : m 0 : v 0 : P 0 : H 0 : 15 mat r 1 X length 3 - (range c 1 Y length 3 - (range Δx 0220115 mat r + : Δy 2200115 mat c + : 12 mat H h el Z Δx h el Δy h el el L nl el - : n 1 H 1 el 0 < (+ 2 H 2 el 0 < (* + 4 H 3 el 0 < (* + 8 H 4 el 0 < (* + eval : ¦AuxIPE n el 0 ≡ continue H 5 el 0 < (n 11 ≡ (n 6 ≡ (| (& n n 11 ≡ (6 * n 6 ≡ (11 * + eval : 0 if P h el X Δx h el el Y Δy h el el 12 mat : k 14 range ¦AuxIPB k n el 0 ≠ v k el P k el P 1 k 4 mod + el P k el - H 1 k 4 mod + el H k el - / H k el * - eval : continue if for k 1 : k 4 ≤ (¦AuxIPT k n el 0 ≠ (& k k 2 + : M m m 1 + : el v ¦AuxIPT k n el el v ¦AuxIPT k 1 + n el el stack eval : for 41 line if 41 line for for CL nl el M : 31 line for CL mat2sys.1 41 line :

11248015 mat H 0 < vectorize * (1 el + eval

x y f B N pGrid

returns a grid G = [X Y Z] with N1xN2 values of f in the box B

G g pFillContour

plots the filled contour of with the data in G with colors in g

B x1 x2 y1 y2 22 mat ≡

N nx ny 21 mat ≡

n α pCmap

creates a Jet colormap of n colors with α transparency

g 200 0.9 pCmapJet :

G ν pContour

plots ν conotur levels with the data in G.

Examples

Select the number of contours:

ν 9 :

x y f x 2^4 y 2^* - : B 5 - 55 - 522 mat : N 606021 mat : 12 mat G x y f B N pGrid : FC G g pFillContour : CC G ν pContour : 51 line

FC CC 2 1 sys FC CC 2 1 sys

x y f x sin y cos + : B 2 π * - 2 π * 04 π * 22 mat : N 808021 mat : 12 mat G x y f B N pGrid : FC G g pFillContour : CC G ν pContour : 51 line

FC CC 2 1 sys FC CC 2 1 sys

Contour of a polar function z(ρ,φ) = sin(2φ)(1-ρ)

x y f ρ φ 12 mat x i y * + abs y x atan 12 mat : 2 φ * sin 1 ρ - (* 21 line : B 1 - 11 - 122 mat : N 10010021 mat : 12 mat G x y f B N pGrid : FC G g pFillContour : CC G ν pContour : 51 line

FC CC 2 1 sys FC CC 2 1 sys

Alvaro