my program runs but keeps on crashing afterwards:(

[crisel12]

i dont know what is wrong with my program. i fixed the errors. it runs but in the middle it stopped working. it's a program for inverse matrix. please help me


#include<string.h>
#include<math.h>
#include<conio.h>
#include<iostream>
#include<stdlib.h>
#include <stdio.h>

using namespace std;
int main() ;

void reduction(float a[15][15], int SIZE, int J, int E, float aa[15][15], int S, int R, int C)


{
int i,j, m, n;
float factor;
for (i = 0; i< SIZE; ++i);
{
if (i != J)
{
factor = (a[i][E])/(a[E][E]);
for (j = 0; j < SIZE; ++j)
{
a[i][j] = a[i][j] - a[J][j] * factor;
aa[i][j] = aa[i][j] - aa[J][j] * factor;
}
}
}
}
int main()
{
char ans;
do{

system("cls");

static float mat[15][15]={0};
static float A[15][15]={1};
static float X[15]={0};
float u,FACTOR ;
int i, j, R, C , Rows, Column, g, h, o;
printf("\n\tEnter Number Rows:\n");
printf("\t");
scanf("%d",&R);
printf("\n\tEnter Number Columns:\n");
printf("\t");
scanf("%d",C);
printf("\n\tEnter the elements of the matrix:\n");

for (i = 0; i< R; ++i)
{
X[i]=0;
for (j = 0; j <= C; ++j)
{
mat[i][j] = 0;
}
}
for (Column = 0; Column < C ; ++Column)
{
}
printf(" \n");
for (i = 0; i< R; ++i)
{
for (j = 0; j < C; ++j)
{
std::cin>>mat[i][j];
}
}
//clears the value of the result matrix 'a' to 1 in diagonal.
// 1 0 0 0
// 0 1 0 0
// 0 0 1 0
// 0 0 0 1
for (i = 0; i< R; ++i)
{
for (j = 0; j < C; ++j)
{
if (j == i)
{
A[i][j] = 1;
}
else
{
A[i][j] = 0;
}
}
}
for (i = 0; i< C; ++i)
{
if (mat[i][i]!=0)
{
reduction(mat, C, i, i, A, C, i, i);
}
else
{
for (i = (C-1); i>= 0; --i)
{
if (mat[i][i]!=0)
{
reduction(mat, C, i, i, A, C, i, i);
}
}
}
}
printf("\nInverse Matrix");
for (i = 0; i< C; ++i)
{
for (j = 0; j < C; ++j)
{
if ((i==j) && ((mat[i][j]) != 1 ) && ((mat[i][j]) != 0))
{
for (h = 0; h < C; h++)
{
A[i][h] = A[i][h] * (1/(mat[i][j]));
}
}
}
}
printf("\n");
for (i = 0; i< C; i++)
{
printf("\n");
for (j = 0; j < C; j++)
{
printf("%8.2f", A[i][j]);
}
}
printf("\n");
for (i = 0; i< C; ++i)
{
printf("\n");
for (j = 0; j < C; j++)
{
X[i] = X[i] + (A[i][j] * mat[j][C]);
}
printf("x%d = %8.2f ", i+1,X[i]);
}
printf("\n \nTRY AGAIN? Y/N :");
ans=getche();
}while((ans=='Y')||(ans=='y'));
return 0;
}

[Manga]

#include<string.h>
#include<math.h>
#include<conio.h>
#include<iostream>
#include<stdlib.h>
#include <stdio.h>

using namespace std;

void reduction(float a[15][15], int SIZE, int J, int E, float aa[15][15], int S, int R, int C)

{
	int i, j;
	float factor;
	for (i = 0; i< SIZE; ++i);
	{
		if (i != J)
		{
			factor = (a[i][E]) / (a[E][E]);
			for (j = 0; j < SIZE; ++j)
			{
				a[i][j] = a[i][j] - a[J][j] * factor;
				aa[i][j] = aa[i][j] - aa[J][j] * factor;
			}
		}
	}
}

int main()
{
	char ans;
	do {

		system("cls");

		static float mat[15][15] = { 0 };
		static float A[15][15] = { 1 };
		static float X[15] = { 0 };
		float u, FACTOR;
		int i, j, R, C, Rows, Column, g, h, o;
		printf("\n\tEnter Number Rows:\n");
		printf("\t");
		scanf_s("%d", &R);
		printf("\n\tEnter Number Columns:\n");
		printf("\t");
		scanf_s("%d", &C);
		printf("\n\tEnter the elements of the matrix:\n");

		for (i = 0; i< R; ++i)
		{
			X[i] = 0;
			for (j = 0; j <= C; ++j)
			{
				mat[i][j] = 0;
			}
		}

		for (Column = 0; Column < C; ++Column)
		{
		}
		printf(" \n");
		for (i = 0; i< R; ++i)
		{
			for (j = 0; j < C; ++j)
			{
				std::cin >> mat[i][j];
			}
		}
		//clears the value of the result matrix 'a' to 1 in diagonal.
		// 1 0 0 0
		// 0 1 0 0
		// 0 0 1 0
		// 0 0 0 1
		for (i = 0; i< R; ++i)
		{
			for (j = 0; j < C; ++j)
			{
				if (j == i)
				{
					A[i][j] = 1;
				}
				else
				{
					A[i][j] = 0;
				}
			}
		}
		for (i = 0; i< C; ++i)
		{
			if (mat[i][i] != 0)
			{
				reduction(mat, C, i, i, A, C, i, i);
			}
			else
			{
				for (i = (C - 1); i >= 0; --i)
				{
					if (mat[i][i] != 0)
					{
						reduction(mat, C, i, i, A, C, i, i);
					}
				}
			}
		}
		printf("\nInverse Matrix");
		for (i = 0; i< C; ++i)
		{
			for (j = 0; j < C; ++j)
			{
				if ((i == j) && ((mat[i][j]) != 1) && ((mat[i][j]) != 0))
				{
					for (h = 0; h < C; h++)
					{
						A[i][h] = A[i][h] * (1 / (mat[i][j]));
					}
				}
			}
		}
		printf("\n");
		for (i = 0; i< C; i++)
		{
			printf("\n");
			for (j = 0; j < C; j++)
			{
				printf("%8.2f", A[i][j]);
			}
		}
		printf("\n");
		for (i = 0; i< C; ++i)
		{
			printf("\n");
			for (j = 0; j < C; j++)
			{
				X[i] = X[i] + (A[i][j] * mat[j][C]);
			}
			printf("x%d = %8.2f ", i + 1, X[i]);
		}
		printf("\n \nTRY AGAIN? Y/N :");
		ans = _getche();
	} while ((ans == 'Y') || (ans == 'y'));

	return 0;
}