The goal of this part was only to check the speed of the
diffusion. We use the same model as the 2011 IGEM team of Grenoble, which is based on the Bangalore 2007 IGEM team model. We first have to add the diffusion terms in the equations thus we get:
Because of the temporal and special derivatives, we couldn’t use a classic matlab solver to solve this set of equations. The approximation we used consisted in dividing the space into a grid:
It is the finite difference method. We thus get:
with h=lx/N, where lx is the length of the grid and N is the number of points of discretization along x. By using the same approximation on yi, and assuming that ly=lx, and that we have the same number of points of discretization, we get:
"
