function out=osc1(t,x,y,a,b,c,w)
out=a*cos(w*t)-b*sin(x)-c*y;