fe m Fri Feb1508:15:072002 1. 693e-57 6 viscosity of air at 3000m mp =917*4*pi/3*a 3: t mass of ice particle 8 Initial conditions u0 -0.001: 8 Use a small initial velocity to avoid problems with Re=o and CD z0- 0: 8 Initial particle locat et time length of integration, and number of steps t Start iterative loop for n 2: N+1 8 Calculate drag at n-l 6/(1+sqrt(Re)) te right-hand sides at n-1 8 Update using Forward Euler :,n-1)+dt·f t Plot results linspace [0, Tmax, N+ subplot (211) abel('time): ylabel[ 'z' subplot (212) xlabel ('time'): ylabel('u')