(→Progamming code) |
(→Progamming code) |
||
Line 24: | Line 24: | ||
{{hidden | {{hidden | ||
|Our computational model | |Our computational model | ||
− | /* | + | |/* |
− | + | | * Program to run a model of bacteria in a fermenter | |
− | + | | * with a 2 plasmid contention system. | |
− | + | | */ | |
− | + | | | |
− | #include <stdio.h> | + | |#include <stdio.h> |
− | #include <stdlib.h> | + | |#include <stdlib.h> |
− | #include <string.h> | + | |#include <string.h> |
− | #include <stdint.h> | + | |#include <stdint.h> |
− | #include <sys/time.h> | + | |#include <sys/time.h> |
− | #include <gsl/gsl_rng.h> | + | |#include <gsl/gsl_rng.h> |
− | #include <gsl/gsl_randist.h> | + | |#include <gsl/gsl_randist.h> |
− | #include <gsl/gsl_sf.h> | + | |#include <gsl/gsl_sf.h> |
− | #include <math.h> | + | |#include <math.h> |
− | + | | | |
− | #define MAX_POPULATION 2E6 | + | |#define MAX_POPULATION 2E6 |
− | + | | | |
− | + | | // Model parameters | |
− | + | | // | |
− | + | | // Fermentation variables | |
− | double Sf = 10.; // Substrate concentration (g/l) in feed solution | + | |double Sf = 10.; // Substrate concentration (g/l) in feed solution |
− | double S = 0.19; // Current substrate concentration in fementer (g/l) | + | |double S = 0.19; // Current substrate concentration in fementer (g/l) |
− | double V = 0.01; // Simulation volume with about 10^6 bacteria in ml | + | |double V = 0.01; // Simulation volume with about 10^6 bacteria in ml |
− | double Vtot = 100.; // Volume of fermenter | + | |double Vtot = 100.; // Volume of fermenter |
double D = 100.; // Dilution rate ml/hr | double D = 100.; // Dilution rate ml/hr | ||
double alpha = 3.4e-11;// growth yield g of substrate needed for 10^6 cells | double alpha = 3.4e-11;// growth yield g of substrate needed for 10^6 cells |