24小时热门版块排行榜    

CyRhmU.jpeg
查看: 586  |  回复: 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 的主题更新
信息提示
请填处理意见