I need to write a script for this problem in Matlab.

For a packed bed, the Ergun equation provides a relationship between the pressure drop per unit length and the properties of the bed. Delta P/L = 150 (eta v_0/D^2_p) (1 - elementof)^2/elementof^3 + 7/4 (rho v^2_0/D_P) 1 - elementof/elementof^3 where eta is the fluid viscosity, v_0 is the superficial velocity, D_p is the particle diameter, rho is the fluid density, and elementof is the void fraction of the bed. Consider a 1.5 m long packed bed containing 5 cm diameter particles. Fluid of density 2 g/cm^3 and viscosity 1 cP flow through this bed at a superficial velocity of 0.1 m/s. The pressure drop in the bed is 416 Pa. (a) Write a program that uses Newton's method to determine the void fraction of the bed and outputs it to screen along with the number of iterations and the value of |f(x)|. Provide the equations needed for Newton's method along with evidence from mum program that the solution has quadratic convergence. (You should pick a bad initial guess and a very small tolerance to see the quadratic convergence.) (b) Now let's use the program from the first part to see how the solution speed is affected by the initial guess. Modify your program from part (a) to compute the number of iterations required for convergence with a tolerance |f(elementof)|