This needs to be programed as a function in matlab

A molecule A decomposes in a stirred tank reactor according to the following mass balance: dC_A/dt = 1/alpha(C_A0 -C_A)-kC_A where: C_AO = concentration of A at the reactor inlet = 20 (dimensionless) for times less than or equal to 10 min and 40 (dimensionless) for any times greater than 10 min C_a= concentration of A at the reactor outlet t = time (min) alpha = the residence time of the reactor = 5 min k = the rate constant for the decomposition of A = 0.12/min The reactor is well mixed but not at steady state. C_A equals zero at time zero. Write a program that prompts the user for a step size (eg, input) and then uses Euler's method to solve the differential equation from t = 0 to t = 20 min. Use your code to solve for C_A, once using a step size of 0.1 min and once using a step size of 1 min. Plot C_A versus time t for these two step sizes values on the same graph. Plot the lines representing each step size in different colors. Add a title, legend, and axis labels to your plot. In your report, comment on how the results vary with the step size