Two things you cannot understand. When building implicit functions using this technology, you need to move all the members to the RHS functions so that your implicit function becomes the following:
0 = -1+ x^2 + 4*(y^2) + x*y
Then using the value of the zero loop will make sense:
x<-seq(-1.1,1.1,length=1000) y<-seq(-1,1,length=1000) z<-outer(x,y,function(x,y) 4*y^2+x^2+x*y -1 ) contour(x,y,z,levels=0)
In the first version, I was wrong. @mnels' was right.

42-
source share