for o = odchylenie;
    for s = szum;
        for r = repeatTimes;
        blank = zeros;
            for r1 = 0:r;
                %generowanie szumu
                zaszumiony = imnoise(I_gray,'gaussian',s/255,o/255);
                %dodawanie szumów obrazu
                blank = blank + zaszumiony;
            end

        %usrednienie szumu dzielac zsumowane wartości przez ilość zaszumień
        srednia = blank/r1;
        %dodanie szumu do szarego obrazu
        zaszumiony = srednia + I_gray;
        %pokaz obraz wynikowy
        %figure, imshow(zaszumiony);

        I_zaszumiony_gray = 255 * uint8(zaszumiony);
        roznica =  (I_zaszumiony_gray - I_gray);

        param1 = num2str(o);
        param2 = num2str(s);
        param3 = num2str(r);

        figure('Name',['histogram: odchylenie ', param1 ,' szum ', param2 ,' ilosc_powtorzen ' , param3])
        subplot(3,1,1)
        imshow(zaszumiony)
        subplot(3,1,2)
        imshow(roznica)
        subplot(3,1,3)
        imhist(roznica)

        end
    end
end