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 ]

回复此楼

» 猜你喜欢

垃圾破二本职称评审标准已经有18人回复
职称评审没过，求安慰已经有53人回复
毕业后当辅导员了，天天各种学生超烦已经有5人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
EST投稿状态问题已经有7人回复

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

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

好好学习，天天向上。

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

已阅回复此楼关注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是奇异的么？