Составители:
Рубрика:
ПРИЛОЖЕНИЕ
Функция обращения матрицы
void omat(float A[10][10],int N)
float s,r;
int i,j,k;
for(i=0;i<N;i++)
{
for(j=N;j<2*N;j++)
{
A[i][j]=0;
}
A[i][i+N]=1;
}
for(k=0;k<N;k++)
{
s=A[k][k];j=k;
for(i=k;i<N;i++)
{
r=A[i][k];
if (abs(r)>abs(s))
{
s=r;j=i;
}
}
if(j!=k)
{
for (i=k-1;i<2*N;i++)
{
r=A[k][i];
A[k][i]=A[j][i];A[j][i]=r;
}
}
for(j=k+1;j<N*2;j++)
{
A[k][j]=A[k][j]/s;
}
for(i=k+1;i<N;i++)
{
r=A[i][k];
for(j=k+1;j<N*2;j++)
{
23