• <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 OrderGaussianElimination()
            {
                            int i,j,k;
                            double T;
                            //消元
                for(k=1;k<=n-1;k++)
                            {
                                if( fabs(A[k][k])<=e ) { printf("Error!\n"); return ;}
                                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);

                OrderGaussianElimination();
                print(X);
                            return 0;
            }

            

             
            

	
            posted on 2009-05-16 16:32 wyiu 閱讀(244) 評論(0)  編輯 收藏 引用  所屬分類: 數值分析 
            

            







            
	    
    




            






            

				

            




    
    







            
	   
	



国产精品久久久久久久久久影院|
中文字幕人妻色偷偷久久|
日本久久久精品中文字幕|
97精品伊人久久久大香线蕉|
国产激情久久久久影院小草|
一个色综合久久|
国产亚洲婷婷香蕉久久精品|
人妻系列无码专区久久五月天|
日本人妻丰满熟妇久久久久久|
99久久精品久久久久久清纯|
成人午夜精品无码区久久|
99久久精品免费看国产一区二区三区|
久久亚洲精品无码aⅴ大香|
久久99精品国产麻豆婷婷|
久久99精品久久久久久动态图|
香蕉久久永久视频|
国产一区二区精品久久凹凸|
亚洲国产精品无码久久一线|
久久99热这里只有精品66|
精品久久香蕉国产线看观看亚洲|
久久久久亚洲AV无码专区首JN|
精品99久久aaa一级毛片|
女人香蕉久久**毛片精品|
久久国产精品无码一区二区三区|
久久伊人精品一区二区三区
|
伊人伊成久久人综合网777|
欧美精品一区二区精品久久|
亚洲AV日韩AV天堂久久|
一本久久知道综合久久|
久久亚洲国产成人精品无码区|
国产AV影片久久久久久|
色综合久久中文色婷婷|
狠狠色丁香婷综合久久|
婷婷久久综合九色综合98|
国产V综合V亚洲欧美久久|
久久被窝电影亚洲爽爽爽|
成人国内精品久久久久影院|
狠狠色婷婷综合天天久久丁香
|
国产L精品国产亚洲区久久|
国内精品久久久久久久久电影网|
亚洲成色999久久网站|