Can you design K value and initial condition in Matlab base on this problem?

2.4.1 Open reaction network. Suppose a reaction network is composed of the following reactions k1Ak2B+CBk3Ck42D2Dk5CDk6 with mass-action rate constants as indicated. a) Construct a differential equation model of the network. b) Determine the steady-state concentrations of all species as functions of the mass-action constants. (b) Based on your model, simulate the dynamics of the system. You can design your own k values and initial conditions. You may need to play around with the k values a little to see any dynamics. ( 3 pts)