matrix N( matrix a)		/* Null space */
{
    matrix v;
    int rank;

    rank = svd( a, :, :, v);		/* don't need u or s */

    if( rank == cols(a))		/* there is no null space! */
	return [];
    else
	return v(:,rank+1:cols(v));
}
