I=imread(saturn.tif);
figure,imshow(I)
[m,n]=size(I);
J(m,n)=0;
for i=1:N
X=imnoise(I,'gaussian');
    figure,imshow(X)
    Y=double(X);
    J=J+Y/N;
end
figure,imshow(mat2gray(J))