• <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            posts - 100,  comments - 15,  trackbacks - 0
            #include<iostream>
            #include
            <math.h>
            using namespace std;
            #define MAX 100

            double A[MAX+1][MAX+1];
            double B[MAX+1];
            double X[MAX+1];
            double e;
            int n;

            void ColGaussianElimination()
            {
                
            int i,j,k,kmi;
                
            double T;
                
            for(k=1;k<=n-1;k++)
                
            {
                    
            //選主元
                    T=0;
                    
            for(i=k;i<=n;i++)
                        
            if( fabs(A[i][k])> T ){ T=A[i][k];kmi=i;}
                    
            if( T<=e) { printf("Error!\n"); return ;}
                    
            if(kmi!=k)
                    
            {
                        T
            =B[k];B[k]=B[kmi];B[kmi]=T; //swap B[k] and B[kmi]
                        for(j=k;j<=n;j++//swap row kmi and k of A
                        {
                            T
            =A[k][j];
                            A[k][j]
            =A[kmi][j];
                            A[kmi][j]
            =T;
                        }

                    }

                    
            //消元
                    for(i=k+1;i<=n;i++)
                    
            {
                        T
            =A[i][k]/A[k][k];
                        B[i]
            -=T*B[k];

                        
            for(j=k;j<=n;j++)
                            A[i][j]
            -=T*A[k][j];
                    }

                }

                
            //回代
                if( fabs(A[n][n])<=e ) { printf("Error!\n"); return ;}

                X[n]
            =B[n]/A[n][n];

                
            double S_Aij_Xj;
                
            for(i=n-1;i>=1;i--)
                
            {
                    S_Aij_Xj
            =0;
                    
            for(j=i+1;j<=n;j++)
                        S_Aij_Xj
            +=A[i][j]*X[j];

                    X[i]
            =(B[i]-S_Aij_Xj)/A[i][i];
                }

            }

            void print(double X[])
            {
                
            int i;
                printf(
            "X[]:\n");
                
            for(i=1;i<=n;i++)
                    printf(
            "%f\n",X[i]);
            }


            int main()
            {
                
            int i,j;

                printf(
            "n:");
                scanf(
            "%d",&n);

                printf(
            "A[][]:\n");
                
            for(i=1;i<=n;i++)
                    
            for(j=1;j<=n;j++)
                        scanf(
            "%lf",&A[i][j]);

                printf(
            "B[]:\n");
                
            for(i=1;i<=n;i++)
                    scanf(
            "%lf",&B[i]);

                printf(
            "e:");
                scanf(
            "%lf",&e);

                ColGaussianElimination();
                print(X);
                system(
            "pause");
                
            return 0;
            }


            posted on 2009-05-16 16:34 wyiu 閱讀(616) 評(píng)論(0)  編輯 收藏 引用 所屬分類: 數(shù)值分析
            怡红院日本一道日本久久| 精品国产一区二区三区久久久狼| 97精品久久天干天天天按摩| 久久精品国产亚洲欧美| 久久93精品国产91久久综合| 东方aⅴ免费观看久久av| 色婷婷综合久久久久中文| 99久久国产综合精品成人影院| 精品久久久久久无码国产| 伊人久久大香线焦AV综合影院| 丰满少妇高潮惨叫久久久| 亚洲国产成人乱码精品女人久久久不卡| 久久久久久免费视频| 国内精品久久久久| 久久婷婷五月综合97色直播| 久久精品无码一区二区三区| 久久精品无码一区二区WWW| 国产午夜精品理论片久久| 人妻无码αv中文字幕久久琪琪布| 国产亚洲精午夜久久久久久| 久久久免费精品re6| 久久久久亚洲AV无码观看| 欧美久久久久久精选9999| 精品久久久久久综合日本| 亚洲AV成人无码久久精品老人| 久久人人爽人人澡人人高潮AV | 国内精品伊人久久久久777| 久久91精品国产91久久户| 亚洲日韩中文无码久久| 国产精品成人久久久| 亚洲精品无码久久久久AV麻豆| 精品久久久久一区二区三区 | 亚洲国产综合久久天堂| 7国产欧美日韩综合天堂中文久久久久 | 久久99热这里只频精品6| 久久国产乱子伦精品免费午夜| 久久九九亚洲精品| 国产国产成人久久精品| 91精品婷婷国产综合久久| 国产精品日韩深夜福利久久| 国产精品内射久久久久欢欢|