24小时热门版块排行榜

返回列表

zyj8119

木虫 (著名写手)

应助: 65 (初中生)
贵宾: 0.003
金币: 915.1
散金: 1440
红花: 35
帖子: 2936
在线: 1329.4小时
虫号: 664177
注册: 2008-11-29
性别: GG
专业: 理论和计算化学

[交流] 【求助】此程序怎么会数组越界？编译没有错误【已完结】已有2人参与

CODE:

PROGRAM main

REAL,DIMENSION(1:5,1:5)::a

REAL::det,determinate

a=reshape((/0.0,-3.0,-1.0,2.0,3.0,1.0,0.0,3.0,7.0,-2.0,3.0,&

-1.0,0.0,5.0,-5.0,2.0,6.0,-4.0,1.0,2.0,-2.0,5.0,0.0,&

-1.0,3.0/),(/5,5/))

det=determinate(a,5)

WRITE(*,'("Value of determinate A=",F8.2)')det

END PROGRAM main

!

FUNCTION determinate(a,n) RESULT(det)

INTEGER::n,rn,i,j,k,m

REAL::a(1:n,1:n),det,eps=1E-8

out:DO i=1,n-1

       rn=i

           IF(ABS(a(i,i))
in:        DO

           rn=rn+1

                   IF(rn>n)THEN

                   det=0.0;EXIT

                   ELSEIF(ABS(a(rn,i))>eps)THEN

                    DO m=i,n

           a(i,m)=a(i,m)+a(m,m)

                        END DO

                        EXIT in

                END IF

        END DO in

END IF

DO k=i+1,n

 DO j=i+1,n

  a(k,j)=a(k,j)-a(i,j)*a(k,i)/a(i,i)

 END DO

END DO 

END DO out

det=1.0

DO i=1,n

det=det*a(i,i)

END DO

END FUNCTION determinate

[ Last edited by 余泽成 on 2010-9-13 at 11:05 ]

回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

很多人会错失爱情，为什么？看完这个文章才觉得需要学习的不仅仅是数理化已经有15人回复
【纯水】值此节日之际想对公虫子说：已经有8人回复
【整理】《提问的智慧》——献给那些不会提问的人已经有78人回复

好好学习，天天向上。

1楼 2010-09-01 00:51:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

有目标有希望

★ ★ ★ ★
zyj8119(金币+1):谢谢参与
zyj8119(金币+2):谢谢回帖。 2010-09-01 10:27:33
nono2009(金币-5):专业区请勿纯表。谢谢！ 2010-09-01 13:08:29

2楼2010-09-01 07:51:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

sxf2012

木虫 (正式写手)

应助: 23 (小学生)
金币: 1133.8
红花: 3
帖子: 313
在线: 148.9小时
虫号: 900812
注册: 2009-11-12
专业: 光学

★
zyj8119(金币+1):谢谢参与
zyj8119(金币+2):是的。 2010-09-01 16:02:41

引用回帖:

Originally posted by zyj8119 at 2010-09-01 00:51:00:
[code]PROGRAM main
REAL,DIMENSION(1:5,1:5)::a
REAL::det,determinate
a=reshape((/0.0,-3.0,-1.0,2.0,3.0,1.0,0.0,3.0,7.0,-2.0,3.0,&
-1.0,0.0,5.0,-5.0,2.0,6.0,-4.0,1.0,2.0,-2.0,5.0,0.0,&
- ...

矩阵a是奇异的么？