I couldn't figure out how to write a code for this question

The differential equation dC B(C/K)" ac, dt 1+(C/K)" is used in systems biology to calculate the concentration of an auto-repressing transcription factor protein produced by a cell. In this equation, C is the concentration of the protein (which varies with time), K is a constant called the 'activation coefficient' for the reaction; B is a second constant, known as the 'maximal expression rate', and a is a constant that determines the rate of degredation or dilution of the protein, and n is also a constant. Write a function in your m-file that will calculate dC / dt given values of t, C. Then, use the MATLAB ode solver (e.g., ode45) to compute, and plot the variation of C with time (as usual, the MATLAB solver won't give you a function for C, but instead will return two vectors, one containing values of time, and the other containing values of C. You can plot these. Use the following parameters: B = 0.8, a=0.6, K = 0.25, n = 4. Plot (on the same graph) the solutions for C = 0.1 at time t = 0 and C = 0.2 at time t = 0, for the interval 0