function [v] = Euler(f,k,m,a)
w0=sqrt(k/m)
Q=sqrt(k*m)/a
u(1)=0;
v(1)=0;
for i=1:1:5000;
u(i+1)=u(i)- 0.001*(w0*w0*f(i)+w0/Q*u(i)+w0*w0*v(i));
v(i+1)= v(i)+0.001*u(i);
end