版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

宇天旅梦

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 338.6
红花: 1
帖子: 134
在线: 57小时
虫号: 2178982
注册: 2012-12-10
性别: MM
专业: 凝聚态物性 II ：电子结构

[求助] 编一个fortran程序

需要一个fortran程序，比较简单的那种，例如f=((5x+6z)*e^(2i) )/(2x-z),其中z=(3y-1)*(x+2)

要求输出的结果是三列数据，分别对应x,y,f .并且每列有20个数。

这三组数据是要导入orgin中进行图像模拟的。
初学fortran，很多东西还不连贯。
请大师们对上面给的式子编个程序示范一下。

回复此楼

» 猜你喜欢

1楼 2013-08-24 18:46:41

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

phy_yo

铁虫 (小有名气)

应助: 3 (幼儿园)
金币: 452.3
散金: 1667
红花: 121
帖子: 255
在线: 60.8小时
虫号: 2414907
注册: 2013-04-12
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

16楼: Originally posted by 宇天旅梦 at 2013-08-25 21:45:24
这一部分明白了，可是自己写了一个和这大同小异的程序，出错了。

部分程序如下：
do 2 k=1,100
m=(KY*50+i*SQRT(KX*KX+KY*KY))
$/(SQRT(KX*KX+KY*KY+KY*KY*50*50)-50*50/(6.582122*10**28))

...

你这部分有点部分过头了吧循环变量k在循环体内根本没看到你应该不是故意删的……所以你的这些表达式不会是直接复制粘贴的吧……还全是大写

KX,KY是是自变量要做成数组形式如KX（20），KY（20）之类的再套入循环变量就行了你自己慢慢琢磨吧我感觉我前面写的那些代码够你做参考了~或者你可以考虑用matlab 作这种计算他比fortran方便

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

宇天旅梦

网络错误无法获取数据

17楼2013-08-25 22:20:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 18 个回答

phy_yo

铁虫 (小有名气)

应助: 3 (幼儿园)
金币: 452.3
散金: 1667
红花: 121
帖子: 255
在线: 60.8小时
虫号: 2414907
注册: 2013-04-12
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

感谢参与，应助指数 +1

open(10,file='data.dat',status='unknown')
real*8 x,y,z,x(20),y(20)
do 1 i=1,20
input (*,*) x,y
x(i)=x
y(i)=y
1 continue
do 2 1,20
z=(3y(i)-1)*(x(i)+2)
f=((5x+6z)*e^(2i) )/(2x-z)
write(*,*) x(i),y(i),f
write(10,*) x(i),y(i),f
2 continue
pause
end

赞一下

回复此楼

网络错误无法获取数据

2楼2013-08-24 20:27:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

宇天旅梦

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 338.6
红花: 1
帖子: 134
在线: 57小时
虫号: 2178982
注册: 2012-12-10
性别: MM
专业: 凝聚态物性 II ：电子结构

引用回帖:

2楼: Originally posted by phy_yo at 2013-08-24 20:27:45
open(10,file='data.dat',status='unknown')
real*8 x,y,z,x(20),y(20)
do 1 i=1,20
input (*,*) x,y
x(i)=x
y(i)=y
1 continue
do 2 1,20
z=(3y(i)-1)*(x(i)+2)
f=((5x+6z)*e^(2i) )/(2x-z)
write(*,* ...

式子中的 i 是虚数