double cond(A)	/* condition number of A in 2-norm */
{
    matrix S;
    int n = min(sizeof(A));

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

    return S(1,1)/S(n,n);
}
