Discussion :

[lyn sheill]

i=2.....256
x(i)=f^-1 (f(x(i-1))+(0.0049/x(i-1)))
avec: f(x)=exp(-x² /2)
x(1)=3.65

[Invité]

Bonjour,

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
% g = fonction réciproque de f
% donc g = sqrt(-2*log(x))
 
x(1)=3.65;
for i = 2:256
    a = exp(-0.5*x(i-1)^2)+(0.0049/x(i-1));
    x(i) = sqrt(-2*log(a));
end

[lyn sheill]

mercii beacoup