Метод Гаусса - C++ - Ответ 1536648

#include <iostream.h>
#include <math.h>
#include <iomanip.h>
#include <windows.h>
 
int main()
{
    int n,i,j,k,p,m,z,y,l;
    char f;
    double **a,**v,*b,*x,*pr,s,c,h,r,*u,*pog,e;
y=1;
 
    
    
               while (y>0)
{
 l=1;
 y=0;
 
      cout<<"Vvedite razmer:";
  cin>>n;
  if (n<2) {cout<<"ne imeet smisla!"<<endl; l=0;}
  a = new double *[n+1];
  v = new double *[n+1];
    for (i=0;i<=n;i++)
       a[i]= new double[n+1];
    for (i=0;i<=n;i++)
       v[i]= new double[n+1];
 b= new double[n+1];
 x= new double[n+1];
 u= new double[n+1];
 pog= new double[n+1];
 
 if (n>1)     
 {
 for(i=1;i<=n;i++)
  for(j=1;j<=n;j++)
   {cout<<"a["<<i<<"]["<<j<<"]="; cin>>a[i][j];v[i][j]=a[i][j];}
 
 
 for(i=1;i<=n;i++)
  {cout<<"Svobpdnii 4len "<<i<<"="; cin>>b[i]; u[i]=b[i];}
 
  cout<<"Vid sistemi:"<<endl;
for(i=1;i<=n;i++)
   for(j=1;j<=n;j++)
   {cout<<setw(3)<<a[i][j]<<" "; if(j==n) cout<<"= "<<setw(3)<<b[i]<<endl;}
 
 
 for(k=1;k<=n-1;k++)
  {
      p=k;
      for(m=k+1;m<=n;m++)
          if(fabs(a[p][k])<fabs(a[m][k])) p=m;
     //usaem ykazateli^_^
          pr=a[k];
          a[k]=a[p];
          a[p]=pr;
     r=b[k];b[k]=b[p];b[p]=r;
    //Vremya gausa!!
     for(m=k+1;m<=n;m++)
     {
         c=a[m][k]/a[k][k];
        b[m]=b[m]-c*b[k];
         for(i=k;i<=n;i++)
             a[m][i]-=c*a[k][i];
     }
  }
 
               }
 
//proverka diagonali
   z=1;
    for (i=1;i<=n;i++)
    {if (a[i][z]==0)  {cout<<"Vashe sistema ne imeet rewenii libo ih beskonechno mnogo"<<endl;l=0;i=n+1; }z++;}
        
//check
   if (l!=0)
     {
 cout<<"Posle preobrazovaniya:"<<endl;
 for(i=1;i<=n;i++)
    for(j=1;j<=n;j++)
    {cout<<setw(3)<<a[i][j]<<" "; if(j==n) cout<<"= "<<setw(3)<<b[i]<<endl;}
  
 //idem nazad!!!! 
  x[n]=b[n]/a[n][n];
   for(k=n-1;k>=1;k--)
   {
      s=0;
      for(i=k+1;i<=n;i++)
           s+=a[k][i]*x[i];
      x[k]=(b[k]-s)/a[k][k];
   }
  cout<<"----------------------------------------"<<endl;
  cout<<"Korni:"<<endl;
  for(i=1;i<=n;i++)  cout<<"x["<<i<<"]="<<x[i]<<endl;
 
 for(i=1;i<=n;i++)
  {   h=0;
    for(j=1;j<=n;j++) h+=x[j]*v[i][j];
    pog[i]=fabs(u[i]-h);
  }  
 e=pog[1];
 for(i=1;i<n;i++)
   if (pog[i+1]<pog[i]) e=pog[i];
 cout<<"Obsolutnaya pogreshnost'="<<e<<endl;
 cout<<"----------------------------------------"<<endl;
   }
 
for(i=0;i<n+1;i++)
  delete[]v[i];
for(i=0;i<n+1;i++)
  delete[]a[i];
delete []a;
delete []v;
delete []b;
delete []u;
delete []pog;
 
cout<<"  -Gelaete li rewit' drygyu sistemy(y/n)?"<<endl;
cin>>f;
          if (f=='y') 
          {  y=1; system("cls"); if(l!=0){ cout<<"Korni prevedyshei sistemi:"<<endl; 
             for(i=1;i<=n;i++)
                   cout<<"x["<<i<<"]="<<x[i]<<"  "; cout<<endl<<endl;
             cout<<"Obsolutnaya pogreshnost'="<<e<<endl;}}
          else 
          {system("cls");  cout<<"Korni prevedyshei sistemi:"<<endl; 
             for(i=1;i<=n;i++)
                   cout<<"x["<<i<<"]="<<x[i]<<"  "; cout<<endl<<endl;
          cout<<"Obsolutnaya pogreshnost'="<<e<<endl;
 
    }
delete []x;
}
 
return 0;
}