**in Matlab. please use a while loop INSTEAD of calling a function. please use comments on lines of code

5.16 Water is flowing in a trapezoidal channel at a rate of Q = 20 m/s. The critical depth y for such a channel must satisfy the equation 0 = 1 - 2B Ac where g = 9.81 m/s, A= the cross-sectional area (m), and B = the width of the channel at the surface (m). For this case, the width and the cross-sectional area can be related to depth y by B = 3 +y and A. = 37t: Solve for the critical depth using (a) the graphical method, (b) bisec- tion, and (c) false position. For (b) and (c) use initial guesses of x, = 0.5 and x, = 2.5, and iterate until the approximate error falls below 1% or the number of iterations exceeds 10. Discuss your results