double norm2(A)	/* 2-norm for matrices */
{
    matrix S;

    (void) svd( A, :, S);

    return S(1,1);
}
