Программа на С++ вылетает, не могу понять почему!?

#include <iostream>
#include <stdlib.h>
#include <stdio.h>
#include <allegro.h>
#include <math.h>
#include <time.h>
#define W 1024
#define H 768
 
BITMAP *bmp;
PALETTE pal;
 
void ed();
float va();
float vb();
float vc();
float vd();
void vpi();
void draw();
void clean();
void preobr();
void matrp();
 
int green=makecol(50,0,0);
int sx=0,sy=0;
int i=0,j=0,kg,pl[30][10],nx=W/2,ny=H/2;
float p[4][4],cd[100][4],cd1[100][4],plosk[30][4];
 
using namespace std;
 
int main(int argc, char *argv[])
{
  char s[100];
  gets(s);
  allegro_init();
  install_keyboard();
  install_mouse();
  bmp=create_bitmap(W,H);  
  set_gfx_mode(GFX_DIRECTX_WIN,W,H,0,0);
  set_color_depth(24);
  //text_mode(-1);
  
  FILE *in;
  in=fopen(s,"r");
  fscanf(in,"%d",&kg);
  /*do
  {
      for(j=0;j<4;j++)
      fscanf(in,"%f",&cd[i][j]);
      i++;
  }  
  while(cd[i-1][3]!=-1);
  for(i=0;i<kg;i++)
  {
      j=0;
      do
      {
         fscanf(in,"%d",&pl[i][j]);
         j++;
      }
      while(pl[i][j-1]>=0);
  }     */ 
  fclose(in);
  kg=6;
  cd[0][0]=-50;cd[0][1]=-50;cd[0][2]=-50;cd[0][3]=1;
  cd[1][0]=-50;cd[1][1]=50;cd[1][2]=-50;cd[1][3]=1;
  cd[2][0]=50;cd[2][1]=50;cd[2][2]=-50;cd[2][3]=1;
  cd[3][0]=50;cd[3][1]=-50;cd[3][2]=-50;cd[3][3]=1;
  cd[4][0]=-50;cd[4][1]=-50;cd[4][2]=50;cd[4][3]=1;
  cd[5][0]=-50;cd[5][1]=50;cd[5][2]=50;cd[5][3]=1;
  cd[6][0]=50;cd[6][1]=50;cd[6][2]=50;cd[6][3]=1;
  cd[7][0]=50;cd[7][1]=-50;cd[7][2]=50;cd[7][3]=1;
  cd[8][0]=0;cd[8][1]=0;cd[8][2]=0;cd[8][3]=1;
  pl[0][0]=0;pl[0][1]=1;pl[0][2]=2;pl[0][3]=3;pl[0][4]=-1;
  pl[1][0]=4;pl[1][1]=5;pl[1][2]=6;pl[1][3]=7;pl[1][4]=-1;
  pl[2][0]=0;pl[2][1]=1;pl[2][2]=5;pl[2][3]=4;pl[2][4]=-1;
  pl[3][0]=2;pl[3][1]=3;pl[3][2]=7;pl[3][3]=6;pl[3][4]=-1;
  pl[4][0]=0;pl[4][1]=3;pl[4][2]=7;pl[4][3]=4;pl[4][4]=-1;
  pl[5][0]=1;pl[5][1]=2;pl[5][2]=6;pl[5][3]=5;pl[5][4]=-1; 
  do
  {
      clean();
      preobr();
      vpi();
      draw();
      blit(bmp,screen,0,0,0,0,W,H);
  }
  while(!key[KEY_ESC]);        
  
  system("PAUSE");  
  return 0;
}
END_OF_MAIN();
 
void ed()
{
    for(i=0;i<4;i++)
      for(j=0;j<4;j++)
      if(i==j)
      p[i][j]=1;
      else
      p[i][j]=0;
}
 
float va()
{
    return((cd[pl[i][0]][1]-cd[pl[i][1]][1])*(cd[pl[i][2]][2]-cd[pl[i][0]][2])-(cd[pl[i][0]][2]-cd[pl[i][1]][2])*(cd[pl[i][2]][1]-cd[pl[i][0]][1]));
}
 
float vb()
{
    return(-((cd[pl[i][0]][0]-cd[pl[i][1]][0])*(cd[pl[i][2]][2]-cd[pl[i][0]][2])-(cd[pl[i][0]][2]-cd[pl[i][1]][2])*(cd[pl[i][2]][0]-cd[pl[i][0]][0])));
} 
 
float vc()
{
    return((cd[pl[i][0]][0]-cd[pl[i][1]][0])*(cd[pl[i][2]][1]-cd[pl[i][0]][1])-(cd[pl[i][0]][1]-cd[pl[i][1]][1])*(cd[pl[i][2]][0]-cd[pl[i][0]][0]));
} 
 
float vd()
{
    return(plosk[i][0]*(-cd[pl[i][0]][0])+plosk[i][1]*(-cd[pl[i][0]][1])+plosk[i][2]*(-cd[pl[i][0]][2]));
}
 
void vpi()
{
    for(i=0;i<kg;i++)
      {
          plosk[i][0]=va();
          plosk[i][1]=vb();
          plosk[i][2]=vc();
          plosk[i][3]=vd();
      }
    for(i=0;i<kg;i++)
    {
        if((plosk[i][0]+plosk[i][1]+plosk[i][2]+plosk[i][3])<0)
        {
            plosk[i][0]=-plosk[i][0];
            plosk[i][1]=-plosk[i][1];
            plosk[i][2]=-plosk[i][2];
            plosk[i][3]=-plosk[i][3];
        }
    }
}
 
void draw()
{
    for(i=0;i<kg;i++)
    {
        if(plosk[i][2]<0)
        {
           j=0;
           while(pl[i][j+1]>=0)
           {
               line(bmp,(int)cd[pl[i][j]][0]+nx,(int)cd[pl[i][j]][1]+ny,(int)cd[pl[i][j+1]][0]+nx,(int)cd[pl[i][j+1]][1]+ny,green);
               sx=sx+(int)cd[pl[i][j]][0]+nx;
               sy=sy+(int)cd[pl[i][j]][1]+ny;
               j++;
           }
           line(bmp,(int)cd[pl[i][0]][0]+nx,(int)cd[pl[i][0]][1]+ny,(int)cd[pl[i][j]][0]+nx,(int)cd[pl[i][j]][1]+ny,green);
           sx=sx+(int)cd[pl[i][j]][0]+nx;
           sy=sy+(int)cd[pl[i][j]][1]+ny;
           floodfill(bmp,sx/(j+1),sy/(j+1),sx/(j+1)+sy/(j+1));
           //floodfill(bmp,50,150,50);
           sx=0;
           sy=0;
           if(getpixel(bmp,0,0)!=0)
           floodfill(bmp,0,0,0);
       }
   }
}
 
void clean()
{
    for(i=0;i<kg;i++)
    {
        if(plosk[i][2]<0)
        {
           j=0;
           while(pl[i][j+1]>0)
           {
               line(bmp,(int)cd[pl[i][j]][0]+nx,(int)cd[pl[i][j]][1]+ny,(int)cd[pl[i][j+1]][0]+nx,(int)cd[pl[i][j+1]][1]+ny,0);
               sx=sx+(int)cd[pl[i][j]][0]+nx;
               sy=sy+(int)cd[pl[i][j]][1]+ny;
               j++;
           }
           line(bmp,(int)cd[pl[i][0]][0]+nx,(int)cd[pl[i][0]][1]+ny,(int)cd[pl[i][j]][0]+nx,(int)cd[pl[i][j]][1]+ny,0);
           sx=sx+(int)cd[pl[i][j]][0]+nx;
           sy=sy+(int)cd[pl[i][j]][1]+ny;
           floodfill(bmp,sx/(j+1),sy/(j+1),0);     
           sx=0;
           sy=0;
       }
   }
}
 
void preobr()
{
      if(key[KEY_UP]&&ny>0)
      ny=ny-1;
      if(key[KEY_DOWN]&&ny<H)
      ny=ny+1;
      if(key[KEY_LEFT]&&nx>0)
      nx=nx-1;
      if(key[KEY_RIGHT]&&nx<W)
      nx=nx+1;    
      if(key[KEY_A])
      {
          ed();
          for(i=0;i<3;i++)
          p[i][i]=1.1;
          matrp();
      }
      if(key[KEY_Z])
      {
          ed();
          for(i=0;i<3;i++)
          p[i][i]=0.9;
          matrp();
      }        
      if(key[KEY_W])
      {
         ed();
         p[0][0]=0.9962;
         p[0][1]=0.0872;
         p[1][0]=-p[0][1];
         p[1][1]=p[0][0];
         matrp();
     } 
     if(key[KEY_E])
     {
         ed();
         p[0][0]=0.9962;
         p[0][1]=-0.0872;
         p[1][0]=-p[0][1];
         p[1][1]=p[0][0];
         matrp();
     }
     if(key[KEY_S])
     {
         ed();
         p[1][1]=0.9962;
         p[1][2]=0.0872;
         p[2][1]=-p[1][2];
         p[2][2]=p[1][1];
         matrp();
     }
     if(key[KEY_D])
     {
         ed();
         p[1][1]=0.9962;
         p[1][2]=-0.0872;
         p[2][1]=-p[1][2];
         p[2][2]=p[1][1];
         matrp();
     }
     if(key[KEY_X])
     {
         ed();
         p[0][0]=0.9962;
         p[2][0]=0.0872;
         p[0][2]=-p[2][0];
         p[2][2]=p[0][0];
         matrp();
     }                                          
     if(key[KEY_C])
     {
         ed();
         p[0][0]=0.9962;
         p[2][0]=-0.0872;
         p[0][2]=-p[2][0];
         p[2][2]=p[0][0];
         matrp();
     }                          
}
 
void matrp()
{
    i=0;
    do{
    for(j=0;j<4;j++)
    cd1[i][j]=cd[i][0]*p[0][j]+cd[i][1]*p[1][j]+cd[i][2]*p[2][j]+cd[i][3]*p[3][j];
    i++;
    }while(cd[i-1][3]!=-1);
    i=0;
    do{
    for(j=0;j<4;j++)
    cd[i][j]=cd1[i][j];
    i++;
    }while(cd[i-1][3]!=-1);    
}

//END_OF_MAIN();