"
Team:Evry/Programming
From 2013.igem.org
Programming methods
This part exposes our technical and algorithmical choices, the implementation and the numerical resolution of our models.
Models implementation
All our models are using the Python language. Python has several major perks:
- Very good floating numbers management
- High level programming
- Object-oriented
- Not an exclusively numerical tool, which allows much more options than Scilab or Matlab
Numerical resolution of ODEs
We chose not to use the resolution functions included in Scipy to have more coding liberties.
To solve our ODEs, we used a simple Euler method. Here is an example:
Let the following be a Cauchy problem:
Considering the following Taylor expansion:
On défini la suite tel que
et
et ainsi on définie tel que
et
Thus, we can define the Euler method:
Which, standardized, becomes .
In the end, all our differential equation systems are implemented like so:
where N is the number of equations in the system.
Display
We used the Matplot library to display 2D and 3D curves.
