% analysis ZF
for i=247:267;  % input images zf247.bmp~zf267.bmp
   fname=sprintf('zf%.0f.bmp',i);
   im=imread(fname);
   A(:,:,:,i-246)=imadjust(im,[],[],2); % Gamma=2 		
   nameout=sprintf('Azf%.0f.bmp',i);
   imwrite(A(:,:,:,i-246),nameout); % output images Azf247.bmp~Azf267.bmp
end;
mov=immovie(A,JET);	%make a movie
movie(mov,3,5);		%play movie for 3 times, 0.2 sec for each image