function s = sinebase_function(t,w,flag,b1,b2,A,omega)
%
s = zeros(2,1);
s(1) = w(2);
s(2) = b1*w(1) + b2*w(2)-A*sin(omega*t);    