function error = freqEqError(omega,k)

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

end

