Please Use MATLAB to solve. (from Numerical Methods for Engineers 7th edition)

28.30 Under a number of simplifying assumptions, the steady- state height of the water table in a one-dimensional, unconfined groundwater aquifer (Fig. P28.30) can be modeled with the follow- ing second-order ODE, K + N = 0 Ground surface Infiltration Water table Aquifer Groundwater flow - Confining bed FIGURE P28.30 An unconfined or "phreatic' aquifer. where x = distance (m), K = hydraulic conductivity (m/d), h = height of the water table (m), h = the average height of the water table (m), and N = infiltration rate (m/d). Solve for the height of the water table for x = 0 to 1000 m where h(0) = 10 m and h(1000) = 5 m. Use the following parameters for the calculation: K = 1 m/d and N = 0.0001 m/d. Set the average height of the water table as the average of the boundary conditions. Obtain your solution with (a) the shooting method and (b) the finite- difference method (Ar = 100 m). 28.30 Under a number of simplifying assumptions, the steady- state height of the water table in a one-dimensional, unconfined groundwater aquifer (Fig. P28.30) can be modeled with the follow- ing second-order ODE, K + N = 0 Ground surface Infiltration Water table Aquifer Groundwater flow - Confining bed FIGURE P28.30 An unconfined or "phreatic' aquifer. where x = distance (m), K = hydraulic conductivity (m/d), h = height of the water table (m), h = the average height of the water table (m), and N = infiltration rate (m/d). Solve for the height of the water table for x = 0 to 1000 m where h(0) = 10 m and h(1000) = 5 m. Use the following parameters for the calculation: K = 1 m/d and N = 0.0001 m/d. Set the average height of the water table as the average of the boundary conditions. Obtain your solution with (a) the shooting method and (b) the finite- difference method (Ar = 100 m)