function errorEq = freqEqError(omega,k)

errorEq = sin(omega) + k*omega.*cos(omega);

end

