[djamel123]

هذه محاولتي ولا أدري اذا كانت صحيحة لأن الرسم الناتج بدى لي غريبا.

clc
clear

x = 0.1; y = 0.22;
niters = input('How many iterations of the Ikeda map');
xall = zeros(niters,1);
yall = zeros(niters,1);

y13 = 1.0/3.0
y23 = 2.0/3.0
tic 
for k = 1:niters
  xp = x; yp=y;
 if y<=0.5 
   y = 2*yp;
   x = y13*xp;
 else 
  x = y13*xp+y23;
  y = 2*yp - 1;
 end
  xall(k) = x;
  yall(k) = y;
end
 

tm = toc()
clf
figure(1)
plot(xall, yall)