If you look at the bottom of the code, how can I make it so if the user enters "y" or "Y" the program runs again? Right now it just quits either way.

Code:

#include <stdio.h>
#include <math.h>
#define N 100
int main()
{
//bounds for the multidimensional array which can handle up to 100 variables
float coeff[N][N+1]; /*={{4,2,1,11},{2,3,4,20},{3,5,3,22},};*/
int i,j,k; // loop counters
int n=3; // Number of equations. Can be more, but for now set at 3
float pivot;
float sum;
float x[N];
char ch;
do
{
printf("To use this program enter all coefficients of variables\n");
printf("Start with the first row and first column, and move on to the next row after each one is entered\n");
printf("\nEnter the number of variables\n");
scanf("%d",&n);
for(i=0;i<n;i++)
{
for(j=0;j<n+1;j++)
{
printf("\nEnter row %d col %d element\n",i+1,j+1);
scanf("%f",&coeff[i][j]);
}
}
for(k=0;k<n;k++)
{
for(i=k+1;i<n;i++)
{
// This determines the varible value we need to pivot
pivot=coeff[i][k]/coeff[k][k];
for(j=k;j<n+1;j++)
coeff[i][j]=coeff[i][j]-pivot*coeff[k][j];
}
}
x[n-1]=coeff[n-1][n]/coeff[n-1][n-1]; //determines value of variables (x's)
for(i=n-2;i>=0;i--)
{
// store the sum of variables in here in order to substitute back into the matrix
sum=0;
for(j=i;j<n+1;j++)
sum=sum+coeff[i][j]*x[j];
x[i]=(coeff[i][n]-sum)/coeff[i][i];
}
for(i=0;i<n;i++)
printf("\nx%d = %g\n",i+1,x[i]);
printf("\nDo you wish to continue[y/n]\n");
fflush(stdin);
scanf("%c",&ch);
}while(ch=='y' || ch=='Y');
getchar();
}