C++中求逆矩阵的方法。。。。

double L;
int q,w,e;
for(q=0;q<N-1;q++)
{ for(w=0;w<N-1;w++)
aug[q][w]=yakebi[q][w];

    for(w=N-1;w<(2*N-2);w++)
    {
        if(w==q+N-1)
        {
            aug[q][w]=1;}
        else {aug[q][w]=0;}}}

        for(q=0;q<N-1;q++)
        {
            for(e=q+1;e<N-1;e++)
            {
                L=-aug[e][q]/aug[q][q];
                for(w=q;w<2*(N-1);w++)
                    aug[e][w]=aug[e][w]+L*aug[q][w];
            }
        }
    for(q=N-2;q>0;q--)
    {
        for(e=q-1;e>=0;e--)
        {
            L=-aug[e][q]/aug[q][q];
            for(w=2*(N-1)-1;w>=0;w--)
                aug[e][w]=aug[e][w]+L*aug[q][w];
        }
    }
    for(q=N-2;q>=0;q--)
    for(w=2*(N-1)-1;w>=0;w--)
        aug[q][w]=aug[q][w]/aug[q][q];
    for(q=0;q<N-1;q++)
    {
        for(w=N-1;w<2*(N-1);w++)
        {
            inv_yakebi[q][w-N+1]=aug[q][w];
        }
    }这段代码能不能求出逆矩阵啊？有点懵逼了。。。

http://blog.sina.com.cn/s/blog_754d32dd01014kee.html