C++. Крамер, Гаусс, БД, текст. лабораторная за 1 семестр / Крамер

#include "iostream.h"
#include "stdio.h"
#include "conio.h"

void diagonal(float**,int);

void treugolnik(float**,int);

void printmatrix(float**,int);

float det(float**,int);

void main()
{
cout<<" \t\t Metod KPAMEPA"<<endl;
cout<<"Chislo uravnenii: ";
int n=0;
cin>>n;

float** a = new float*[n];

for(int i=0;i<n;i++)
a[i]=new float[n];
float** t=new float*[n];
for(i=0;i<n;i++)
t[i]=new float[n];

float* b=new float[n];



cout<<"Vvedite matritsu koeffitsientov: "<<endl;
for(i=0;i<n;i++)
for(int j=0; j<n; j++)
{
cout<<"A["<<i<<','<<j<<"]=";
cin>>a[i][j];
t[i][j]=a[i][j];
}

cout<<"Vvedite vector svobodnih chlenov: "<<endl;

for(i=0;i<n;i++)
{
cout<<"B["<<i<<"]=";
cin>>b[i];

}

float Delta=det(t,n);
if(Delta==0)
{
cout<<"Dennaya sistema ne reshaetsya metodom KPAMEPA:)";
return;
}

for(int j=0;j<n;j++)
{
for(int i=0;i<n;i++)
for(int k=0;k<n;k++)
t[i][k]=a[i][k];


for(i=0;i<n;i++)
t[i][j]=b[i];

float x=det(t,n)/Delta;
printf(" %s %d %s %3.5f \n","X",j, "=", x);
cout<<endl;
}

for(i=0;i<n;i++)
delete *(t+i);
delete []t;

for(i=0;i<n;i++)
delete *(a+i);
delete []a;

delete []b;
getch();
}

void diagonal(float** a,int n)
{
for(int j=0;j<n;j++)
if(*(*(a+j)+j)==0)
for(int i=j+1;i<n;i++)
if(*(*(a+i)+j)!=0)
for(int k=0;k<n;k++)
{ float t=*(*(a+j)+k);
*(*(a+j)+k)=*(*(a+i)+k);
*(*(a+i)+k)=t;
}
if(*(*(a+n-1)+n-1)==0)
for(int i=0;i<n;i++)
if( (*(*(a+n-1)+i)!=0)&&(*(*(a+i)+n-1)!=0))
{
for(int k=0;k<n;k++)
{float t=*(*(a+i)+k);
*(*(a+i)+k)=*(*(a+n-1)+k);
*(*(a+n-1)+k)=t;
}
return;
}

}

void treugolnik(float** a,int n)
{
for(int j=0;j<n;j++)
{
for(int i=j+1;i<n;i++)
if(*(*(a+i)+j)!=0)
{ float q=(*(*(a+i)+j))/(*(*(a+j)+j));
for(int k=0;k<n;k++)
(*(*(a+i)+k))-=*(*(a+j)+k)*q;
}

diagonal(a,n);
}
}

void printmatrix(float** a,int n)
{
for( int i=0;i<n;i++)
{
for(int j=0;j<n;j++)
printf("%3.2f \t",*(*(a+i)+j));
printf("\n");

}
}
float det(float** a,int n)
{
diagonal(a,n);
treugolnik(a,n);
float result=*(*(a+0)+0);
for(int i=1;i<n;i++)
result*=(*(*(a+i)+i));
return result;

}