3D Gasket.c

#include <stdio.h>
#include <GL/glut.h>

float v[4][3]={{0.0,0.0,1.0},{0.0,0.9,-0.3},{-0.8,-0.4,-0.3},{0.8,-0.4,0.3}};
int n;

void draw_triangle(float *a,float *b,float *c)
{
    glBegin(GL_POLYGON);
    glVertex3fv(a);
    glVertex3fv(b);
    glVertex3fv(c);
    glEnd();
}

void midpoint(float *save,float *a,float *b){
    save[0]=(a[0]+b[0])/2;
    save[1]=(a[1]+b[1])/2;
    save[2]=(a[2]+b[2])/2;
}

void divide_tetrahedron(float *a,float *b,float *c,float *d,int m){
    float ab[3],ac[3],ad[3],bc[3],bd[3],cd[3];
    if(m>0){
        midpoint(ab,a,b);
        midpoint(ac,a,c);
        midpoint(ad,a,d);
        midpoint(bc,b,c);
        midpoint(bd,b,d);
        midpoint(cd,c,d);
        divide_tetrahedron(a,ab,ac,ad,m-1);
        divide_tetrahedron(ab,b,bc,bd,m-1);
        divide_tetrahedron(ac,bc,c,cd,m-1);
        divide_tetrahedron(ad,bd,cd,d,m-1);
    }
    else{
        glColor3f(1.0,0.0,0.0);
        draw_triangle(a,b,c);
        glColor3f(0.0,0.0,1.0);
        draw_triangle(a,c,d);
        glColor3f(0.0,0.0,0.0);
        draw_triangle(c,b,d);
        glColor3f(0.0,1.0,0.0);
        draw_triangle(a,b,d);
    }
}
void display(void){
    glClear(GL_COLOR_BUFFER_BIT|GL_DEPTH_BUFFER_BIT);
    glLoadIdentity();
    divide_tetrahedron(v[0],v[1],v[2],v[3],n);
    glFlush();
}
void myReshape(int w,int h){
    glViewport(0,0,w,h);
    glMatrixMode(GL_PROJECTION);
    glLoadIdentity();
    if(w<=h)
        glOrtho(-2.0,2.0,-2.0*(GLfloat)h/(GLfloat)w,2.0*(GLfloat)h/(GLfloat)w,-10.0,10.0);
    else
        glOrtho(-2.0*(GLfloat)w/(GLfloat)h,(2.0)*(GLfloat)w/(GLfloat)h,-2.0,2.0,-10.0,10.0);
    glMatrixMode(GL_MODELVIEW);
    glutPostRedisplay();
}
int main(int argc,char **argv){
    printf("Enter number of division\n");
    scanf("%d",&n);
    glutInit(&argc,argv);
    glutInitDisplayMode(GLUT_SINGLE|GLUT_RGB|GLUT_DEPTH);
    glutInitWindowSize(500,500);
    glutCreateWindow("3D Gasket");
    glutReshapeFunc(myReshape);
    glutDisplayFunc(display);
    glEnable(GL_DEPTH_TEST);
    glClearColor(1.0,1.0,1.0,1.0);
    glutMainLoop();
}