for i=1:142485

    for j=1:142485

       d(i,j)=image1(i)-image2(j)

    end

end

N = numel(image1);
d = images(:,ones(1,N)) - image2(ones(N,1),:);