“…[x, y, z] = sphere (4) x = x * radius (5) y = y * radius (6) z = z * radius (7) xinit(i) = rand*length (8) yinit(i) = rand*width (9) zinit(i) = rand*height (2*uf))^2)*uf) (27) Dworadius = 6*pi*vis (28) cubaf = cuba(f) (29) Qfgfacaf3 = Qfgfac*cuba(f) (30) Q1Bupufaf3 = Q1Bupuf*cuba(f) (31) delx(f,g) = xsys(t,f)-xsys(t,g) (32) dely(f,g) = ysys(t,f)-ysys(t,g) (33) delz(f,g) = zsys(t,f)-zsys(t,g) (34) crit = a(f)+a(g) (35) Fx(f,g) = Fmx(f,g)+Fpx(f,g) (36) Fy(f,g)=Fmy(f,g)+Fpy(f,g) (37) Fz(f,g)=Fmz(f,g)+Fpz(f,g)+dvis(f)*U*(xsys(t,f)-a(f))/length (38)…”