метод прогонки решения линейных уравнений

[Laplas]

привет всем. я новенький) уважаемые, посмотрите прогу, она коректная??? и вообще она правильные вещи считает??

#include<stdafx.h>
#include<stdio.h>
#include<conio.h>
#include<math.h>
#define n 7

int main(int argc, char* argv[])
{
	double a[n]={0,1,2,3,2,1,1},b[n]={2,3,5,7,5,3,2},c[n]={1,1,2,3,2,1,0},d[n]={4,10,27,46,27,10,4};
	double A[n],B[n],e,x[n],D;
	int i,k=0;
	for (i=0;i<=n-1;i++)
	
		printf("a[%d]=%2.2f  ",i+1,a[i]);

	printf("\n\n");
	for (i=0;i<=n-1;i++)
	
		printf("b[%d]=%2.2f  ",i+1,b[i]);
	
	printf("\n\n");
	for (i=0;i<=n-1;i++)
	
		printf("c[%d]=%2.2f  ",i+1,c[i]);

	printf("\n\n");
	
	for (i=0;i<=n-1;i++)
	
		printf("d[%d]=%2.1f  ",i+1,d[i]);
	
	printf("\n\PROVERKA DIAGONAL'NUX ELEMENTOV: ");
	for(i=0;i<=n-1;i++)
	{
		if(fabs(b[i])<(fabs(a[i])+fabs(c[i])))
			k++;
			
	}
	if(k==0)
			printf(" VUPOLNENO \n");
		else
			printf(" NE VIPOLNENO!!!!\n");
	

	A[0]=-c[0]/b[0];
	B[0]=d[0]/b[0];
	printf("\nA[1]=%1.3f  ",A[0]);
	printf("\nB[1]=%1.3f\n  ",B[0]);
	for(i=1;i<n-1;i++)
	{	
		e=a[i]*A[i-1]+b[i];
		A[i]=-c[i]/e;
		B[i]=(d[i]-a[i]*B[i-1])/e;
		printf("\nA[%d]=%1.3f  ",i+1,A[i]);
		printf("\nB[%d]=%1.3f\n  ",i+1,B[i]);

	}
	x[n-1]=(d[n-1]-a[n-1]*B[n-2])/(b[n-1]+a[n-1]*A[n-2]);
    for(i=n-2;i>=0;i--)
	{
		x[i]=x[i+1]*A[i]+B[i];
    }
	printf("\n\n REWENIE SISTEMU:\n\n");
    for(i=0;i<=n-1;i++)
       printf("x[%d]=%2.2f  ",i+1,x[i]);
	
	getch();
	return 0;
}

[ISergeyN]

вообще она правильные вещи считает

а проверить нельзя?

[Laplas]

ну она впринципе работает, но я же не знаю правильные ли она ответы выдает..

[MAKEDON]

Частный случай реши на листке, а потом проверь его программой. Или при помощи Debug просмотри по шагу что происходит. Там и поймешь: верно/неверно.

[Laplas]

ок..благодарю)))