Код:
library(plot3D)
test<-function(x,y){
exp(-((x-4)^2+(y-4)^2)^2/1000) + exp(-((x +4)^2+(y+4)^2)^2/1000) + 0.1*exp(-((x +4)^2+(y+4)^2)^2)+0.1*exp(-((x -4)^2+(y-4)^2)^2)
}
x <- seq(-10,10,length=100)
y <- x
z <- outer(x,y,test)
persp3D(x,y,z,theta=-30,phi=30,colkey=F)
title("С новым 2017 годом!","Всем добра, борща и сисек!")