function f=dy1(t,x,a)
f=-a*x+sin(t);