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 ]

回复此楼

» 猜你喜欢

天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有11人回复
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了已经有9人回复
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有10人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有10人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有11人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有11人回复

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

很多人会错失爱情，为什么？看完这个文章才觉得需要学习的不仅仅是数理化已经有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是奇异的么？