I need to write a program that uses newton's method to output the void fraction of the bed

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_p^2) (1 - epsilon)^2/epsilon^3 + 7/4(rho v_0^2/D_p) 1 -epsilon/epsilon^3 where eta is the fluid viscosity, lo is the superficial velocity, D_p is the particle diameter, rho is the fluid density, and epsilon 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 flows 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 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 your 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(epsilon)|