function [A, n] = ddbucky()
% DDBUCKY - diagonally dominant modification of bucky()
    A = bucky();
    [n, ~] = size(A);
    A = A + 5*speye(n);  % ensure diagonal dominance
end
