Помогите пожалуйста найти ошибку!!

#include <iostream.h>
#include <math.h>
#include <conio.h>
int sign(double x)
{
if (x < 0) return -1;
else return 1;
}

void main()
{
double dd = -2;
double q = -3.23;
const int n = 5;

double a[n][n] = {
{q, 1, 0, 0, 0},
{1, -2, 1, 0, 0},
{0, 1, -2, 1, 0},
{0, 0, 1, -2, 1},
{0, 0, 0, 1, q}
};
double b[n] = {0, dd, dd, dd, 0};
double x[n];

int d[n];
double delta, s[n][n];

for (int k = 1; k <= n; k++)
{
delta = a[k - 1][k - 1];
if (k != 1)
for (int i = 1; i <= k - 1; i++)
delta -= d[i - 1]*s[i - 1][k - 1]*s[i - 1][k - 1];

d[k - 1] = sign(delta);
s[k - 1][k - 1] = sqrt(fabs(delta));

for (int j = k + 1; j <= n; j++)
{
delta = a[k - 1][j - 1];
if (k != 1)
for (int i = 1; i <= k - 1; i++)
delta -= d[i - 1]*s[i - 1][k - 1]*s[i - 1][j - 1];

s[k - 1][j - 1] = delta/(s[k - 1][k - 1]*d[k - 1]);
}
}

double y[n];
y[0] = b[0]/(s[0][0]*d[0]);

for (int i = 2; i <= n; i++)
{
delta = b[i - 1];
for (int k = 1; k <= i - 1; k++)
delta -= d[k - 1]*y[k - 1]*s[k - 1][i - 1];
y[i - 1] = delta/(s[i - 1][i - 1]*d[i - 1]);
}

x[n - 1] = y[n - 1]/s[n - 1][n - 1];
for (i = n - 1; i >= 1; i--)
{
delta = y[i - 1];
for (int k = i + 1; k <= n; k++)
delta -= x[k - 1]*s[i - 1][k - 1];
x[i - 1] = delta/s[i - 1][i - 1];
}
for (k = 0; k < n; k++) cout << " " << x[k];
cout << endl;
for(n=1,n<=k,n++) cout<<" "<<delta;
double maxi=-100;
{
if(maxi<=delta)maxi=delta;
}
cout<<"maxi="<<maxi<<endl;

return;
}





ошибка в том, что получается 5 одинаков чисел, т.е 5 одинаковых невязок!!!
мне надо тут вычислить 5 невязок для пяти Х и вывести на экран эти 5 пять невязок и макс невязку!!!
помогите пожалуйста буду очень благодарен!!!!

0