| 查看: 181 | 回复: 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 ] |
回复此楼» 猜你喜欢
0856求调剂
已经有13人回复
-
0703化学
已经有11人回复
-
329求调剂
已经有7人回复
-
321求调剂
已经有8人回复
-
0856材料化工调剂 总分330
已经有12人回复
-
309求调剂
已经有3人回复
-
070300化学354求调剂
已经有4人回复
-
330分求调剂
已经有3人回复
-
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
zhougc
木虫 (正式写手)
黑社会老大
- 应助: 0 (幼儿园)
- 金币: 2590.9
- 散金: 120
- 帖子: 982
- 在线: 43.8小时
- 虫号: 393957
- 注册: 2007-06-06
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2008-10-09 16:13:16
oytxtu
铁杆木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 5912.3
- 散金: 20
- 红花: 1
- 帖子: 536
- 在线: 177.1小时
- 虫号: 484191
- 注册: 2007-12-26
- 专业: 凝聚态物性I:结构、力学和
3楼2008-10-09 16:14:06
