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];
}
}这段代码能不能求出逆矩阵啊?有点懵逼了。。。