Please use matlab to solve the following problem.

Please include function file or script for function file.

For the main script, a while loop and a for loop (within the while loop) should be used.

The terminal velocity of a spherical particle falling through a fluid is given by the equation: vt=3CD4g(p)Dp Where The coal particle density p=1800kg/m3 Density of water =994.6kg/m3 Gravitational constant g=9.80665m/s2 CD= the coefficient of drag Dp= the diameter of the coal particle (m) The viscosity =8.931104kg/ms The coefficient of drag is a function of the particle Reynolds number and is given by: CD=(Re24)(1+0.14Re0.7)Re=Dpvt For a given particle of coal falling through water, find the terminal velocity vt(m/s) using the Matlab fzero function for 4 separate particle diameters. Set up a while loop in which the particle diameter is varied as follows: 0.15mm,0.20mm,0.25mm&0.30mm. (Hint you must convert from mm to m ) Use the fzero function within the loop to find vt by passing the particle diameter and an initial guess (use vt =0.03m/s ) to the function file for each iteration. Define and save the function file as terminal.m to solve the problem. Within the loop save the values of terminal velocity vt, and particle diameter Dp in a new array and plot the terminal velocity as a function of particle diameter for each value of the diameter without using the hold on command. Also print out the terminal velocity using an fprintf statement in the following format: The terminal velocity is . for amm particle