function f=Geschwindigkeit(m,Fa,d,v,rho,cw,A,dT)
N=10000;

for i=1:N
    v(i+1)=v(i)+(Fa/m-d/m*v(i)-rho/2/m*cw*A*(v(i))^2)*dT;
end
f=v;
end