[afen]

for i=1:n
fout=sprintf('%03d.fits',i);
load(fout);
end