Интерполяция в каноническом виде. Возможна ошибка в системе Гаусса

#include <string.h>
#include <conio.h>
#include <stdio.h>
#include <windows.h>
#include <stdlib.h>
#include <iostream>
#include <io.h>
#include <math.h>
#define razm 3
#define N razm
 
 
using namespace std;
int main(void)
{ setlocale(LC_ALL, "Russian");
 float massx[razm];
 float massy[razm];
 float mas[razm][razm+1];
  int i, j, k, n=0;
 printf("Введите x и y %i точек:\n",razm);
 for(i=0;i<razm;i++)
 {
     scanf("%f%f",&massx[i],&massy[i]);
 }
  for(i=0;i<razm;i++)
 {
     printf("%f %f\n",massx[i],massy[i]);
 }
for(i=0;i<razm;i++)
{
    for(j=0;j<razm;j++)
    mas[i][j]= (pow(massx[i],j));
    mas[i][razm]=massy[i];
}
  for(i=0;i<razm;i++)
 {
     for(j=0;j<razm+1;j++)
     printf("%f ",mas[i][j]);
     printf("\n");
 }
 
  double x[N]; //Корни системы
  int otv[N]; //Отвечает за порядок корней
 
 
  for ( i = 0; i < n + 1; i++ )
    otv[i] = i;
  //Прямой ход метода Гаусса
  for ( k = 0; k < n; k++ )
  {
 
    if ( fabs( mas[k] [k] ) < 0.0001 )
    {
      printf( "Система не имеет единственного решения" );
      return ( 0 );
    }
    for ( j = n; j >= k; j-- )
      mas[k] [j] /= mas[k] [k];
    for ( i = k + 1; i < n; i++ )
      for ( j = n; j >= k; j-- )
        mas[i] [j] -= mas[k] [j] * mas[i] [k];
  }
  //Обратный ход
  for ( i = 0; i < n; i++ )
    x[i] = mas[i] [n];
  for ( i = n - 2; i >= 0; i-- )
    for ( j = i + 1; j < n; j++ )
      x[i] -= x[j] * mas[i] [j];
  //Вывод результата
  printf( "Ответ:\n" );
    for ( j = 0; j < n; j++ )
 
        printf( "%f \n", x[j] );
 
  return ( 0 );
}