版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

hznu2007

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 33.2
散金: 22
帖子: 47
在线: 29.1小时
虫号: 658352
注册: 2008-11-20
专业: 凝聚态物性 II ：电子结构

[求助] FORTRAN程序DO WHILE出现语句错误，请哪位高手修改一下（紧急求助）已有1人参与

program main
  use IMSL
  implicit none
  external FCN

  real, parameter :: ERRREL=0.0001
  integer, parameter :: N=2
  integer, parameter :: ITMAX=100
  real :: XGUESS(N)=(/ 0.5, 0.5 /)
  real X(N) , FNORM

  CALL NEQNF(FCN, ERRREL, N, ITMAX, XGUESS, X, FNORM)
  write(*,*) X
  stop
end

subroutine FCN(XA, F, N)
implicit none
integer N
external FK, FE
double precision a, b, c, Romberg
real, target :: XA(N)
real F(N)
real, pointer :: x, ita
x=> XA(1)
ita=> XA(2)

c=Romberg(a, b, ita, FK)
F(1)=2*c/3.14159265+0.5/sqrt(1-ita**2)-0.5
c=Romberg(a, b, ita, FE)
F(2)=2*c/3.14159265+0.5*sqrt(1-ita**2)-2.5+x*(ita**2)

return
end subroutine

double precision function FE(g,ita)
  double precision g, ita
  FE=( 1-(ita*cos(g))*(ita*cos(g))) **0.5
end

double precision function FK(g,ita)
  double precision g, ita
  FK=( 1-(ita*cos(g))*(ita*cos(g))) **0.5
  FK=1/FK
end

double precision function Romberg(a, b, ita, f)
  double precision a, b, ita, f, eps
  integer m, i, k, point
  parameter(m=30, eps=1e-9)
  double precision t(m), er, h, tt, t2, l

  h=b-a
  k=1
  er=1
  point=1
  t(1)=h/2*(f(a,ita)+f(b,ita))

do 10 while( (er.gt.eps).and.(k.lt.n))
   tt=t(1)

   t(1)=0
   do 20 i=1,point
   t(1)=t(1)+f(a+i*h-h/2,ita)
   continue
   t(1)=tt/2+t(1)*h/2

   l=1
   do 30 i=2,k
      t2=t(i)
   l=l*4
   t(i)=(t(i-1)*l-tt)/(l-1)
   tt=t2
   continue

   point=point*2
   k=k+1
   h=h/2
   l=l*4
   t(k)=(t(k-1)*1-tt)/(l-1)
   er=abs(t(k)-tt)

   continue
   Romberg=t(k)
end do

end
--------------------Configuration: 201401 - Win32 Debug--------------------
Compiling Fortran...
D:\compaq visual fortran\MSDEV98\MyProjects\hello\01.f90
D:\compaq visual fortran\MSDEV98\MyProjects\hello\01.f90(70) : Error: This DO variable has already been used as an outer DO variable in the same nesting structure. [I]
do 30 i=2,k
-----------------^
D:\compaq visual fortran\MSDEV98\MyProjects\hello\01.f90(86) : Error: This ENDDO statement requires a label. [30]
end do
^
D:\compaq visual fortran\MSDEV98\MyProjects\hello\01.f90(70) : Error: There is an extra DO construct label. [30]
do 30 i=2,k
--------------^
D:\compaq visual fortran\MSDEV98\MyProjects\hello\01.f90(64) : Error: An unterminated block exists.
do 20 i=1,point
^
D:\compaq visual fortran\MSDEV98\MyProjects\hello\01.f90(60) : Error: An unterminated block exists.
do 10 while( (er.gt.eps).and.(k.lt.n))
^
D:\compaq visual fortran\MSDEV98\MyProjects\hello\01.f90(60) : Error: This label is undefined. [10]
do 10 while( (er.gt.eps).and.(k.lt.n))
-------^
D:\compaq visual fortran\MSDEV98\MyProjects\hello\01.f90(64) : Error: This label is undefined. [20]
do 20 i=1,point
--------------^
Error executing df.exe.

01.obj - 7 error(s), 0 warning(s)

回复此楼

» 猜你喜欢

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

fortran 权威指南英文版以及fortran的教学ppt 已经有235人回复
急，急，急，求一fortran语言高手，帮解决程序问题已经有8人回复
fortran输出求助已经有6人回复
急求fortran运行错误原因，在线等已经有7人回复
请教 fortran 运行错误的原因已经有13人回复
FORTRAN新手求助主程序循环问题已经有10人回复
请热心人能提供一个带imsl库的intel visual fortran编译器下载链接已经有6人回复
大家都用什么fortran编辑器，推荐一下，最好是能加亮或修改语句颜色的已经有15人回复
大家帮我看一个fortran的程序，我总是计算不出正确的结果已经有20人回复
请教一个fortran小程序编译出错的问题，谢谢已经有9人回复
Fortran的格式化输入输出问题已经有14人回复
经典书籍《数值算法大全（第3版）》：英文原版，含C、C++和Fortran源代码已经有303人回复
有没有人关注fortran语言的GPU并行呢？已经有12人回复
fortran两个问题，文件中的空格，屏幕同一行覆盖输出已经有23人回复
写了一个fortran90的小程序，编译通不过，请大侠帮忙已经有59人回复
【分享】Fortran语言简介与Fortran语句功能分析（PPT) 已经有147人回复
【求助】Fortran编译目标文件后的链接问题已经有4人回复

1楼 2014-01-12 22:00:46

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

jpchou

专家顾问 (著名写手)

专家经验: +20
应助: 231 (大学生)
金币: 2342
散金: 4367
红花: 103
帖子: 1061
在线: 757.3小时
虫号: 585915
注册: 2008-08-08
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

【答案】应助回帖

★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2014-01-14 01:49:13

D:\compaq visual fortran\MSDEV98\MyProjects\hello\01.f90(70) : Error: This DO variable has already been used as an outer DO variable in the same nesting structure. [I]
do 30 i=2,k
-----------------^
D:\compaq visual fortran\MSDEV98\MyProjects\hello\01.f90(86) : Error: This ENDDO statement requires a label. [30]
end do
^

你 do 30 i=2,k 給了個 label "30"
可是 end do 忘了給
把 end do 改成 30 end do
就可以
其他地方也要做同樣的修改

赞一下