##### Reference no: EM131305498

**Problem 1: The Schlogl Model**

This model is famous for its bistable steady-state distribution. The reactions are

c1

B1 + 2X ? 3X,

c2

c3

B2 ? X,

c4

where B_{1} and B_{2} denote buffered species whose respective molecular populations N_{1} and N_{2} are assumed to remain essentially constant over the time interval of interest. There is only one time-varying species, X; the state change vectors are ν_{1 }= ν_{3} = 1, ν_{2} = ν_{4} = -1; and the propensity functions are

a_{1}(x) = c_{1}/2 N_{1}x(x - 1),

a_{2}(x) = c_{2}/6 x(x - 1)(x - 2),

a_{3}(x) = c_{3}N_{2},

a_{4}(x) = c_{4}x.

For some values of the parameters this model has two stable states, and that is the case for the parameter values we have chosen here:

c_{1} = 3 × 10^{-7}, c_{2} = 10^{-4}, c_{3} = 10^{-3}, c_{4} = 3.5,

N_{1} = 1 × 10^{5}, N_{2} = 2 × 10^{5}. (3)

Please write a matlab code and run it on octave online. Please plot the trajectories of x when you start from x = 100, 250 and 500 respectively. Then please modify my SSA code to simulate this system with initial value x = 250, plot four trajectories from your simulation to demonstrate that in the stochastic simulation, the system will get quite different results because of randomness. In your simulation, always use the final simulation time T = 6.

**Problem 2: Oscillating Gene Expression Model**

Use your ODE and SSA codes to generate and plot a trajectory for the following oscillating system shown in Figure 1 with the initial conditions: g (gene) = 1; m (mRNA) = 10; p (protein) = 1000; p2 (protein dimer) = 100; gi (inhibited gene) = 0; E (enzyme) = 100; Ep (enzyme protein complex) = 0; Ep2 (enzyme protein dimer complex)= 0. and final time T = 500.

Figure 1: Diagram and Reaction list of module 2.

**Attachment:-** Assignment.rar