équation différentielle

[invitea536a859]

Bonjour tout le monde, je veux résoudre cette équation avec matlab mais j'arrive ps à trouver mon erreur :
E1u'(x)=gu1
d'abord j'ai fait un programme pr calculer le E1:
clear all
clc
%Simplification
muw=6.5e-4;mu0=1.793e-3;J0=1e-4;gamma=0.5;kp=0.2;km=0.1;epsi =0.1;
syms up u x;
k1=kp;
k2=((kp-km)/(2*epsi))*x+(kp+km)/2;
k3=km;
mw=(1/muw)*(1-u)^(3+(2/gamma));
m0=(1/mu0)*(u^2)^(1-((1-u)^(1+(2/gamma))));
lambda=m0/(m0+mw);
Jp=(J0/gamma)*(1-u)^((-1-gamma)/gamma);
g=lambda*mw*m0;
gu1=0.0695;
u1=0.9;
E1=k1*lambda*mw*Jp
puis un autre pour resoudre l'équation:
function du=solanal(x,u);
%solution analytique
clear all
clc
muw=6.5e-4;mu0=1.793e-3;J0=1e-4;gamma=0.5;kp=0.2;
syms up u x;
mw=(1/muw)*(1-u)^(3+(2/gamma));
m0=(1/mu0)*(u^2)^(1-((1-u)^(1+(2/gamma))));
lamda=m0/(m0+mw);
Jp=(J0/gamma)*(1-u)^((-1-gamma)/gamma);
gu1=0.0695;
du=gu1/E1;
end
après
[x,u]=ode45('solanal',0,1,0.7);
-----

[invitea536a859]

salut j'ai fait une modification
function du=solanal(x,u);
%solution analytique
clear all
clc
muw=6.5e-4;mu0=1.793e-3;J0=1e-4;gamma=0.5;kp=0.2;
syms u x;
mw=(1/muw)*(1-u)^(3+(2/gamma));
m0=(1/mu0)*(u^2)^(1-((1-u)^(1+(2/gamma))));
lambda=m0/(m0+mw);
Jp=(J0/gamma)*(1-u)^((-1-gamma)/gamma);
gu1=0.0695;
E1=kp*lambda*mw*Jp;
du=gu1/E1;
end
après
[x,u]=ode45('solanal',0,1,0.7);

[invite06622527]

Salut ,
tu as écrit :
<<Bonjour tout le monde, je veux résoudre cette équation avec matlab mais j'arrive ps à trouver mon erreur :
E1u'(x)=gu1 >>
Quels sont les significations (ou les définitions) de E1, de g, de u1 ?
Qu'est ce qui est connu et qu'est ce qui est inconnu ?
Dans ton programme, il y a des tas de variables et de constantes mw, m0, gamma, J0 et beaucoup d'autres qui tombent du ciel sans que tu dises de quoi il s'agit. Comment veux-tu qu'on trouve une éventuelle erreur dans ces conditions ?