Вычислить Модуль Вектора |x|

Ruslan737

New member
04.04.2012
2
0
#1
Найти вектор X как решение системы уравнений A(n,n)X(n)=B(n). Вычислить мо-
дуль вектора |X|.
Задание: Необходимо написать и отладить программу с обязательным применением динамических матриц.
Есть какие-то наброски,помогите исправить;)
C++:
float* res(float* *a1, float* b1, int n, float tol)
{
float* *a=new float*[n];
float* b=new float[n];
for(int i=0; i<n; i++)
{
b[i]=b1[i];
a[i]=new float[n];
for(int j=0; j<n; j++)
a[i][j]=a1[i][j];
}
for(int i=0; i<n-1; i++)//Прямой проход; a1 -- матрица, tol -- точность, b1 -- свободные члены
{
if(abs(a[i][i])>tol)
{
for(int j=i+1; j<n; j++)
{
float k1=a[j][i]/a[i][i];
for(int k=0; k<n; k++)
a[j][k]-=k1*a[i][k];
b[j]-=k1*b[i];
}
}
else
{
bool flag=false;
for(int j=i+1; j<n; j++)
if(abs(a[j][i])>tol)
{
flag=true;
float k1;//Буффер
for(int k=0; k<n; k++)
{
k1=a[i][k];
a[i][k]=a[j][k];
a[j][k]=k1;
}
k1=b[i];
b[i]=b[j];
b[j]=k1;
break;
}
if(flag==false)
{
return NULL;//Решений нет
}
i--;
}
}
if(abs(a[n-1][n-1])<tol)
return NULL;
for(int i=n-1; i>0; i--)//Обратный проход
{
for(int j=i-1; j>=0; j--)
{
float k1=a[j][i]/a[i][i];
for(int k=0; k<n; k++)
a[j][k]-=k1*a[i][k];
b[j]-=k1*b[i];
}
}
for(int i=0; i<n; i++)//Привели матрицу к единичной
b[i]/=a[i][i];
delete a;
return b;
}