Bonjour tout le monde,
Alors voilà, je dois résoudre cet exercice sur matlab :
http://spiral.univ-lyon1.fr/ensccf/f...ts_2011-12.pdf
Projet numéro 26, c'est à dire le tout dernier.
Mon soucis est que je n'arrive pas à bien entrer les equa diff dans matlab, je reçois toujours une erreur. Voici mon code:
Ma Fonction crée:
function z = ff(t,y)
z(t,y) = -k*y(1)*y(2)*Ca0 - y(1)/tau
z(t,y) = -k*y(1)*y(2)*Ca0 + y(1) + (Cb0-y(2))/tau
z(t,y) = k*y(1)*y(2)*Ca0 - y(3)/tau
Mon programme:
%% Données de départ
global Ca0 Ca0p Cb0 k tau Cc y0
Ca0 = 0.1; % Concentration initiale en A [mole/L]
Ca0p = 1; % Concentration ??
Cb0 = 0.4;
k = 2;
tau = 100;
Cc=0; % juste pour fixer valeur
y0 = [Ca0 Cb0 Cc]; % on définit le vecteur avec les concentrations initiales en A et B, et la conc finale en C
% graphique
[t,y]=ode45(@ff,tspan,y0)
Si quelqu'un saurai m'aider, car je suis complètement perdu !
Merci
-----