Bonjour a toutes et a tous,
j'éssaye d'avoir la résolution d'un systéme d'equation non linaire avec 5 équation en utilisant ODE45. la fonction et la suivante


Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
function dx=systfun(t,x)
 
 
dx(1)=((g*A)/x(4))*[x(3)-Zmax_CE-(p0/(rho*g))-(alpha*x(4)*x(1)^2)];
dx(2)=((g*A)/(L_CFco-x(4)))*[(ps/(rho*g))-x(3)+(1/(2*g*A^2))-(alpha*(L_CFco-x(4)))*x(2)^2];
dx(3)= ((-a_CF)^2/(g*a*x(4)))*[x(1)-x(2)+x(5)*sqrt(x(3)-Zi-(po/(rho*g)))]
dx(4)=0;
dx(5)=0;
 
[ti,xi]= ode45('systfun',[0 4],[20 20  0 0 0]):
plot(ti,xi,(:,1);'r')
mais des que je lance il m'affiche un message d'erreur 
Expression or statement is incomplete or incorrect.
auriez vous une idée du problème merci