[讨论]高斯消元法的程序对比

god

下面我收集的高斯消元法的程序。还是请大家看看，并提点意见。
2 L: \9 _, K  w. j& T  W2 ^  |

弘道

谢谢楼主……辛苦啦！………………

0.9清1.8清2.7清

wuyuwuyul 自己程序拿出来看看呗！干嘛还要体力A！

Rekcahpu

刚学高斯消元

deven1985

，不错，前几天的数值分析课上也是做的这个！

xiaozhangwh

谢谢分享，来试一下

hduscimc

#include<iostream>
* }) ^- c& T: Y6 }+ `#include<cmath>
; x- S7 o: ]* L( kusing namespace std;
# e  r  Z: X! Z+ g4 C1 Nint main()
/ E2 k- S9 s0 Y. L, a{
        int i;
        float *x;
        float c[3][4] ={0.101,2.304,3.555,1.183,
                                        -1.347,3.712,4.623,2.137,
                                        -2.835,1.072,5.643,3.035};
  S" R0 a4 B5 c7 X        float *ColPivot(float *,int);
6 K: G7 q# |8 b5 ?        x=ColPivot(c[0],3);
5 z4 G1 [( {3 [, e3 n        for( i=0;i<=2;i++ )
                cout<<"x("<<i<<")="<<x[i]<<endl;
! ^# k" n) R4 g4 B" O3 `* |}

float *ColPivot( float *c,int n )
{
. \. N7 h. ^8 X8 O        int i,j,t,k;
$ Y# I$ Y3 U4 b' \  a        float *x,p;
        x=new float[n*sizeof(float)];
5 _/ G7 s7 j8 n+ D. A8 _' k* b: J) ~        for( i=0;i<=n-2;i++)
( @3 `$ G4 E9 t* F! M        {
                k=i;
                for(j=i+1;j<=n-1;j++)
0 S# F5 Z4 N8 w/ e0 N) B5 P% }                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
                                k=j;
                        if(k!=i)
& a. i( y  |( k3 x                                for( j=i;j<=n;j++ )
  O' z8 T7 t6 i5 S; n9 `8 W                                {
                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
% A1 X& t; R( K9 R' J                                        *(c+k*(n+1)+j)=p;
                                }
                for( j=i+1;j<=n-1;j++ )
8 }3 J# n4 n; J5 A                {
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
                        for( t=i;t<=n;t++ )
                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
                }
        }
        for( i=n-1;i>=0;i--)
        {
                for( j=n-1;j>=i+1;j--)
                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
4 h/ v  u- G3 U                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
        }
        return x;
) {/ G' g2 X, {5 ~) R/ v( o}

shuainaruto

感谢分享，下来研究一下

水木年华zzu

总的来说写的不错!继续努力

LH611

谢谢分享，要又AS写的就更好了！

fansmiling

8错8错