Complete this problem using Matlab

13.21 An object can be projected upward at a specified velocity. If it is subject to linear drag, its altitude as a function of time can be computed as mg z = zo + . (vo + mi) (1 e-cm') where z = altitude (m) above the earth's surface (defined as z=0), zo = the initial altitude (m), m = mass (kg), c = a linear drag coefficient (kg/s), vo = initial velocity (m/s), and t = time (s). Note that for this formulation, positive velocity is considered to be in the upward direction. Given the following parameter values: g = 9.81m/s2, zo = 100 m, vo = 60 m/s, m = 90 kg, and c = 14 kg/s, the equation can be used to calculate the object's altitude. Determine the time and altitude of the peak elevation (a) graphically, (b) analytically, and (c) with the golden-section search until the approximate error falls below as = 1% with initial guesses of t1 = 0 and tu = 10 s. All problems must be completed with the use of Matlab 13.21 An object can be projected upward at a specified velocity. If it is subject to linear drag, its altitude as a function of time can be computed as mg z = zo + . (vo + mi) (1 e-cm') where z = altitude (m) above the earth's surface (defined as z=0), zo = the initial altitude (m), m = mass (kg), c = a linear drag coefficient (kg/s), vo = initial velocity (m/s), and t = time (s). Note that for this formulation, positive velocity is considered to be in the upward direction. Given the following parameter values: g = 9.81m/s2, zo = 100 m, vo = 60 m/s, m = 90 kg, and c = 14 kg/s, the equation can be used to calculate the object's altitude. Determine the time and altitude of the peak elevation (a) graphically, (b) analytically, and (c) with the golden-section search until the approximate error falls below as = 1% with initial guesses of t1 = 0 and tu = 10 s. All problems must be completed with the use of Matlab