Fri Feb1508:15:122002 rhog =0.9:8 Density of air at 3000m mu =1. 693e-5: t Viscosity of air at 3000m b Check on conver 9.8: t Gravi 7*4pi/3a; t mass of ice particl fprintf( Maximum number of Newton sub-iterations occurred\n')i vI:, n)=w z0= 0: 6 Initial particle locatio t Set time length of integration, and number of steps B Plo t- linspace 10, Tax, N+l) dt = Tmax/N ewton-Raphson convergence parameters ubplot (211) t Initialize vector for oDE integration title(Ba Euler integration' 212); xlabel('time'l: ylabel('u 8 Start iterative loop for n =2: N+1 t Begin sub-iteration loop ng this forces at least one sub-iteration to occur while((m Mmax)&(curres restol)) u"w(} rhog,u·2*a}/mu; pRe_pu rhog* (2*a)/mu =24/Re+6/1+sqrt(Re))+0.4 pCD pRe=-24"Re^(-2)-6*(1+sqrt{Re)^(-2}*0.5/sqrt(Re) 0.5·hogu^2pi"a^2·cD; pDpu=0,5*rhog*u^2p1·a^2· pCD_ pRepRe_pu+rhog"pia^2“uCD; late right-hand sides us linearization of f with respect to w 1)=-pD_pu/mp: i this is pfu dual with respect to w 5 Calculate dw from pres_dwdw = -res