Матрицы. Ступенчатый вид - C++ - Ответ 5593576

#include <string>
#include <iostream>
using namespace std;
 
int prog(int x)
{return (x-1)*x/2;}
 
void pokaz(int n,int l,double a[][20] )
{for(int i=0;i<n;i++){for(int j=0;j<l;j++){cout<<a[i][j]<<"    ";}cout<<endl<<"                   ";}}
 
void pokaz_1(int n,int l,double a[][20] )
{for(int i=0;i<n;i++){for(int j=0;j<l;j++){cout<<a[i][j]<<"    ";}cout<<endl;}}
 
void vvod(int n,int l,double a[][20] )
{for(int i=0;i<n;i++){for(int j=0;j<l;j++){cin>>a[i][j];}cout<<endl;}}
 
int def(int n,int l,double a[][20])
{int t=0,k=0;
for(int i=0;i<n;i++)
{t=0;
for(int j=0;j<l;j++){if(a[i][j]==0){t++;}}
    if(t==l){k++;}
}
 
if(l-n>=0){k=l-n+k;}
if(l-n<0){if(n-k>=l){k=0;}else{k-n+l;}}
return k;
} 
 
double my_det(int n,double a[][20])
{double x,q,p=1;
    for(int k=0;k<prog(n);k++)
{
for(int j=1+k;j<=n-1;j++)
{
x=a[k][k]/a[j][k];
p=p/(-x);
for(int i=0+k;i<n;i++){a[j][i]=a[k][i]-(x*a[j][i]);}
 
}
}
 
 
 
for(int i=0;i<n;i++){p=p*a[i][i];}
 
 
 
return p;
}
 
void chang(int l,int k,int j,double a[][20])
{
    for(int i=0;i<l;i++)
    {
        double x = a[k][i];
    a[k][i]=a[j][i];
    a[j][i]=x;
    }
}
 
void stop()
{system("pause");}
 
void clear()
{system("cir");}
 
 
void main()
{
    /*double c[20][20];
    vvod(3,3,c);
    pokaz_1(3,3,c);
    chang*/
 
 
 
setlocale(0,"");
newmatrix:
string b;
int n,l,g;
double a[20][20];
double x,q,p=1;
cout<<endl<<"Введите кол-во строк : ";     cin>>n;  cout<<endl;
cout<<endl<<"Введите кол-во столбцов : ";  cin>>l;  cout<<endl;
 
 
cout<<"Введите матрицу : "<<endl;vvod(n,l,a);
cout<<endl<<endl<<"Посмотреть шаги построения ? Ввести 'yes' или 'no' ! : ";
cin>>b;
system("cls");
 
cout<<"Исходная матрица : ";pokaz(n,l,a);
cout<<endl;
 
//if(a[0][0]!=0){q=a[0][0];for(int i=0;i<n;i++){for(int j=0;j<l;j++){a[i][j]=a[i][j]/q;}p=p*a[0][0];}}
 
 
for(int k=0;k<prog(n);k++)
{
for(int j=1+k;j<=n-1;j++)
{
am:if(j<n)
{
    if(a[j][k]!=0)
    {
        if(a[k][k]!=0)
        {  x=a[k][k]/a[j][k];p=p/(-x);}
        
        if(a[k][k]==0)
        {
            for(g=k;g<n;g++){if(a[g][k]!=0){break;}}
            chang(l,k,g,a);p=-p;if(b[0]=='y'){cout<<"Замена строки "<<k+1<<" на строку "<<g+1<<endl;
        }
        }
    }
    
    
    if(a[j][k]==0){j++;goto am;}
    if((b[0]=='y') && (j+1<=n)){cout<<"Из строки "<<k+1<< " вычитаем строку "<<j+1<<" умноженную на "<<a[k][k]<<"/"<<a[j][k]<<",записываем в строку "<<j+1<<endl<<endl;}
    for(int i=0+k;i<l;i++){a[j][i]=a[k][i]-(x*a[j][i]);}
    if(b[0]=='y'){pokaz_1(n,l,a);cout<<endl<<endl;}
}
 
}}
 
 
 
for(int i=0;i<n;i++){p=p*a[i][i];}
 
 
 
cout<<endl<<"Ранг матрицы         : "<<l-def(n,l,a);
cout<<endl<<"Дефект матрицы       : "<<def(n,l,a)<<endl;
cout<<endl<<"Детерминант матрицы : ";
 
 
    if(n!=l){cout<<" МАТРИЦА НЕ КВАДРАТНАЯ "<<endl;}
    if (n==l){cout<<p;}
cout<<endl<<endl<<"Ступенчатый вид матрицы : ";    for(int i=0;i<n;i++){for(int j=0;j<l;j++){cout<<a[i][j]<<"      ";}cout<<endl;cout<<"\t"<<"\t"<<"\t"<<"  ";}
 
 
cout<<endl<<"________________________________________________________________________________"<<endl;
 
 
 
 
 
 
 
 
 
goto newmatrix;
}