plot((1+cos(x-1))^2 + (1+cos(y-1))^2)