решение СЛАУ методом квадратного корня. найти ошибку в коде.

[slevin_45]

решение СЛАУ методом квадратного корня.Выводит неверные корни.Вот блок-схема и формулы рассчёта массивов d и s

#include<iostream>
#include<math.h>
using namespace std;
int main (void)
{
	double a[6][6],b[6],x[6],y[6],delta,delta2=0,delta3=0,del2,del3,delta4=0,del4,s[6][6],d[6],H;
	int i,j,n,k;
	cout<<"vvedite poryadok matrici";
	cin>>n;
	n=n+1;
	a[1][1]=-4.88; a[1][2]=1;  a[1][3]=0;a[1][4]=0; a[1][5]=0;
	a[2][1]=1;     a[2][2]=-2; a[2][3]=1; a[2][4]=0;a[2][5]=0;
	a[3][1]=0;     a[3][2]=1;  a[3][3]=-2;a[3][4]=1; a[3][5]=0;
	a[4][1]=0;      a[4][2]=0;  a[4][3]=1; a[4][4]=-2;a[4][5]=1;
	a[5][1]=0;      a[5][2]=0;  a[5][3]=0; a[5][4]=1; a[5][5]=-4.88;
	b[1]=0;b[2]=-8;b[3]=-8;b[4]=-8;b[5]=0;
	for(k=1;k<n;k++)//нахождение элементов S и d;
	{
		for(i=1;i<=k-1;i++)
		{
			del2=d[i]*pow(fabs(s[i][k]),2);
			delta2=delta2+del2;
		}
		if((a[k][k]-delta2)>0)
			d[k]=1;
			if((a[k][k]-delta2)<0)
				d[k]=-1;
				s[k][k]=pow(fabs(a[k][k]-delta2),1/2);
				for(j=k+1;j<n;j++)
				{
					for(i=1;i<k-1;i++)
					{
						del4=d[i]*s[i][k]*s[i][j];
						delta4=delta4+del4;
					}
					s[k][j]=(a[k][j]-delta4)/(s[k][k]*d[k]);
				}
	}
	for(k=1;k<n;k++)
	{
		delta=a[k][k];
		if(k!=1)
		{
			for(i=1;i<k-1;i++)
			{
				delta=delta-d[i]*pow(s[i][k],2);
			}
		}
		if(delta>0)
		d[k]=1;
		if(delta<0)
			d[k]=-1;
			s[k][k]=pow(fabs(delta),1/2);
				for(j=k+1;j<n;j++)
				{
					delta=a[k][j];
					if(k!=1)
					{
						for(i=1;i<k-1;i++)
						{
							delta=delta-d[i]*s[i][k]*s[i][j];
						}
					}
					s[k][j]=delta/(s[k][k]*d[k]);
				}
	}
	y[1]=b[1]/(s[1][1]*d[1]);
	for(i=2;i<n;i++)
	{
		delta=b[i];
		for(k=1;k<i-1;k++)
		{
			delta=delta-d[k]*y[k]*s[k][i];
		}
		y[i]=delta/(s[i][i]*d[i]);
	}
	x[n-1]=y[n-1]/s[n-1][n-1];
	for(i=n-2;i>=1;i--)
	{
		delta=y[i];
		for(k=i+1;k<n;k++)
		{
			delta=delta-x[k]*s[i][k];
		}
		x[i]=delta/s[i][i];
	}
	cout<<endl;
	for(i=1;i<n;i++)
	{
		cout<<x[i];
		cout<<endl;
	}
	cout<<"nevyazka";
	 for(i=1;i<n;i++)
				  {
                      double z=0;
                        for(int j=1;j<n;j++)
                        {
                        z=z+(a[i][j]*x[j]);
                        }
							H=fabs(b[i]-z);
							cout<<endl;
							 cout<<"Nevyazka"<<H;
	 }
	system("pause");
	return(0);
}

[Smogg]

можно лишь посоветовать на листочке просчитать алгоритм, а потом пошагово смотреть, где начинаются различия между тем, что должно быть и тем, то есть на самом деле.

[Sciv]

Один дурацкий вопрос: а почему вы сделали программу не точно по представленной блок-схеме, а переставили местами расчеты конкретных величин? В итоге и сами потерялись, и нам задачу только усложнили. Я честно пытался полчаса вникнуть в решение, однако большую часть времени потратил только на сопоставление блоков схемы с кусками кода вместо того, чтобы искать непосредственно проблему.

Это не претензия, упаси боже, простой интерес.