• <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 Z[MAX+1];
            int D[MAX+1]; //未知變量位置的變化
            int n;
            int e;

            void input()
            {
                
            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);
            }


            void SwapE(double a,double b)    //swap elements
            {
                
            double T;
                T
            =a;
                a
            =b;
                b
            =T;
            }


            void  SwapR(int k,int kmi)
            {
                
            int j;
                
            for(j=k;j<=n;j++)
                    SwapE(A[k][j],A[kmi][j]);
            }


            void SwapC(int k,int kmj)
            {
                
            int i;
                
            for(i=k;i<=n;i++)
                    SwapE(A[i][k],A[i][kmj]);
            }


            void AllGaussianElimination()
            {
                
            int kmi,kmj;//the i and j of the max(abs) element when k
                int i,j,k;
                
            double T;

                
            for(i=1;i<=n;i++)//初始化位置變量位置
                    D[i]=i;

                
            for(k=1;k<=n-1;k++)
                
            {
                    
            //選主元
                    T=0;
                    
            for(i=k;i<=n;i++)
                        
            for(j=k;j<=n;j++)
                            
            if(fabs(A[i][j])>T) { T=fabs(A[i][j]); kmi=i;kmj=j;}

                    
            if(T<=e) {printf("Error!\n"); return ;}

                    
            if(kmi!=k) { SwapR(k,kmi);  SwapE(B[k],B[kmi]); }
                    
            if(kmj!=k) { SwapC(k,kmj);  SwapE(D[k],D[kmj]); }
                    
            //消元
                    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(A[n][n]<=e) {printf("Error!\n");return ;}
                    Z[n]
            =B[n]/A[n][n];
                    
                    
            double S_Aij_Zj;
                    
            for(i=n-1;i>=1;i--)
                    
            {
                        S_Aij_Zj
            =0;
                        
            for(j=i+1;j<=n;j++)
                            S_Aij_Zj
            +=A[i][j]*Z[j];

                        Z[i]
            =(B[i]-S_Aij_Zj)/A[i][i];
                    }


                    
            for(j=1;j<=n;j++)
                        X[D[j]]
            =Z[j];
                }

            }


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


            int main()
            {
                input();
                AllGaussianElimination();
                print(X);
                system(
            "pause");
                
            return 0;
            }
            posted on 2009-05-16 16:37 wyiu 閱讀(805) 評論(0)  編輯 收藏 引用 所屬分類: 數值分析
            久久se精品一区精品二区| 亚洲午夜久久久久久久久久| 久久不射电影网| 久久久91人妻无码精品蜜桃HD | 久久最近最新中文字幕大全 | 久久精品国产亚洲AV高清热| 国产情侣久久久久aⅴ免费| 99久久精品无码一区二区毛片 | 香蕉久久夜色精品国产2020 | 久久香蕉国产线看观看乱码| 久久精品无码专区免费| 人妻精品久久久久中文字幕69| 中文字幕久久欲求不满| 国内高清久久久久久| 久久国产精品免费一区| 人妻无码久久一区二区三区免费 | 曰曰摸天天摸人人看久久久| 国产香蕉久久精品综合网| 久久婷婷综合中文字幕| 99久久无色码中文字幕人妻| 久久国产精品99精品国产987| 色综合久久88色综合天天 | 91精品国产综合久久香蕉| 国产激情久久久久久熟女老人| 国产福利电影一区二区三区久久老子无码午夜伦不 | 久久久久久国产a免费观看黄色大片 | 久久精品国产久精国产果冻传媒| 色综合合久久天天综合绕视看| 中文字幕乱码久久午夜| 亚洲色欲久久久久综合网 | 亚洲综合久久夜AV | 久久国产精品一区| 国产一区二区三精品久久久无广告 | 久久99精品久久久久久hb无码| 精品无码久久久久国产动漫3d| 日日狠狠久久偷偷色综合0| 久久99精品久久久久久9蜜桃| 国产亚洲精午夜久久久久久 | 合区精品久久久中文字幕一区| 国产精品成人无码久久久久久| 99久久成人18免费网站|