function f = dispersion(k,h,w,g)
f = k*tanh(k*h)-w^2/g;
