for n = 1:nMax
    guess=(n-0.5)*pi;
    omega=fzero(@(omega) freqEq(omega,k),guess);
    fprintf('Frequency %2i: guess: %6.3f; exact: %6.3f\n',...
        n,guess,omega)
end
