resolution numerique equation derivé partielle

bonjour tous le monde
quelqu'un a une idée sur la méthode de résolution de ce système :
I0/dt=Ci*(t)Ci(t)
dI1/dt=Ci(t)*I0(t) dCi/dt=1-0.190*Ci(t)*Ci(t)-2.305*Ci(t)*Cv(t)-0.058*Ci(t)*I1(t)-1.87*10-13*d2Ci/d2z. dCv/dt=1-2.305*Ci(t)*Cv(t).
dCb/dt=0.145*Ci(t)*Ci(t).
dr/dt=0.145*Ci(t)
les variable sont Ci(t) ,Cv(t),Cb(t),r(t) ,I0(t),I1(t),et les conditions initiales sont tous à 0 pour Ciz(=0)=0 et Ci(z=h)=0 qlq soit le temps avec h et l’épaisseur de la lamme c à d z varie de 0 à h.
la resolution numerique soit avec C ou matlab.

Bonjour.

Je rappelle que l'on aide volontier à corriger un code, mais que l'on ne fait ni le code ni la réflexion à la place du demandeur.
Une petite lecture de ceci, te sera très profitable:
http://www.gnurou.org/writing/smartquestionsfr

________________________________________________________________________
Historique de mes créations, et quelques articles:
http://0217021.free.fr/portfolio
Merci d'utiliser Réponse acceptée si un post répond à votre question