function error = freqEq1Error(omega)

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

end

