Считать координаты из файлов txt и нарисовать кривую по точкам OpenGL

#pragma comment (lib, "opengl32.lib")
#pragma comment (lib, "glu32.lib")
#pragma comment(lib,"glut32.lib")
#pragma comment (lib, "glaux.lib")
#pragma comment(lib, "glew32.lib")
#pragma comment (lib, "legacy_stdio_definitions.lib")
 
#include <C:\Users\User\Desktop\include\glew.h>
//Для Windows - приложения
#include <windows.h>
#include <winuser.h>
 
 
 
//#define UNICODE
//#define WIN32_NO_STATUS
//#define auxDIBImageLoad auxDIBImageLoadA
//#define AUX_RGB 
 
 
 
#include <C:\Users\User\Desktop\include\glaux.h>
#include <C:\Users\User\Desktop\include\glut.h>
#include <GL\GL.h>
#include <GL\GLU.h>
 
//________________________________________________________________________________
 
#include <tchar.h> //библиотека для использования макроса _T() -
//подстановка UNICODE или ASCII символов в зависимости от натсроек приложения
#define WIN_CLASS_NAME _T("winclass1")
//__________________________________________________________________________________
#define GLUT_DISABLE_ATEXIT_HACK
 
 
#include <stdlib.h>
#include <stdio.h>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <math.h>  
#include <iostream>
#include <conio.h>
#include <fstream> 
#include <iomanip>
#include <string.h> 
#include <process.h>
#include <stdarg.h>
 
 
using namespace std;
 
int Window;
 
 
#define WIDTH 700 
#define HEIGHT 650
 
static int scrwidth = WIDTH;
static int scrheight = HEIGHT;
 
 
int n, m;
GLfloat M[160];
GLfloat N[160];
GLfloat O[160];
float x, y, z;
 
 
    void INCL()
    {
        ifstream in1("C:/Users/User/Desktop/X.txt");
        ifstream in2("C:/Users/User/Desktop/Y.txt");
        ifstream in3("C:/Users/User/Desktop/Z.txt");
 
        for (int i = 0; !in1.eof(); i++)
        {
            in1 >> M[i];
            //cout << "X[" << i << "] = " << M[i] << endl;
            M[i] = x;
        }
        for (int i = 0; !in2.eof(); i++)
        {
            in2 >> N[i];
            //cout << "Y[" << i << "] = " << N[i] << endl;
            N[i] = y;
        }
        for (int i = 0; !in3.eof(); i++)
        {
            in3 >> O[i];
            //cout << "Z[" << i << "] = " << O[i] << endl;
            O[i] = z;
        }
 
        glutSwapBuffers();
        glBegin(GL_POINTS);
        glVertex3f(x, y, z);
        for (int i = 0; i < 160; ++i)
        {
            glVertex3f(M[i], N[i], O[i]);
        }
        glEnd();
        glEndList();
        glFlush();
 
 
        
    }
 
    void Initialize()
    {
        glClearColor(0.0, 0.0, 0.0, 0.0); // Background color
        glMatrixMode(GL_PROJECTION);
        glLoadIdentity();
        glOrtho(0.0, 1.0, 0.0, 1.0, -1.0, 1.0);
 
    }
 
    int main(int argc, char** argv)
    {
        glutInit(&argc, argv);
        glutInitDisplayMode(GLUT_DEPTH | GLUT_DOUBLE | GLUT_RGBA | GLUT_ALPHA);
        glutInitWindowPosition(750, 380);
        glutInitWindowSize(300, 300);
        Window = glutCreateWindow("INCL");
        Initialize();
        glClearColor(1.0, 1.0, 1.0, 0.0);
        glMatrixMode(GL_PROJECTION);
        glLoadIdentity();
        glFrustum(-100, 100, -100, 100, 100, 2000); //Создает матрицу перспективной проекции
        //glOrtho(-600,600,-600,600,-4000,4000); //Ортогональная проекция
        glMatrixMode(GL_MODELVIEW);
        glTranslatef(0, 0, -1200);
        glRotatef(-60, 1, 0, 0);//сначала угол поворота 30град, а затем ось Х (1,0,0)
        //glRotatef(90, 0, 0, 1);
        glutDisplayFunc(INCL);
        
 
    }