| 查看: 2127 | 回复: 1 | ||
hznu2007铁虫 (初入文坛)
|
[求助]
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) |
回复此楼» 猜你喜欢
本人考085602 化学工程 专硕
已经有15人回复
-
283求调剂
已经有3人回复
-
274求调剂
已经有5人回复
-
302求调剂
已经有6人回复
-
11408 一志愿西电，277分求调剂
已经有3人回复
-
化学工程321分求调剂
已经有13人回复
-
341求调剂
已经有5人回复
-
材料与化工304求B区调剂
已经有7人回复
-
0703化学336分求调剂
已经有4人回复
-
[复试调剂]西南科技大学国防/材料导师推荐
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
jjdg: 金币+1, 感谢参与 2014-01-14 01:49:13
感谢参与，应助指数 +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 就可以 其他地方也要做同樣的修改 |
2楼2014-01-13 20:23:47
50
