Помогите найти ошибку в коде С++

[Aleksandr89]

Помогите найти ошибку в коде С++

Надо решить СЛАУ методом Якоби
вот код:

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

int main()
{
int n = 4;
int A[4][4] = {{2, 0, 7, 1},{0, 2, 0, 0},{1, 2, -1, -1},{-1, -2, 1, 2}};
int C[4] = {11, 6, 2, 0};
float X[4] = {1, 1, 1, 1}, X1[4];
float delta = 0.05;
float maxdelta = 1;
while ( maxdelta >= delta )
{
for (int i = 0; i < n; i++)
{
X1 = -C;
for (int j = 0; j < n; j++)
{
if (i != j)
{
X1 += A[j] * X[j];
}
}
X1 /= -A;
}
maxdelta = fabs(X[0] - X1[0]);
for (int i = 0; i < n; i++)
{
if (maxdelta < fabs(X - X1))
maxdelta = fabs(X - X1);
X = X1;
}
}
for (int i = 0; i < n; i++)
printf("X[%i] = %5.2f\n", i+1, X);
system ("pause");
}

никак не могу найти ошыбку.......

 

[Matrix]

напиши задание подробнее...
так то, код правильный, компилируется и работает!

 

[Aleksandr89]

он компилируется, но работает не правильно (выдает результат: -1.#J). Вот задание: решить СЛАУ методом Якоби:

2*x1+7*x3+x4=11
2*x2=6
x1+2*x2-x3-x4=2
-x1-2*x2+x3+2*x4=0

 

[Long Cat]

int A[4][4] = {{2, 0, 7, 1},{0, 2, 0, 0},{1, 2, -1, -1},{-1, -2, 1, 2}};
int C[4] = {11, 6, 2, 0};
float X[4] = {1, 1, 1, 1}, X1[4];
float delta = 0.05;
float maxdelta = 1;

заменить на

double A[4][4] = {{2, 0, 7, 1},{0, 2, 0, 0},{1, 2, -1, -1},{-1, -2, 1, 2}};
double C[4] = {11, 6, 2, 0};
double X[4] = {1, 1, 1, 1}, X1[4];
double delta = 0.05;
double maxdelta = 1;

 

[Aleksandr89]

не помогает, теперь выдает результат: 0.

 

[Long Cat]

printf("X[%i] = %5.2f\n", i+1, X[i]);

заменить на

printf("X[%d] = %5.2lf\n", i+1, X[i]);

Настоятельно рекомендую код отформатировать и выложив, заключив в блок CODE. Тогда его хотя бы можно будет читать и разбирать.

 

[Aleksandr89]

Отформатировал вот результат:

#include <stdio.h>
#include <stdlib.h>
#include <math.h>

int main()
{
  double n = 4;
  double A[4][4] = {{2, 0, 7, 1},{0, 2, 0, 0},{1, 2, -1, -1},{-1, -2, 1, 2}};
  double C[4] = {11, 6, 2, 0};
  double X[4] = {1, 1, 1, 1}, X1[4];
  double delta = 0.05;
  double maxdelta = 1;
  while ( maxdelta >= delta )
  {
         for (int i = 0; i < n; i++)
         {
             X1[i] = -C[i];
             for (int j = 0; j < n; j++)
             {
                     if (i != j)
                     {
                        X1[i] += A[i][j] * X[j];
                     }
             }
             X1[i] /= -A[i][i];
         }
         maxdelta = fabs(X[0] - X1[0]);
         for (int h = 0; h < n; h++)
         {
             if (maxdelta < fabs(X[h] - X1[h]))
             maxdelta = fabs(X[h] - X1[h]);
             X[h] = X1[h];
         }
  }
   for (int i = 0; i < n; i++)
   printf("X[%d] = %lf\n", i+1, X[i]);
  system ("pause");
}

Результат: -1.#IND00

 

[Long Cat]

if (maxdelta < fabs(X[h] - X1[h]))
             maxdelta = fabs(X[h] - X1[h]);
             X[h] = X1[h];

по определению алгоритма должно быть так

if (maxdelta < fabs(X[h] - X1[h]))
{
             maxdelta = fabs(X[h] - X1[h]);
             X[h] = X1[h];
}

или так

if (maxdelta < fabs(X[h] - X1[h]))
{
             maxdelta = fabs(X[h] - X1[h]);
}
X[h] = X1[h];

?
и да, верните обратно int n=4;

 

[Aleksandr89]

if (maxdelta < fabs(X[h] - X1[h]))
{
             maxdelta = fabs(X[h] - X1[h]);
             X[h] = X1[h];
}

X[h] = X1[h] нужно выполнить обязательно посколько дальше программа работать не сможет

if (maxdelta < fabs(X[h] - X1[h]))
{
             maxdelta = fabs(X[h] - X1[h]);
}
X[h] = X1[h];

и

if (maxdelta < fabs(X[h] - X1[h]))
             maxdelta = fabs(X[h] - X1[h]);
             X[h] = X1[h];

одинаковы посколько {} не обязательно если выполняется одно дейсвие

 

[Long Cat]

Метод Якоби — Википедия
- готовая реализация метода Якоби. Обратите внимание на то, что перед осуществлением итераций требуется выполнить преобразования матрицы в вектора правой части.

 

[Aleksandr89]

У меня тот же алгоритм только масивы статические.