
>n=40; X=random(n); Y=random(n); >function map dist(x,y) := log(prod((x-X)^2+(y-Y)^2)) >x=linspace(0,1,500); y=x'; >z=dist(x,y); levels=linspace(totalmin(z),totalmax(z),100); >fullwindow(); plot2d(z,levels=levels,a=0,b=1,c=0,d=1,<grid):

>n=40; X=random(n); Y=random(n); >function map dist(x,y) := log(prod((x-X)^2+(y-Y)^2)) >x=linspace(0,1,500); y=x'; >z=dist(x,y); levels=linspace(totalmin(z),totalmax(z),100); >fullwindow(); plot2d(z,levels=levels,a=0,b=1,c=0,d=1,<grid):