function s=fracPuiss(mi,ma,m)
% production d'une image fractale: on part d'un nombre complexe
% on calcule An+1=i^An. On teste si après un certain temps on
% trouve la cas limite: 0.43828+0.36059i. Si oui on colorie le
% point en fonction du nombre d'itération, sinon, on met le point 
% à noir.

limR=0.43828;
limI=0.36059;
li=0;
s=zeros(m+1,m+1);
for j=mi:2/m:ma
	li=li+1;
	co=0;
	for k=mi:2/m:ma
		co=co+1;
		d=j+k*i;	% point de départ
		n=0;
		while (abs(real(d)-limR)>1e-5) & (abs(imag(d)-limI)>1e-5) & (n < 256)
			n=n+1;
			d=i^d;
		end
		s(li,co)=n;
	end
end
