function [B,i] = diagdom(A,tol)
%DIAGDOM:  modify the matrix A
% [B,i] = diagdom(A,tol) returns a
% diagonally dominant matrix B by
% modifying the diagonal of A.  i is a
% list of modified diagonal entries.
error('diagdom mexFunction not found');