clear all;
N = 500; Tfinal = 5; n = 0;
step = Tfinal/N;
dbkp = zeros(ceil(Tfinal/step),1);
Tbkp = zeros(ceil(Tfinal/step),1);
for T = step:step:Tfinal,
    d = 0.5; 
    for k = 1:1000,
        d = tanh(2*d/T);
    end;
    n = n+1
    dbkp(n) = d;
    Tbkp(n) = T;
end;
close all; % hold on;
plot(Tbkp,dbkp,'b.'); grid on;