function S = StochExpression(zeit,OMEGA) °/„ file: StochExpression.m
% Stochastische Simulation der Gen Expression I
if (nargin<2)  , DMEGA=1; if (nargin<l), zeit=100; end;end;
% Parameter b = 5.0; AO = 100; dm = 0.35; da = 0.2;
ka = b*dm; km = A0*da/b;
% Anfangsbedingungen
M=0;
A=0;
t=0; ix=l;
% Simulation while (t<zeit);
S(ix,:)=[t M/OMEGA A/GMEGA]; ix=ix+l;
% Berechnung der Raten fuer alle Uebergaenge
wl	= 0MEGA*km
w2	= dm*M;
w3	= ka*M;
w4	= da*A;
rate=wl+w2+w3+w4;
% Ziehe Zeitschritt dt */ dt = -log(1-rand)/rate; t = t + dt;
rate=rate*rand;
rate	= rate	- wl	if	(rate	<	0	0),	M	=M+1;	continue;	end
rate	= rate	- w2	if	(rate	<	0	0),	M	=M-1;	continue;	end
rate	= rate	- w3	if	(rate	<	0	0),	A	=A+1;	continue;	end
rate	= rate	- w4	if	(rate	<	0	0),	A	=A-1;	continue;	end
end;
Dabei ist Q = 1 ein Parameter der Systemgrösse (System-size parameter), der die Fluktuationen bestimmt. Im Limit Q —> oo ergeben sich für M/Q. und A/Q, die deterministischen Konzentrationen [M] und [A] (siehe Vorlesung). Für die gewählten Paremeter befinden sich in der Zelle im Mittel 1.5 mRNA Moleküle und 100 Proteine.
38