• <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ù)值分析 
            

            







            
	    
    




            






            

				

            




    
    







            
	   
	



综合久久久久久中文字幕亚洲国产国产综合一区首
|
婷婷伊人久久大香线蕉AV|
日产精品久久久久久久|
久久精品aⅴ无码中文字字幕不卡|
狠狠精品久久久无码中文字幕|
狠狠色狠狠色综合久久|
亚洲一区二区三区日本久久九|
青青青青久久精品国产|
久久久久久久久波多野高潮|
av国内精品久久久久影院|
久久国产三级无码一区二区|
香蕉久久夜色精品升级完成|
久久久久久极精品久久久|
久久久av波多野一区二区|
国产精品午夜久久|
色综合久久久久久久久五月|
久久99精品久久久久久9蜜桃|
亚洲人成网亚洲欧洲无码久久
|
伊人久久大香线蕉综合热线|
人妻少妇久久中文字幕一区二区|
国产午夜精品久久久久九九|
久久亚洲AV成人出白浆无码国产|
亚洲欧洲精品成人久久曰影片|
国产精品国色综合久久|
久久成人小视频|
亚洲国产成人精品女人久久久|
伊人丁香狠狠色综合久久|
久久精品国产精品亚洲毛片
|
亚洲人成无码网站久久99热国产|
国产99久久久久久免费看
|
漂亮人妻被中出中文字幕久久|
久久精品中文字幕久久|
精品国产VA久久久久久久冰|
97精品依人久久久大香线蕉97|
香港aa三级久久三级老师2021国产三级精品三级在
|
精品国产乱码久久久久软件
|
青青青国产精品国产精品久久久久|
99精品国产99久久久久久97|
亚洲中文字幕无码久久综合网|
日本五月天婷久久网站|
亚洲人成精品久久久久|