function Yprime = equa_diffexo2(t,Y)
global k

y = Y(1);
yprime = Y(2);
yseconde = - y + cos(k*t);
Yprime = [yprime;yseconde];