Aide programme matlab

[invite320d6e59]

bonsoir tout le monde
je débute en Matlab, et je dois résoudre l'équation de dispersion suivante (mécanique des fluides) par la methode de newton-raphson :
w^2=g*kn*tanh(kn*h)
B/X=tanh(X) avec B=hw²/g X=kh
-----

[JPL]

À tout hasard lis http://forums.futura-sciences.com/lo...ces-forum.html.

[invite320d6e59]

Merci JPL , je savais pas voici un programme que j'ai fait mais il y a un probleme
x(1)=input(' donner une valeur initiale de x(1) : \n ');
omega = input('omega= ');
h = input('h= ');
g=9.8 ;
A=(h*omega^2)./g ;
e=1e-10;
n=5000;
for i=2:n
f=tanh(x(i-1))-(A/(x(i-1)));
diff=(1/cosh(x(i-1))^2)+(A/(x(i-1))^2);
x(i)=x(i-1)-f/diff;
if abs(x(i)-x(i-1))<=e
xp=x(i);
fprintf('xp=%f\n',x(i));
break;
end
end
j=1:i;
figure(2);
plot(j,x(j),'*r',j,x(j));
xlabel('Nombre d''itérations');
title('Convergence de la solution de la relation de dispersion par la Méthode de Newt.-Raph.');
disp('Les valeurs successives de x(i) sont :');
x'