Please, I need help writing a program for this problem using MATLAB. Thanks for the help.

NOTE: You must submit the listing of any function that you write and a log of your session in MATLAB with publish feature. 1. Colebrook equation and plot The frictional loses in pipe flows are computed using the Colebrook equation given below. It is also often presented in the form of a graph called Moody Chart. Your task is to write a MATLAB function and scripts to generate a plot of f us Re for various values of Take the range of Re E [ 1000, 10000000 and E.-|0 0.001 0.002 0.005 0.010 0.020 0.050] The lines for every value of should be on the same graph. The Colebrook equation is given below. f=|-410g, le/D) , 1.256]1-2 -4log10 3.7 Re 3.7 Revf fsolve Note that f appears implicitly in this equation and so you must use fsolve to solve for various values of f for a give set of Re,e/D. The graph should appear as shown below with all of the labels and titles as shown in the graphs. In doing this you will have to learn about the following commands on your own: xlim, xlabel, ylabel, title, grid, text