Are you looking for a plot with normalized output?
syms x y
x0 = -.5:0.001:.5;
y0 = x0;
w0 = 0.05;
[X0,Y0] = meshgrid(x0,y0);
omegax = 2;
omegay = omegax;
scale = 1/(sinh((omegax*0.5 + omegay*0.5)/w0));
f = @(x0, y0) scale*sinh((omegax*x0 + omegay*y0)/w0);
surf (X0,Y0,abs(f(X0,Y0)) )
shading interp