function w = omega2(k)

%  linearized Korteweg-deVries  
%   alpha = 10 and beta = 1/3

      w = 10*k - k.^3/3;