不明白也,能不能详细说一说啊
$ Q; w' H8 c( B: M这是我的代码:& y$ c2 ?7 h7 E* ]" Q- w3 F0 M
错在哪里????" c) Q: R( Q {. e0 b+ Q' x
void main()/ J- }2 n. b+ L( e5 F. ^
{n=4;' I+ l+ ~: ?0 O8 E$ J6 M7 U0 @
cout<<"请输入基向量 bb["<) J/ @& R5 |' l7 \( y5 G" xfor(i=1;i<=n-1;i++)
; ~/ \" d) i# g1 t{for(j=1;j<=n-1;j++) / a2 [3 W5 k9 b; i4 V1 w3 H
cin>>bb[j];( c! P: C6 ~# ~" I, D7 ?
cout<<'\n';}0 z m0 s5 L# d' h) ]+ D* h+ n
void computertb() //求逆矩阵tb[]
" p# i0 k& M0 J. \2 O{ for(i=1;i<=n-1;i++)
0 l ~# _5 C* y) J' T7 h/ I( G- \{for(j=n;j<=2*n-2;j++)
$ R2 p- D* J" z+ @{if (i==j-n+1) bb[j]=1;
" x/ v" i+ D: F2 U" S7 ^/ @else bb[j]=0;
2 C5 A2 \" B6 c}}* n* i6 G) F) Z. H1 |1 b
for(i=1;i<=n-1;i++)8 P4 W3 k- j0 D- l4 S" r; j; U
{if(bb!=1)# O$ f. b Y' Y0 x% ~! _% c
{for (j=i+1;j<=2*n-2;j++) bb[j]=bb[j]/bb;
) ]* ~" K* c. G( l( p/ \bb=1;}# x4 G5 ?; l7 k, V2 [6 x, k# }7 E
for(k=i+1;k" \/ D5 b4 W$ I! ]/ Z% q
{if(bb[k]!=0)& r R: w3 d. M
{for(j=i+1;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];1 a, K/ _5 \( ?; Z6 ]
bb[k]=0;
- ^( k) F; y8 T; L1 Z1 P3 A# p) G3 W}}}
$ d+ L6 T! o$ A$ A2 hfor(i=n-1;i>=1;i--)
% k" Q& o+ r6 |1 e{for(k=i-1;k>=1;k--)
- l5 e0 Y2 f1 d0 a& g( qif(bb[k]!=0)
4 ~1 Y4 E$ N# X{for(j=n;j<=2*n-2;j++) bb[k][j]=bb[k][j]-bb[k]*bb[j];
! F; @6 M$ J1 W5 Fbb[k]=0;
1 e, T3 A2 v" W# L" z, \: a% \6 `: ^}( U6 w4 A- f* c
}
, \! d; @) ^% y/ i7 `for(i=1;i2 }! o7 H" C5 g; V' qfor(j=n;j<=2*n-2;j++) tb[j-n+1]=bb[j];8 Y4 [' k K* Z- m2 G% G
for(i=1;i<=n-1;i++)
5 v* j! q5 r0 d1 L" F{for(j=1;j<= n-1;j++)
! g0 t* v4 r( t, d; _6 {cout<[j];7 W. G8 \ m2 ?
cout<<'\n';% X0 S9 J; s B0 b
}}' c5 _* Z4 W% r* k7 U# P
代码如上:
6 @ X( V ~; ]输入: 8 g1 S8 U- Q' C# |
1 2 -1) T/ a8 m9 a0 s3 j6 {
2 -2 1 {7 \5 y4 W" s5 O7 n
1 -1 2
" |. v5 t3 _! E; ?6 e S运行结果:& c, M% W) p% x# ^; o" \4 k
1 0 0 0.33 1 -1.3/ ^( D6 A3 O/ K& ^2 B$ h5 Y
0 1 0 0.33 -1 1.7
4 _: L! m; R k" r& S0 0 1 6e-008 -1 2
- n% u+ N0 m6 F/ A+ w5 l3 g) S+ r9 d" L2 |6 g7 d
# b* h4 {5 O1 M9 C, T1 B4 H! }
+ ]; Y4 G- h+ p- W# q5 T
|