function f = wave_num(h,w,g)
f = fzero(@(k) dispersion(k,h,w,g),w^2/g);
