24小时热门版块排行榜    

查看: 888  |  回复: 11
当前主题已经存档。

woshilsh

荣誉版主 (职业作家)

优秀版主

[交流] 【求助】一个简单的fortran程序也错误?(答谢500金币)【已解决】

用的版本compaq visual fortran v 6.6


程序如下:
program woshilsh
        parameter (N=10000,pi=3.1415927)
        real B(0:N),u1(0:N),u2(0:N)

        end       
错误就那一个,都是相同的,高手帮忙解答下了!

[ Last edited by nono2009 on 2009-11-27 at 09:20 ]
回复此楼
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了![/img][/url][/center]
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
woshilsh(金币+1,VIP+0):代码贴出来了,可老出错误,晕的不行了! 11-26 12:33
你这个错误怎么会发生在第19、23、35行?你帖出来到是ccc.f90文件的全部code吗?
2楼2009-11-26 12:13:31
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主

这个破工具,晕死我了
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了![/img][/url][/center]
3楼2009-11-26 12:30:33
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

Done

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
波不动(金币+3,VIP+0):找到问题,嘉奖! 11-26 13:06
before
program woshilsh
        parameter (N=10000,pi=3.1415927)
        real B(0:N),u1(0:N),u2(0:N)
        open(unit=1,file='u1.txt')
        open(unit=2,file='u2.txt')

after (program前要缩进或空格)
        program woshilsh
        parameter (N=10000,pi=3.1415927)
        real B(0:N),u1(0:N),u2(0:N)
        open(unit=1,file='u1.txt')
        open(unit=2,file='u2.txt')


[ Last edited by nono2009 on 2009-11-26 at 13:03 ]
4楼2009-11-26 13:00:47
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
nono2009(金币+1,VIP+0):谢谢。请转成jpg格式通过附件上传。 11-26 13:04
nono2009(金币+2,VIP+0):加奖! 11-26 13:04
我测试了,画了origin图片,不知道是否是这个结果??



       i=3*N/4 when t=dt*real(i)   :   u1   u2
      

       t=0 when i=0         :   u11   u22
   


      t=pi/4 when i=N/4 :   u111   u222

[ Last edited by maomao1210 on 2009-11-26 at 13:23 ]
5楼2009-11-26 13:01:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主

引用回帖:
Originally posted by maomao1210 at 2009-11-26 13:01:
我测试了,画了origin图片,不知道是否是这个结果??

[ Last edited by maomao1210 on 2009-11-26 at 13:04 ]

i=3*N/4             !t=0 when i=0; t=pi/4 when i=N/4; ......
利用这行注释,获得不同数据,一起画图,maomao再弄弄,

另外问下,这个在fortran里能画出来不?
还是一定要弄出数据在 orgin里画?谢谢
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了![/img][/url][/center]
6楼2009-11-26 13:08:43
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
jjdg(金币+1,VIP+0):画图用origin是个不错的选择! 11-26 13:27
woshilsh(金币+3,VIP+0):谢谢了,把你画的方法说下,发我邮箱吧,我给你发个邮件 11-26 15:04
引用回帖:
Originally posted by woshilsh at 2009-11-26 13:08:


i=3*N/4             !t=0 when i=0; t=pi/4 when i=N/4; ......
利用这行注释,获得不同数据,一起画图,maomao再弄弄,

另外问下,这个在fortran里能画出来不?
还是一定要弄出数据在 orgin里画?谢谢

fortran可以实现这个功能,不过很复杂,最好是计算出数据,然后再用专业画图软件进行画图,呵呵。
7楼2009-11-26 13:11:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主

引用回帖:
Originally posted by maomao1210 at 2009-11-26 13:01:
我测试了,画了origin图片,不知道是否是这个结果??

[ Last edited by maomao1210 on 2009-11-26 at 13:08 ]

第24行,i=0时对应t=0, i=N/4对应t=pi/4, i=N/2 对应t=pi/2, i=3N/4对应t=3pi/4, i=N 对应t=pi

分别设不同的i值,然后运行程序,就能得到每个x值下对应的u1,和u2,保存在文本文件内。这样就可以画图了

即第一次i=0 -> 运行 ->去生成的2个文本文件分别读取 u1,u2(这个得到的是纵坐标),至于横坐标,你可以用我给你发的EXCEL文件里面第一列的数据,横坐标是一直不变的

你每次改变i的值,只是得到不同的u1,u2

我这个是用最简单的2维画图,但已经很清楚反应问题的物理意义了。如果要生成三维图,即u1,u2为纵轴,x,t分别为2个横轴,程序需要改一点点。我建议找个稍微懂FORTRAN的人,当面一问就很清楚了


这是另外一个人的解答,后面的他让我继续问问别人,所以有劳maomao美女了哈!呵呵
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了![/img][/url][/center]
8楼2009-11-26 13:12:54
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

nono2009

超级版主 (文学泰斗)

No gains, no pains.

优秀区长优秀区长优秀区长优秀区长优秀版主

★
小木虫(金币+0.5):给个红包,谢谢回帖交流
还是把fortran翻译成matlab,用matlab画图可能效率更高些。
9楼2009-11-26 13:27:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
nono2009(金币+3,VIP+0):辛苦了! 11-26 14:17
见附件。。。。。。。。

[ Last edited by maomao1210 on 2009-11-26 at 14:05 ]
10楼2009-11-26 14:02:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 woshilsh 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗? 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[硕博家园] 中科院海洋所招收海洋材料相关领域博士研究生及研究助理 +3 sujiner 2026-10-08 4/200 2026-10-14 20:40 by qicaiyun
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见