4314.sagews

var('x y')

(x, y)

plot3d( log(1 + (x*y)^2) - 2*cos(x - 2*y), (x,-2*pi,2*pi), (y,-2*pi,2*pi) )