Need a matlab code for this assignment. thank you

Consider a biochemical reactor where biomass and a desired product are produced by microbial growth from a single, rate limiting substrate. Denote the biomass concentration as X, the substrate concentration as S, and the product concentration as P. The substrate is continuously supplied to the bioreactor at volumetric flow rate F and concentration S,. A liquid stream is continuously removed from the bioreactor at volumetric flow rate F such that the liquid volume V remains constant. The specific growth rate has the form: A (1 - As (SP) Km+S+ where the constants K, and Pm reflect growth inhibition by high concentrations of sub- strate and product, respectively. The bioreactor is described by the following mass balance equations: P PM 1 dx dt ds dt dP dr -DX + (S, P) X = 1 (X,S,P) D(S, -S) - -(SP)X = 12(X,S,P) Yxis -DP + (op(S, P) + B)X = S(X,S,P) where D = F/V. Consider the following parameter values: Yx/s = 0.4 g/g, a = 2.5 g/g, B = 0.15 hr-. .= 0.80 hr-, Pm = 52 g/L, Km 1.7 g/L, K, = 20 g/L, D = 0.20 hr-, and S, = 10 g/L. 1. Develop a M-file to solve the steady-state model equations with the Matlab code fsolve. Find the washout steady state and the non-trivial steady state. 2. Lincarize the nonlinear differential equation model about the non-trivial steady state to obtain the linearized model dx/dt = Ax. Determine if the steady state is stable by using the Matlab function eig to calculate the cigenvalues of the A matrix. 3. Use the M-file in Part a to solve the nonlinear differential equation model using one of the Matlab integration codes. Make the program sufficiently general such that different integration codes can be used by simply changing the function call. 4. Test your Matlab code by integrating the nonlinear differential equation model with the initial condition set equal to the steady state in Part a. Plot the dynamics responses X(), S(0), and P(t). What do you expect these dynamic responses to look like? 5. Integrate the nonlinear differential equation model with the same initial condition (as Part d) for the dilution rate D = 0.15 hr-1. Repeat this test for D = 0.25 br-1. Plot the dynamics responses X(t), S(t), and Plt) for these two dilution rates. What happens to the biomass concentration for these two cases