function [Matrix] = cumsuminverse(A,dimension)
%[outputmatrix] = cumsuminverse(inputmatrix,dimension)
%dimension: -1 == down to up
%           -2 == right to left
[zeile, spalte] = size(A);
if dimension == -1
    for k = 1:1:spalte
        for j = zeile:-1:2
            A(j,k) = A(j,k) - A(j-1,k);
        end
    end
    Matrix = A;
elseif dimension == -2
    for j = 1:1:zeile
        for k = spalte:-1:2
            A(j,k) = A(j,k) - A(j, k-1);
        end
    end
    Matrix = A;
else
    disp('wrong dimension');
end

end