C / C++ / MFC

public:
	void SetBitmapDataSize(const int x, const int y);

private:
	int m_PixelNumberX;
	int m_PixelNumberY;
	vector<vector<COLORREF> > BitmapData; 

void CBitmapOp::SetBitmapDataSize(const int x, const int y)
{
	m_PixelNumberX = x;
	m_PixelNumberY = y;
	if((x>=0) && (y>=0))
	{		
		BitmapData.resize(m_PixelNumberX); 
		for(int x = 0; x < m_PixelNumberX; x++) 
			BitmapData[x].resize(m_PixelNumberY); 
	}
}

private:
	int m_PixelNumberX;
	int m_PixelNumberY;
	CArray< CArray<COLORREF, COLORREF&>, CArray<COLORREF, COLORREF&> > m_arrBitmapData;

float determinant(float **input,int order)
{

int j,p,q,t,i;
float pr;
float d;
float *c;
//c=memory_alloc_2D1(row,column);
float **answer;
answer=memory_alloc_2D1(order,order);
float **b;
b=memory_alloc_2D1(order,order);

for(j=1;j<=order;j++)
   {       
      int r=1,s=1;
      for(p=1;p<=order;p++)
       	 {
         
           for(q=1;q<=order;q++)
           		 {
            		        if(p!=1&&q!=j)
              				{
   				
		answer[r][s]=input[p][q]; // error as segmentation fault
                				s++;
                				if(s>order-1)
                 				{
                   					r++;
                   					s=1;
                  				}   
              				 }
            		 }
        }

     	for(t=1,pr=1;t<=(1+j);t++)
     		pr=(-1)*pr;
     		c[j]=pr*determinant(b,order-1);
     }
     		for(j=1,d=0.0;j<=order;j++)
     		{
      			 d=d+(input[1][j]*c[j]);
     		}

return(d);

}

float determinantmatrix;
determinantmatrix=determinant(covariance1,waveframesize);