Code: Select all
all 200
comp xgrid = %seqa(-3,.01,601)
comp ygrid = %seqa(-3,.01,601)
decl rec f(601,601)
ewise f(i,j) = 1./((2*%pi)*%det(||1,.9|.9,1||)^.5)*exp(%scalar(-1*tr(||xgrid(i)|ygrid(j)||)*inv(||1,.9|.9,1||)*(||xgrid(i)|ygrid(j)||)/2))
spgraph(hfi=1,vfi=1)
gcontour(x=xgrid,y=ygrid,f=f,number=10)
grtext(x=-2,y=-2,notransperent) "X"
grtext(x=-2,y=2,notransperent) "X"
grtext(x=2,y=-2,notransperent) "X"
grtext(x=2,y=2,notransperent) "X"
spgraph(done)