hw4 - Run the simulation[t result = [email protected]..

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
% This program calls a function which implements the basic model of the repressilator, % Close all windows and clear all variables clear all; close all; % Set the simulation to run over 3 hours (3*60*60 seconds) tspan = [0 12*60*60]; % Set the initial conditions of the 5 species in molecules per cell % [mRNA_cI, mRNA_lacI, mRNA_tetR, protein_cI, protein_lacI, protein_tetR] IC = [0, 0, 0, 0, 100, 0];
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: % Run the simulation!! [t result] = ode23s(@repressilator, tspan, IC, ); % Plot the concentrations of the 2 protein species figure(2); plot(t/60, result(:,4), 'b-', t/60, result(:,5), 'r-', t/60, result(:,6), 'g-'); title('Oscillation of the 3 repressor proteins'); legend('cI', 'lacI', 'tetR'); xlabel('Time (min)'); ylabel('Proteins per cell');...
View Full Document

This note was uploaded on 11/11/2011 for the course BIO 20.410j taught by Professor Rogerd.kamm during the Spring '03 term at MIT.

Ask a homework question - tutors are online