Решение СЛАУ методом Гаусса

#include <stdio.h>
#include <conio.h>
#include <tchar.h>
#include <iostream>
#include <stdlib.h>
#include <time.h>
#define eps 0.0000000001
 
class CMatrix
{private: int m;
int n;
double*p;
public: CMatrix (){};
CMatrix(int,int);
friend CMatrix operator+(CMatrix&,CMatrix&);
friend CMatrix operator*(CMatrix&,CMatrix&);
friend CMatrix gauss(CMatrix&,CMatrix&);
void input();
void output(int,int);
void output1(int,int);
void input1();
double det();
};
 
void gauss();
void adddouble ();
void multdouble ();
void det();
 
void main ()
{int k;
clrscr();
cout<<» add 1\n multiply 2\n schutat’ opredelitel 3\n metod gaussa 4\n»<>k;
if(k==1) adddouble();
else if (k==2) multdouble();
else if (k==3) det();
 
else if (k==4) gauss();
}
 
CMatrix::CMatrix(int a,int b)
{m=a;
n=b;
p=new double [m*n];
}
 
void CMatrix::output(int x, int y)
{int i;
int j;
 
for(i=0;i<m;i++)
{
for(j=0;j<n;j++)
{
gotoxy(5*j+x,2*i+y);
cout<<p[i*n+j];
}
}
}
 
void CMatrix::output1(int x, int y)
{int i;
int j;
 
for(i=0;i<m;i++)
{
for(j=1;j<n;j++)
{
gotoxy(5*j+x,2*i+y);
cout<<p[i*n+j];
}
}
}
 
void det()
{int n;
double d;
cout<<»vvedite razmer matrici n» <>n;
CMatrix A(n,n);
clrscr();
A.input();
A.output(31,9);
d=A.det();
gotoxy(28,15);
cout<<»opredel raven «<<d;
getch();
}
 
void adddouble ()
{int m,n;
clrscr();
cout<<»vvedite razmeri matriz m i n»<>m>>n;
CMatrix A(m,n);
CMatrix B(m,n);
CMatrix C(m,n);
A.input();
B.input();
clrscr();
C=A+B;
getch();
}
 
void multdouble ()
{int m,n,l;
cout<<»vvedite m, n, l»<>m>>n>>l;
CMatrix A(m,n);
CMatrix B(n,l);
CMatrix C(m,l);
A.input();
B.input();
clrscr();
C=A*B;
A.output(1,1);
B.output(41,1);
C.output(21,9);
getch();
}
 
void CMatrix::input()
{int i, j;
cout<<»vvedite info»<<endl;
for(i=0;i<m;i++)
for(j=0;j>p[i*n+j];
 
}
 
void CMatrix::input1()
{int i, j;
cout<<»vvedite info»<<endl;
for(i=1;i<m;i++)
for(j=1;j>p[i*n+j];
 
}
 
CMatrix operator+(CMatrix&obj1,CMatrix&obj2)
{int i,j;
CMatrix temp(obj1.m,obj1.n);
for(i=0;i<obj1.m;i++)
{
 
for(j=0;j<obj1.n;j++)
{
temp.p[i*obj1.n+j]=obj1.p[i*obj1.n+j]+obj2.p[i*obj1.n+j];
gotoxy(5*j+1,2*i+1);
cout<<obj1.p[i*obj1.n+j];
gotoxy(5*j+41,2*i+1);
cout<<obj2.p[i*obj1.n+j];
gotoxy(5*j+21,2*i+8);
cout<<temp.p[i*obj1.n+j];
}
}
return temp;
 
}
 
CMatrix operator*(CMatrix&obj1,CMatrix&obj2)
{
int i,j,k;
double s;
CMatrix temp(obj1.m,obj2.n);
for(i=0;i<obj1.m;i++)
{
for(j=0;j<obj2.n;j++)
{
s=0;
for(k=0;k2)
 
{
CMatrix OBR(n-1,n-1);
for(i=0;i<n;i++)
{ for(j=1;j<n;j++)
{for (k=0;k<n;k++)
if(ki)
OBR.p[(j-1)*(n-1)+(k-1)]=p[j*n+k];
}
znak*=-1;
d+=(p[i]*OBR.det()*znak);
 
}
 
}
return d;
}
 
void gauss()
{
int n;
cout<<»vvedite n «<>n;
CMatrix a (n,n);
CMatrix b (n,1);
CMatrix x (1,n);
a.input1();
b.input1();
x=gauss(a,b);
gotoxy (28,15);
cout<<»otvet = «;
x.output1(38,15);
getch;
 
}
 
CMatrix gauss(CMatrix&a,CMatrix&b)
{
CMatrix x(1,a.n);
int i;
int j;
int k;
int n=a.n;
for (k=1;k<n;k++)//nomer obnul stolbcha
{
/*if(fabs(a.p+k*(a.n+1)+k)<eps )
{
 
}*/
for(i=k;i<n;i++)//i nomer obr str
{ clrscr();
a.output1 (1,1);
b.output1 (37,1);
 
for(j=k;j<n;j++)// nomer stolbcha v obr stroke
{
a.p[i*n+j]=a.p[i*n+j]-a.p[i*n+k]/a.p [k*n+k] * a.p[k*n+j];
 
}
 
b.p[i]=b.p[i]-a.p[i*n+k]/a.p[k*n+k]*b.p[k];
a.output1(40,1);
b.output1(70,1);
getch;
}
 
}
return x;
}