Assume the complete model for negative feedback, dX/dt = beta/(1+(X/K)^n)-alpah x. Use the provided quasi-stochastic solver (MATHEMATICA) or an equivalent script MATLAB and run a simulation several times to produce a distribution of steady states (for fixed parameter values). Plot the output X versus time. Compare the negative autoregulation steady state values distribution to that of a simple regulation. Note that you have to adjust the simple regulation parameters so that the mean values between the two distributions overlap. Plot the coefficient of variation of the negative autoregulation distribution versus K, beta, alpha, and n (i.e. 4 different plots)