Système équation différentielle [empty sym]

Fermé
gxxx Messages postés 5 Date d'inscription samedi 23 août 2014 Statut Membre Dernière intervention 18 juin 2016 - 18 juin 2016 à 15:19
Bonjour,
Je suis débutant sur matlab et je rencontre un problème en essayant de programmer un système d'équa diff, voici ce que j'ai essayé:

clear all
close all


syms x(t) y(t) z(t) g m h p wz wx wy

Dx=diff(x)
Dy=diff(y)
Dz=diff(z)
D2x=diff(x,2)
D2y=diff(y,2)
D2z=diff(z,2)
eq1=diff(x,2)==(-h*sqrt(diff(x).^2+diff(y).^2+diff(z).^2)*diff(x)+p*(wy.*diff(z)-wz*diff(y)))/m
eq2=diff(y,2)==(-h*sqrt(diff(x).^2+diff(y).^2+diff(z).^2)*diff(y)+p*(wz.*diff(x)-wx*diff(z)))/m
eq3=diff(z,2)==(-m*g-h*sqrt(diff(x).^2+diff(y).^2+diff(z).^2)*diff(z)+p*(wx.*diff(y)-wy*diff(x)))/m

%S=dsolve(eq1, eq2, eq3)
g=9.81
m=2.5*10^-3
h=7.3*10^-4
p=1.8^-5
wz=-20*2*pi
wx=0.2*pi
wy=-20*2*pi
c1=x(0)==-0.5
c2=y(0)==-0.5
c3=z(0)==0.3
c4=Dx(0)==8
c5=Dy(0)==2
c6=Dz(0)==3
S=dsolve(eq1, eq2, eq3,c1,c2,c3,c4,c5,c6,t)
xSol(t)=S*x
ySol(t)=S*y
zSol(t)=S*z


fplot(xSol)
hold on
fplot(ySol)
fplot(zSol)

mais ça met mets S= [empty sym]... quelqu'un peut m'aider? Merci d'avance .
Je précise que j'ai déjà cherché sur ce forum et sur d'autre mais rien de m'a aidé.