double norm1( matrix a)		/* 1-norm is maximum absolute column sum */
{
    return max(+/abs(a'));
}
