function ID=inarnold(IC,a,b,n)
[h w]=size(IC);
imgn=zeros(h,w);
N=h;
for i=1:n
for y=1:h
for x=1:w
xx=mod((a*b+1)*(x-1)-b*(y-1),N)+1;
yy=mod(-a*(x-1)+(y-1),N)+1 ;
imgn(yy,xx)=IC(y,x);
end
end
IC=imgn;
end
ID=imgn;