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