function  dy_dt = Fun_Feder_3( t, y, w2, delta, F0_m, we )
   dy_dt(1,1) = y(2);
   dy_dt(2,1) = - w2*y(1) - 2*delta * y(2) + F0_m * sin(we*t);
end