Upload New File

Ex_3_3.py

+#!/usr/bin/env python3
+# -*- coding: utf-8 -*-
+"""
+Created on Wed May 10 08:07:11 2023
+
+@author: delvigne
+"""
+
+'''
+Exercice 3.3 Simple gene expression model dX/dt = beta - alpha*X (constitutive gene)
+'''
+
+import numpy as np
+import matplotlib.pyplot as plt
+from scipy.integrate import odeint
+
+
+def SimpleGeneExpression(X,t):
+        alpha  = 0.1  #Removal rate in h-1
+        beta = 0.14    #Production rate in h-1
+        dXdt = beta - alpha*X
+        return dXdt
+
+t0   = 0
+tmax = 50
+dt   = 0.1
+t = np.arange(t0,tmax,dt)
+
+X0   = 0.1 #Initial protein concentration (g/L)
+X  = odeint(SimpleGeneExpression,X0,t)
+
+plt.plot(t,X)
+plt.xlabel('Time (h)')
+plt.ylabel('Protein X')
+plt.show()