disp("reading image");
image=imread("comparagram.bmp");

for i=1:256
    for j = 1:256
       newMatrix(i,j)=log(image(i,j)*(i-4*j));
    end

end
disp("computing result");
    fflush(1);fflush(2);
vec1 =1:256;

vec2=4*ones(256,1);
newMatrix\vec2
%imshow(newMat)
%imshow(Z)