| 查看: 189 | 回复: 2 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
oytxtu铁杆木虫 (正式写手)
|
[交流]
【求助】哪位程序高手帮我看看，这个程序（已获得助）
|
||
|
这是一个计算龙格-库塔的程序 老是出现： Linking... rk.obj : error LNK2001: unresolved external symbol _GRKT1@36 Debug/rk.exe : fatal error LNK1120: 1 unresolved externals Error executing link.exe. 这个错误 program rk dimension y(3),d(3),z(3,11),b(3) double precision y,d,z,t,h,b,m,n integer i,j t=0.0 y(1)=-1.0 y(2)=0.0 y(3)=1.0 h=0.01 m=3 n=11 call grkt1(t,y,m,h,n,z,f,d,b) write(*,*) do 10 i=1,n t=(i-1)*h write(*,50) t write(*,100) (z(j,i),j=1,m) write(*,*) 10 continue 50 format(1x,'t=',f7.3) 100 format(1x,'y(1)=',d13.6,3x,'y(2)=',d13.6,3x,'y(3)=',d13.6) end subroutine grekt1(t,y,m,h,n,z,f,d,b) dimension y(m),d(m),z(m,n),a(4),b(m) double precision y,d,z,a,b,t,h,x,tt a(1)=h/2.0 a(2)=a(1) a(3)=h a(4)=h do 5 i=1,m 5 z(i,1)=y(i) x=t do 100 j=2,n call f(t,y,m,d) do 10 i=1,m 10 b(i)=y(i) do 30 k=1,3 do 20 i=1,m y(i)=z(i,j-1)+a(k)*d(i) b(i)=b(i)+a(k+1)*d(i)/3.0 20 continue tt=t+a(k) call f(tt,y,m,d) 30 continue do 40 i=1,m 40 y(i)=b(i)+h*d(i)/6.0 do 50 i=1,m 50 z(i,j)=y(i) t=t+h 100 continue t=x return end subroutine subroutine f(t,y,m,d) dimension y(m),d(m) double precision y,d,t d(1)=y(2) d(2)=-y(1) d(3)=-y(3) return end subroutine [ Last edited by kuhailangyu on 2008-12-9 at 10:45 ] |
回复此楼» 猜你喜欢
河北省自然科学基金
已经有8人回复
-
西安交大新媒学院副院长用撤稿论文结题
已经有5人回复
-
论文撤稿了
已经有5人回复
-
某211大学教师把个人教师官方主页改成:我跑了我跑了我跑了!官宣跑路！
已经有5人回复
-
26/27申博自荐
已经有9人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有7人回复
-
揭秘青基评审内幕：几个A才能顺利中标
已经有4人回复
-
青B发送上会通知了吗
已经有7人回复
-
博士申请
已经有3人回复
oytxtu
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5912.3
- 散金: 20
- 红花: 1
- 帖子: 536
- 在线: 177.1小时
- 虫号: 484191
- 注册: 2007-12-26
- 专业: 凝聚态物性I:结构、力学和
3楼2008-10-09 16:14:06
zhougc
木虫 (正式写手)
黑社会老大
- 应助: 0 (幼儿园)
- 金币: 2590.9
- 散金: 120
- 帖子: 982
- 在线: 43.8小时
- 虫号: 393957
- 注册: 2007-06-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2008-10-09 16:13:16
10