#include<iostream.h>
#include<iomanip.h>
#include<conio.h>
#include<stdio.h>
#include<stdlib.h>
#include<math.h>
void show(float **a,int n);
float det(float **b,int n);
void minor(float **a,int o,int i,int j);
void main()
{
clrscr();
int o,i,j;
float **a,ans;
cout<<"Enter The Order Of Matrix:";
cin>>o;
for(i=0;i<o;i++)
a[i]=new float[o];
for(i=0;i<o;i++)
for(j=0;j<o;j++)
cin>>a[i][j];
ans=det(a,o);
cout<<ans;
getch();
}
void show(float **a,int n)
{
cout<<"\n";
for(int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
{cout<<setw(8)<<a[i][j]<<" ";}
cout<<"\n\n";
}
}
float **b;
float det(float **c,int n)
{
float ans=0;
float *d;
d=new float[n];
for(int i=0;i<n;i++)
d[i]=c[0][i];
if(n>2)
{
for(int i=0;i<n;i++)
{
minor(c,n,0,i);
if(i%2==0)
ans+=d[i]*det(b,n-1);
else
ans-=d[i]*det(b,n-1); 
}
}
if(n==2)
{
return(c[0][0]*c[1][1]-c[0][1]*c[1][0]);
}
return ans;
}
void minor(float **e,int o,int i,int j)
{
int k1=0,l1;
for(int k=0;k<o;k++)
{
l1=0;
for(int l=0;l<o;l++)
{
if(k==i)
{k1--;break;}
if(l==j)
continue;
[b][color=red]b[k1][l1]=e[k][l];:-O [/color][/b]
l1++;
}
k1++;
}
}

a = new float*[o]; // allocate first dimension
for(int i = 0; i < o; i++)
{
    a[i] = new float(o);
    // initialize each element
    for(int j = 0; j < o; j++)
         a[i][j] = 0.0F;
}