double norminf( matrix a)	/* infinite-norm is maximum absolute row sum */
{
    return max(+/abs(a));
}
