% function B = perturbtril(A,lb,ub)
%
% Input: A square sparse matrix, 
%        scalars lb, ub
%
% Output: lower triangular B containing the lower triangular part of A 
%         with every entry  perturbed independently by a random factor in [lb ub]

error ('perturbtril mexFunction not found') ;
