function e = calcerror(n)
    e = eye(n);
    for i=1:n
        exact = pascal(n) / eye(n)(i,:);
        est = linsolve(pascal(n), eye(n)(i,:));

        e(i,:) = abs(exact-est') ./ abs(exact);
    end
end

calcerror(10)
calcerror(15)
calcerror(20)