Keywords: Devils curve a=0.8 b=1.svg Illustration of en Devil's curve self-made with en Matlab 03 15 14 July 2007 UTC Oleg Alexandrov v An illustration of Devil's curve function main linewidth and font size lw 40; fs 250; colors red 0 867 0 06 0 14; blue 0 129 205/256; green 0 200 70/256; black 0 0 0; white 0 99 1 1 1; Set up the plotting window figure 1 ; clf; set gca 'fontsize' fs 'linewidth' lw/4 ; hold on; axis equal; grid on; N 500; number of points don't make it big code will be slow Lx 2; Ly 2; X Y meshgrid linspace -Lx Lx N linspace -Ly Ly N ; X and Y coordinates a 0 8; b 1; Z Y 2 Y 2 - a 2 - X 2 X 2 - b 2 ; graph the curves using 'contour' in figure 2 figure 2 ; c stuff contour X Y Z 0 0 ; extract the curves from c and graph them in figure 1 using 'plot' need to do this kind of convoluted work since plot2svg can't save the result of 'contour' but can save the result of 'plot' m n size c ; while n > 0 l c 2 1 ; x c 1 2 l+1 ; y c 2 2 l+1 ; get x and y of contours figure 1 ; plot x y 'color' blue 'linewidth' lw ; c c l+2 n ; m n size c ; end figure 1 ; axis equal; axis -Lx Lx -Ly Ly ; plot2svg sprintf 'Devils_curve_a 0 2g_b 0 2g svg' a b ; saveas gcf sprintf 'Devils_curve_a 0 2g_b 0 2g eps' a b 'psc2' ;</nowiki> fl Rational curves Quartic curves Files by User Oleg Alexandrov from en wikipedia Images with Matlab source code |