Помогите с методом ГАУССА (С, не ++)

[bazik]

Написал программку, с тестовой матрицей, но она не хочет работать, уже не знаю, что делать. Есть подозрения, что намудрил со скобками где-то. Спасайте!

#include "stdafx.h"
#include "math.h"

main()
{
	const int N = 4;
	int i, j, s, k;
	double c, d = 1;
	double x[N];
	double b[N] = { 1,2,3,4 };
	double a[N][N] = { 5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20 };

	printf("\n  Начальный массив А:\n");
	for (i = 0; i<N; i++)
	{
		printf("\n");
		for (j = 0; j<N; j++)
		{
			printf("%1.f", a[i][j]);
		}
	}
	printf("\n Нчальный массив В:\n");

	for (i = 0; i<N; i++)
	{
		printf("%1.f\n", b[i]);
	}
	for (s = 0; s<N - 1; s++)
	{
		k = s;
		for (i = s + 1; i<N; i++)
		{
			if (fabs(a[k][s])<fabs(a[i][s]))
				k = i
				if (k != a)
				{
					d = -d;
					c = b[k]; b[s] = b[k]; b[s] = c;
					for (j = s; j<N; j++)
					{
						c = a[s][j]; a[s][j] = a[k][j]; a[k][j] = c;
					}
				}
			for (i = s + 1; i<N; i++)
			{
				c = a[i][s] / a[s][s];
				b[i] = b[i] - b[s] * c;
				for (j = s + 1; j<N; j++)
				{
					a[i][j] = a[i][j] - a[s][j] * c;
				}
			}
			printf("\n Полученный массив А:\n");

			for (i = 0; i<N; i++)
			{
				printf("\n");
				for (j = 0; j<N; j++)
				{
					printf("%1.f", a[i][j]);
				}
			}
			printf("\n\n Полученный массив В:\n");
			for (i = 0; i<N; i++)
			{
				printf("%1.f\n", b[i]);
			}
			for (i = N - 1; i>-1; i--)
			{
				c = b[i];
				for (j = i + 1; j<N; j++)
				{
					c = c - a[i][j] * x[j];
				}
				x[i] = c\a[i][i];
			}
			printf("\n Корни массива: \n");
			for (i = 0; i<N; i++)
			{
				d = d*a[i][i];
			}
		}
		printf("\n det A =%1.f \n", d);
		return 0;
	}
}

[lightfrag]

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

#define N 4

int main()
{

	int i, j, s, k;
	double c, d = 1;
	double x[N];
	double b[] = { 1.,2.,3.,4. };
	double a[][N] = { 5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,15.,16.,17.,18.,19.,20. };

	printf("\n  Начальный массив А:\n");
	for (i = 0; i<N; i++)
	{
		printf("\n");
		for (j = 0; j<N; j++)
		{
			printf("%1.lf", a[i][j]);
		}
	}
	printf("\n Нчальный массив В:\n");

	for (i = 0; i<N; i++)
	{
		printf("%1.lf\n", b[i]);
	}
	
	
	for (s = 0; s<N - 1; s++)
	{
		k = s;
		for (i = s + 1; i<N; i++)
		{
			if (fabs(a[k][s])<fabs(a[i][s]))
				k = i;
				if (k != a)
				{
					d = -d;
					c = b[k]; b[s] = b[k]; b[s] = c;
					for (j = s; j<N; j++)
					{
						c = a[s][j]; a[s][j] = a[k][j]; a[k][j] = c;
					}
				}
			for (i = s + 1; i<N; i++)
			{
				c = a[i][s] / a[s][s];
				b[i] = b[i] - b[s] * c;
				for (j = s + 1; j<N; j++)
				{
					a[i][j] = a[i][j] - a[s][j] * c;
				}
			}
			printf("\n Полученный массив А:\n");

			for (i = 0; i<N; i++)
			{
				printf("\n");
				for (j = 0; j<N; j++)
				{
					printf("%1.lf", a[i][j]);
				}
			}
			printf("\n\n Полученный массив В:\n");
			for (i = 0; i<N; i++)
			{
				printf("%1.lf\n", b[i]);
			}
			for (i = N - 1; i>-1; i--)
			{
				c = b[i];
				for (j = i + 1; j<N; j++)
				{
					c = c - a[i][j] * x[j];
				}
				x[i] = c/a[i][i];
			}
			printf("\n Корни массива: \n");
			for (i = 0; i<N; i++)
			{
				d = d*a[i][i];
			}
		}
		printf("\n det A =%1.lf \n", d);
		return 0;
	}
}

На самом деле практически все лишь ошибки синтаксиса были.

[bazik]

Сообщение от lightfrag

На самом деле практически все лишь ошибки синтаксиса были.

Спасибо! Кстати нам запрещают делать с define, беспощадные численные методы