x = -9.0 : 0.5 : 9.0; y = -9.0 : 0.5 : 9.0; z = zeros(length(x), length(y)); for i = 1 : length(x) for j = 1 : length(x) z(i, j) = 2 * sin(x(i)) * sin(y(i)) / (x(i) * y(i)); end end figure; surf(x,y,z);