|
.
A B.
X, AX = B.
, xi.
O(N3).
C++
//
double* SLE (matrix & a, double b[])
{
matrix B = matrix (b, a.N (), 0);
matrix c = a.concateMatrix (B);
c = gauss (c);
if (detTriang (c) == 0)
return 0;
double *x = new double[a.N ()];
for (int i = a.N () - 1; i >= 0; -- i)
{
double sum = c[i][a.N ()];
for (int j = i + 1; j < a.N (); ++ j)
sum -= c[i][j] * x[j];
x[i] = sum / c[i][i];
}
return x;
}
17:01 04.01.2011
17:03 04.01.2011
По всем вопросам обращаться: rumterg@gmail.com