Разобрать код рисования

#include <cstdlib>
#include <ctime>
#include <iostream>
#include <GL/glut.h>
 
int width;
int height;
 
float rotation = 0.0;
 
float ambient[] = {0.5,0.5,0.5,1.0};
 
double field[8][8] = {{15,20,48,35},{13,24,38,17},{17,35,40,25},{24,30,40,20}};
int steps_x = 30;
int steps_y = 30;
 
double interpolateCubic(double* field, double x) {
    x -= 1.5;
 
    double f_1 = 6.0 * (field[0] - 2.0 * field[1] + field[2]);
    double f_2 = 6.0 * (field[1] - 2.0 * field[2] + field[3]);
 
    double m_1 = (4.0 * f_1 - f_2) / 15.0;
    double m_2 = f_1 - 4.0 * m_1;
 
    double m[] = {0.0, m_1, m_2, 0.0};
 
    int i;
    double o;
 
    if (x <= -0.5) {
        // phi 1
        i = 0;
        o = -1.5;
    }
    else if (x <= 0.5) {
        // phi 2
        i = 1;
        o = -0.5;
    }
    else {
        // phi 3
        i = 2;
        o = 0.5;
    }
 
    double a = (m[i + 1] - m[i]) / 6.0;
    double b = m[i] / 2.0;
    double c = (field[i + 1] - field[i]) - (m[i + 1] + 2.0 * m[i]) / 6.0;
    double d = field[i];
 
    double p = (x - o);
 
    return a * p * p * p + b * p * p + c * p + d;
}
 
/*double interpolateLinear(double* field, double x) {
    x -= 1.5;
 
    int i;
    double o;
 
    if (x <= -0.5) {
        i = 0;
        o = -1.5;
    }
    else if (x <= 0.5) {
        i = 1;
        o = -0.5;
    }
    else {
        i = 2;
        o = 0.5;
    }
 
    return field[i] + (field[i + 1] - field[i]) * (x - o);
}*/
 
double interpolate2D(double x, double y, double (*interpol)(double*, double)) {
    double field_y[] = {interpol(field[0], y),interpol(field[1], y),interpol(field[2], y),interpol(field[3], y)};
    return interpol(field_y, x);
}
 
void reshape(int w, int h) {
    width = w;
    height = h;
}
 
void draw(double x, double y, double z, double (*interpol)(double*, double))
{
    glPushMatrix();
    glTranslated(x,y,z);
    glRotatef(rotation, 0.0,0.5, 0.0);
        // Рисуем точки
    glColor3f(1.0, 0.0, 0.0);
    for (int x = 0; x < 4; x++) {
        for (int y = 0; y < 4; y++) {
            glPushMatrix();
            glTranslated((double) x - 1.5, field[x][y] / 48, (double) y - 1.5);
            glutSolidSphere(0.08, 10, 10);
            glPopMatrix();
        }
    }
    // Рисуем сетку
    
            glColor3f(0.5, 0.5, 0.5);
    glBegin(GL_LINES);
 
    for (int x = 0; x <4; x++)
    { 
        glVertex3d((double) x - 1.5, -1.5, -1.5);
        glVertex3d((double) x - 1.5, -1.5, 1.5);
        
        for (int y = 0; y <4; y++)
        {
            glVertex3d(-1.5, -1.5, (double) y - 1.5);
        glVertex3d(1.5, -1.5, (double) y - 1.5);
            
            glVertex3d((double) x - 1.5, 1.5, (double) y - 1.5);
            glVertex3d((double) x - 1.5, -2, (double) y - 1.5);
        }}
 
    glEnd();
 
    
 
    glColor3f(0.0, 1.0, 0.0);
    for (int x = 0; x < 4; x++) {
        for (int j = 0; j < steps_y; j++) {
            double y = (double) j / (double) (steps_y - 1) * 3.0;
 
            glPushMatrix();
            glTranslated((double) x - 1.5, interpol(field[x], y) / 48, y - 1.5);
            glutSolidSphere(0.04, 5, 5);
            glPopMatrix();
        }
    }
 
    glColor3f(0.0, 1.0, 0.0);
    for (int i = 0; i < steps_x; i++) {
        for (int y = 0; y < 4; y++) {
            double x = (double) i / (double) (steps_x - 1) * 3.0;
            double temp[] = {field[0][y], field[1][y], field[2][y], field[3][y]};
 
            glPushMatrix();
            glTranslated(x - 1.5, interpol(temp, x) / 48, (double) y - 1.5);
            glutSolidSphere(0.04, 5, 5);
            glPopMatrix();
        }
    }
 
    glColor3f(0.0, 0.0, 1.0);
    for (int i = 0; i < steps_x; i++) {
        for (int j = 0; j < steps_y; j++) {
            double x = (double) i / (double) (steps_x - 1) * 3.0;
            double y = (double) j / (double) (steps_y - 1) * 3.0;
 
            glPushMatrix();
            glTranslated(x - 1.5, interpolate2D(x, y, interpol) / 48, y - 1.5);
            glutSolidSphere(0.03, 4, 4);
            glPopMatrix();
        }
    }
    glPopMatrix();
    
}
 
 
 
 
void display() {
    glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT);
 
    glViewport(0, 0, width, height);
 
    glMatrixMode(GL_PROJECTION);
    glLoadIdentity();
    gluPerspective(60, (double) width / (double) height, 1.0, 100.0);
 
    glMatrixMode(GL_MODELVIEW);
    glLoadIdentity();
    gluLookAt(0, 3, 6, 0.0, 0.0, 0.0, 0.0, 1.0, 0.0);
    
    rotation += 0.5;
 
    // Меняем в случайном порядке точки
    for (int x = 0; x <4; x++)
    {
        for (int y = 0; y < 4; y++)
        {
            if (!(rand() % 1000))
            {
                field[x][y] += ((rand() % 3) - 1) * (rand() % 20);
            }
        }
    }
 
    // Рисуем сцены
    //draw(-2.5,0,0,interpolateLinear);
    draw(2.5,0,0,interpolateCubic);
 
    glutSwapBuffers();
}
 
int main(int argc, char** argv) {
    srand(time(NULL));
 
    glutInit(&argc, argv);
    glutInitWindowSize(600, 400);
    glutInitWindowPosition(100, 100);
    glutInitDisplayMode(GLUT_RGB | GLUT_DOUBLE);
 
    int window = glutCreateWindow("3D Интерполяция");
 
    glEnable(GL_DEPTH_TEST);
    glEnable(GL_LIGHTING);
    glEnable(GL_LIGHT0);
    glEnable(GL_COLOR_MATERIAL);
    
    glLightModelfv(GL_LIGHT_MODEL_AMBIENT, ambient);
    glShadeModel(GL_SMOOTH);
    
    glClearColor(0.0, 0.0, 0.0, 0.0);
 
    glutReshapeFunc(reshape);
    glutDisplayFunc(display);
    glutIdleFunc(display);
    
    glutMainLoop();
 
    return 0;
}