

function v=vel2d(x,t)

global omega gamma f

v=zeros(size(x));

%Armonikos talantotis
v(1,:)=x(2,:);
v(2,:)=-omega^2*x(1,:)-2*gamma*x(2,:)+f*sin(omega*t);






