Реализовать метод Гаусса для расчета СЛАУ - C (СИ)

Формулировка задачи:

Всем доброго времени суток) Уважаемые знатоки,помогите разобраться ) нужно реализовать метод Гаусса для расчета СЛАУ на Си (пока только прямой ход) что я делаю не так ? Почему-то программа в некоторых местах считает неправильно,хотя в остальных всё верно

//ввод расширенной матрицы с клавиатуры 
/*прямой ход*/
int n=0,m=0;
n=a[0][0];
for(int i=0;i<size+1;i++)  //получаем единицу в ведущем элементе первой строки
 a[0][i]/=n; 
n=1;
for(int j=0;j<size+1;j++)   
 {
   for(i=j+1;i<size;i++)
   {
     float v=a[i][j];             //получаем 0 по столбцам 
     for(k=0;k<size+1;k++)
      a[i][k]=a[j][k]*v-a[i][k];
   }
 for(m=i;m<size+1;m++)   //получаем единицу в ведущем элементе следующей строки 
  a[n][m]/=a[n][n];
 }

и еще вопрос: когда меняю тип матрицы на double, при вводе любого значения с клавиатуры

double *p,*q;
//циклы по стокам и столбцам 
scanf("%f",q);

в элемент записывается +INF или -INF,а не число ,которое я ввел. Как так?

Решение задачи: «Реализовать метод Гаусса для расчета СЛАУ»

textual
                            Листинг программы 
                        
                            scanf("%lf",q);