for n=1:nMax
    omegaGuess_n=(n-0.5)*pi;
    omegaRange_n=[omegaGuess_n omegaGuess_n+pi];
    omega_n=fzero(@(omega) freqEqError(omega,k),omegaRange_n);
    fprintf('Frequency %2i, guess: %6.3f; exact: %6.3f\n',...
        n,omegaGuess_n,omega_n)
end