function  xp =cap7_se1(t,x)
global qpp u
u=sin(t); a=6; b=2;
xp=-a*x+b*u;
qpp=xp;
end