[w1,k1]=size(obrazek);
suma=0;
for m=1:w1
    for n=1:k1
        suma=suma+(double(obrazek(m,n))-double(obrazek2(m,n)))^2;
    end
end


mse = suma/ (w1*k1)