Advertisement
Jgug

Gayss_01

Mar 12th, 2012
415
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C++ 3.44 KB | None | 0 0
  1. #include <math.h>
  2. #include <iostream.h>
  3. #include <stdio.h>
  4.  
  5. int main()
  6. {
  7.  
  8.     double A[101][101];
  9.     double B[101][1];
  10.     double X[101][1];
  11.     double koeff, koeff2, temp1, temp2, temp3, temp4;
  12.     int size, num=0;
  13.  
  14. //-------------------------------------------------|
  15.  
  16.     cout<<"Vvedite razmer matricsi sistemi: ";
  17.     cin>>size;
  18.  
  19.     cout<<"Vvedite element matricsi A: "<<endl;
  20.     for ( int a1=0; a1<size; a1++)
  21.     {
  22.         for ( int a2=0; a2<size; a2++)
  23.         {
  24.             cin>>A[a1][a2];
  25.         }
  26.     }
  27.  
  28.     cout<<"Vvedennaya matricsa: "<<endl;
  29.     for ( int a3=0; a3<size; a3++)
  30.     {
  31.         for ( int a4=0; a4<size; a4++)
  32.         {
  33.             cout<<A[a3][a4]<<" ";
  34.         }
  35.         cout<<endl;
  36.     }
  37. //-------------------------------------------------|
  38.  
  39.     cout<<"Vvedite stolbecs svobodnich chlenov: ";
  40.     for ( int b1=0; b1<size; b1++)
  41.     {
  42.         cin>>B[b1][0];
  43.     }
  44.  
  45.     cout<<"Vvedennii stolbecs svobodnich chlenov: "<<endl;
  46.     for ( int b2=0; b2<size; b2++)
  47.     {
  48.         cout<<B[b2][0]<<endl;
  49.     }
  50.  
  51. //-------------------------------------------------|
  52.    
  53.     double max=0;
  54.     for ( int c1=0; c1<size; c1++)
  55.     {
  56.         if ( fabs(A[c1][0])>max)
  57.         {
  58.             max=A[c1][0];
  59.             num=c1;
  60.         }
  61.     }
  62.  
  63.     cout<<max<<" "<<num+1<<endl;
  64.  
  65. //-------------------------------------------------|
  66.    
  67.     for ( int d1=0; d1<size; d1++)
  68.     {
  69.         temp1=A[0][d1];
  70.         A[0][d1]=A[num][d1];
  71.         A[num][d1]=temp1;
  72.     }
  73.  
  74.     temp2=B[0][0];
  75.     B[0][0]=B[num][0];
  76.     B[num][0]=temp2;
  77.  
  78.     cout<<"Vvedennaya matricsa: "<<endl;
  79.     for ( int d2=0; d2<size; d2++)
  80.     {
  81.         for ( int d3=0; d3<size; d3++)
  82.         {
  83.             cout<<A[d2][d3]<<" ";
  84.         }
  85.         cout<<endl;
  86.     }
  87.  
  88.     for ( int d4=0; d4<size; d4++)
  89.     {
  90.         cout<<B[d4][0]<<endl;
  91.     }
  92.    
  93. //-------------------------------------------------|
  94.  
  95.     for ( int e1=0; e1<size; e1++)
  96.     {
  97.         for ( int e2=e1; e2<size-1; e2++)
  98.         {
  99.             koeff=A[e2+1][e1]/A[e1][e1];
  100.             B[e2+1][0]=B[e2+1][0]-B[e1][0]*koeff;
  101.             for ( int e3=e1; e3<size; e3++)
  102.             {
  103.                 A[e2+1][e3]=A[e2+1][e3]-A[e1][e3]*koeff;
  104.             }
  105.         }
  106.     }
  107.  
  108.     cout<<"Vvedennaya matricsa: "<<endl;
  109.     for ( int e4=0; e4<size; e4++)
  110.     {
  111.         for ( int e5=0; e5<size; e5++)
  112.         {
  113.             cout<<A[e4][e5]<<" ";
  114.         }
  115.         cout<<endl;
  116.     }
  117.  
  118.     for ( int e6=0; e6<size; e6++)
  119.     {
  120.         cout<<B[e6][0]<<endl;
  121.     }
  122.  
  123. //-------------------------------------------------|
  124. /*
  125.     for ( int f1=size; f1>0; f1--)
  126.     {
  127.         for ( int f2=f1; f2>0+1; f2--)
  128.         {
  129.             koeff=A[f2-1][f1]/A[f1][f1];
  130.             B[f2-1][0]=B[f2-1][0]-B[f1][0]*koeff;
  131.             for ( int f3=f1; f3>0; f3--)
  132.             {
  133.                 A[f2-1][f3]=A[f2-1][f3]-A[f1][f3]*koeff;
  134.             }
  135.         }
  136.     }
  137.    
  138.     cout<<"Vvedennaya matricsa: "<<endl;
  139.     for ( int f4=0; f4<size; f4++)
  140.     {
  141.         for ( int f5=0; f5<size; f5++)
  142.         {
  143.             cout<<A[f4][f5]<<" ";
  144.         }
  145.         cout<<endl;
  146.     }
  147.  
  148.     for ( int f6=0; f6<size; f6++)
  149.     {
  150.         cout<<B[f6][0]<<endl;
  151.     }
  152. */
  153. //-------------------------------------------------|
  154.  
  155.     for ( int f1=size; f1>0; f1--)
  156.     {
  157.         for ( int f2=f1; f2>0; f2--)
  158.         {
  159.             koeff2=A[f2-1][f1]/A[f1][f1];
  160.             B[f2-1][0]=B[f2-1][0]-B[f1][0]*koeff2;
  161.             for ( int f3=f1; f3>0; f3--)
  162.             {
  163.                 A[f2-1][f3]=A[f2-1][f3]-A[f1][f3]*koeff;
  164.             }
  165.         }
  166.     }
  167.  
  168.     cout<<"Vvedennaya matricsa: "<<endl;
  169.     for ( int f4=0; f4<size; f4++)
  170.     {
  171.         for ( int f5=0; f5<size; f5++)
  172.         {
  173.             cout<<A[f4][f5]<<" ";
  174.         }
  175.         cout<<endl;
  176.     }
  177.  
  178.     for ( int f6=0; f6<size; f6++)
  179.     {
  180.         cout<<B[f6][0]<<endl;
  181.     }
  182.  
  183. //-------------------------------------------------|
  184. /*
  185. for(r=1;r<n;r++)
  186.     {
  187.         for(i=1;i<=r;i++)
  188.         {
  189.             t=mas[i-1][r]/mas[r][r];
  190.             B[i-1]=B[i-1]-B[r]*t;
  191.             for(j=r;j<n;j++)
  192.             {
  193.                 mas[i-1][j]=mas[i-1][j]-mas[r][j]*t;
  194.             }
  195.         }
  196.     }
  197. */
  198.  
  199.     return 0;
  200. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement