function repressilator_paramater_dependency()
	% Repressilatorモデル
	% ODEを数値的に解く
	%% パラメータの指定
	s0 = [4, 5, 8, 4, 5, 6]; % 各mRNA、タンパク質の初期値
	param = [100, 1, 5, 2]; % パラメータ a, a0, n
	time = 0.01:0.1:100; %シミュレーションを行う時間
	a0 = -5:1:1; % パラメータを振る
	a0 = 10.^a0; % パラメータを振る
	col = jet; % 色指定に使用
	figure(1);
	clf;
	figure(2);
	clf;
	%% パラメータを振り分けて計算
	for i = 1:numel(a0);
		param(2) = a0(i); %% 振ったパラメータを代入
		
		%% 数値計算
		[t, time_course] = ode45(@(t, s) ODE(t, s, param), time, s0); % ODEを解く
		
		%% figure 1 を描く
		figure(1);
		hold on;
		subplot(numel(a0), 1, i);
		plot(t, time_course(:,4), 'r', t, time_course(:,5), 'g', t, time_course(:,6), 'b'); % 横軸:時間、縦軸:濃度でプロット
		%         axis([0 time(end) -2 2]);
		
		%% figure 2 を描く
		figure(2);
		%         axis([0 2 -1.2 1.2]);
		ph(i)=plot3(time_course(:,4), time_course(:,5), time_course(:,6), 'color', col(floor(53/numel(a0)*i),:)); % 横軸:TetR、縦軸:lacI、高さ軸:λcIでプロット
		hold on;
		leg1(i).name = ['a = ', num2str(a0(i), '%+6.1e')];
		
	end
	
	%% figure 2 の凡例
	figure(2);
	legend([ph],leg1.name,'Location','NorthEastOutside');
	xlabel('x');
	ylabel('y');
	zlabel('z');
end

function dsdt = ODE(t, s, param)
	x = 　　　　　　　　　; % tetR-lite
	y = 　　　　　　　　　; % lacI-lite
	z = 　　　　　　　　　; % λcI-lite
	X = 　　　　　　　　　; % TetR
	Y = 　　　　　　　　　; % LacI
	Z = 　　　　　　　　　; % λcI
	
	alpha = 　　　　　　　　　;
	alpha0 = 　　　　　　　　　;
	beta = 　　　　　　　　　;
	n = 　　　　　　　　　;
	
	dsdt(1, :) = 　　　　　　　　　; % tetR-lite
	dsdt(2, :) = 　　　　　　　　　; % lacI-lite
	dsdt(3, :) = 　　　　　　　　　; % λcI-lite
	dsdt(4, :) = 　　　　　　　　　; % TetR
	dsdt(5, :) = 　　　　　　　　　; % LacI
	dsdt(6, :) = 　　　　　　　　　; % λcI
end

