我正在尝试绘制轮廓图。图像是用于创建矩阵z并为x和y值的每个组合填充f（x，y）值的信息。 Click for Image。这是我正在努力的部分，无论我如何尝试将f（x，y）函数放入代码中，它都根本无法工作。任何帮助将不胜感激。
以下是我到目前为止的代码。
y <- seq(-3.00,3.00,0.05)
z <- matrix(0,121,121)
p <- 3.1415
s <- sqrt(3)
for(i in 1:121) {
for(j in 1:121) {
z[i,j] <- 1/p*s*exp[-2*(x^2+y^2-x*y)/3]
}}
contour(x,y,z,col="red")