Résoudre équation du type y=f(y)

Bonjour,

Je souhaite trouver la longueur d'onde L qui rempli la condition suivante :

L=(9.81*T^2)*tanh(h*(2*pi/L));

Pour cela je crée deux L : L_1 et L_2 que je viens comparer les deux valeurs avec une boucle IF

L_1=(9.81*T^2)*tanh(h*(2*pi/L_2));

Pour ce faire je crée une boucle FOR qui vient tester différente valeur de L_2 (entre 0.1 et 200, mais dans l'idéal entre 0.1 et + l'infini).

Ensuite je teste l'égalité entre L_1 et L_2, si les deux valeurs sont égales alors L=L_2

sinon la boucle FOR recommence avec une autre valeur de L_2 et ainsi de suite.

Ce principe fonctionne très bien avec un tableau Excel mais je n'y arrive pas sur Matlab. J'ai l'impression que ma boucle FOR et IF ne fonctionne pas et je n'arrive pas à comprendre pourquoi.

Avez-vous déjà fait ce genre de boucle ? Si oui pouvez-vous m'aider s'il vous plait ?

Voici mon code :

close all %Ferme les figures
clear     %Efface les variables en mémoire
clc       %Efface le contenu de la fenêtre de commande

%Input 
h = 90; %[m] water depth
H = 1; %[m] wave height
T = 10; %[s] period 

%longueur d'onde. On prend la forme complète

for L_2=0.1:200 %0.1 est la limite basse car tanh(0) = erreur. 200 est une limite mais dans l'idéal il faudrait que la boucle aille jusqu'à l'infinie. 
        L_1=(9.81*T^2)*tanh(h*(2*pi/L_2));
        if L_1==L_2%Problème car L_1 ne sera jamais parfaitement = à L_S. Est-il possible d'ajuster les chiffres significatifs ?
           L=L_2;
        else
           L_2=L_2+0.1;%0.1 représente le pas auquel L_S progresse.
        end
end  

C = ((9.81*T)/(2*pi))*tanh((2*pi)/L);