function w = omega1(k)

%  linearized Klein-Gordon 

     w = sqrt( k.^2 + 4);