Help with this and especially the MATLAB portion, too. Thanks!

Consider the ODE system y'_1(t) = alpha - y_1 - 4y_1 y_2/1 + y_1^2, y'_2(t) = beta y_1 (1 - y_2/1 + y_1^2) where alpha and beta are parameters, and y_1 and y_2 are scaled mass fractions of two species taking part in a chemical reaction. The initial conditions are y_1(0) = 0, y_2(0) = 2. The solution is to he obtained on the interval 0 lessthanorequalto t lessthanorequalto 30 with step size h = 0.01. Take alpha = 10. Write a MATLAB code to find the solution by using the RK2 algorithm discussed in class. (a) Run the code for beta = 2. Plot labelled graphs of y_1 against t and y_2 against y_1. Describe the behavior of the solution. (b) Repeat part (a) for beta = 4. In what inspects does the solution behavior differ? Do not print out the numerical values of the solution for either part