Résolution numérique avec MATLAB de l'équation de la thermocline (océan pacifique)

Pièce jointe 602030 Aidez moi à résoudre numériquement ce système avec c++ ou Matlab.

Bonjour,
Résolution par différences finies (schéma explicite) :

Code:

---

1 2 3 4 5 6 7 8 9 10 11

tau = 0.01; dt = 0.01; k = 100; b = 1; t = [-tau:dt:10]; h(1) = 1; h(2) = 1; for i = 2:length(t)-1 h(i+1) = dt*(-tanh(k*h(i-1))+b*cos(2*pi*t(i)))+h(i); end plot(t, h)

---