1.22 e) Is refering to the book problem from Applied Numerical Methods with MATLAB for Engineers and Scientists (4th Edition)

Create an M-file which IS a function M-file, and which solves Problem 1.22 e) except for N steps rather than 8 steps. The sole input should be N, the number of steps, and the output should be the last 5 velocity values of this approximation and a nice graph of these velocities over the entire time interval.

As depicted in Fig. P1.22, a spherical particle settling through a quiescent fluid is subject to three forces: the downward force of gravity (FG), and the upward forces of buoyancy (Fa) and drag (Fo). Both the gravity and buoyancy forces can be computed with Newton's second law with the latter equal to the weight of the displaced fluid. For laminar flow, the drag force can be computed with Stoke's law, where ? the dynamic viscosity of the fluid (N s/m2), d : the particle diameter (m), and u-the particle's settling velocity (m/s). The mass of the particle can be expressed as the product of the particle's volume and density, ps (kg/m), and the mass of the displaced fluid can be computed as the product of the particle's volume and the fluid's density, ? (kg/m3) The volume of a sphere is Td /6. In addition, laminar flow corresponds to the case where the dimensionless Reynolds number, Re, is less than 1, where Re-?dulu (a) Use a force balance for the particle to develop the differential equation for duldt as a function of d, ?.ps, and ?. (b) At steady-state, use this equation to solve for the particle's terminal velocity. (c) Employ the result of (b) to compute the particle's terminal velocity in m/s for a spherical silt particle settling in water d: 10 ?m, ? : 1 g/cm-Ps-2.65 g/cm3, and ? : 0 014 g/cm s) (d) Check whether flow is laminar. (e) Use Euler's method to compute the velocity from t-0 to 2-5 s with ?-2-18 s given the nitial condition: u(0) 0 Figure. P1.22 Fo FB As depicted in Fig. P1.22, a spherical particle settling through a quiescent fluid is subject to three forces: the downward force of gravity (FG), and the upward forces of buoyancy (Fa) and drag (Fo). Both the gravity and buoyancy forces can be computed with Newton's second law with the latter equal to the weight of the displaced fluid. For laminar flow, the drag force can be computed with Stoke's law, where ? the dynamic viscosity of the fluid (N s/m2), d : the particle diameter (m), and u-the particle's settling velocity (m/s). The mass of the particle can be expressed as the product of the particle's volume and density, ps (kg/m), and the mass of the displaced fluid can be computed as the product of the particle's volume and the fluid's density, ? (kg/m3) The volume of a sphere is Td /6. In addition, laminar flow corresponds to the case where the dimensionless Reynolds number, Re, is less than 1, where Re-?dulu (a) Use a force balance for the particle to develop the differential equation for duldt as a function of d, ?.ps, and ?. (b) At steady-state, use this equation to solve for the particle's terminal velocity. (c) Employ the result of (b) to compute the particle's terminal velocity in m/s for a spherical silt particle settling in water d: 10 ?m, ? : 1 g/cm-Ps-2.65 g/cm3, and ? : 0 014 g/cm s) (d) Check whether flow is laminar. (e) Use Euler's method to compute the velocity from t-0 to 2-5 s with ?-2-18 s given the nitial condition: u(0) 0 Figure. P1.22 Fo FB