1. Use Matlab to solve the following problem. Please explain your steps using comments.

1. Compute Q using the composite trapezoidal rule" means DIY with a step size of, say, n=100. In addition, compare youre results to that returned by a built-in integral routine (e.g., the integral sign in Mathcad or "integral" in Matlab). Find numerically the radius at which the flow is half of that at r₀=3cm.

20.17 If the velocity distribution of a fluid flowing through a pipe is known (Fig. P20.17), the flow rate Q (ie., the vol- ume of water passing through the pipe per unit time) can be computed by QvdA, where v is the velocity, and A is the pipe's cross-sectional area. (To grasp the meaning of this relationship physically, recall the close connection between summation and integration.) For a circular pipe, A = 2 and dA = 2 r dr. Therefore, 0 A I FIGURE P20.17