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

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