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，很多东西还不连贯。
请大师们对上面给的式子编个程序示范一下。

回复此楼

» 猜你喜欢

药学硕士，第一、第二作者已发表6 篇 SCI，药理方向及相关方向2026年/2027年博士申请已经有6人回复
一篇MDPI论文改变了学习工作和生活已经有5人回复
26年博士申请自荐-电催化已经有3人回复
中国地质大学（北京）博士招生补录，数理学院材料科学与工程专业和材料与化工专业已经有6人回复
收到国自然专家邀请后几年才会有本子送过来评已经有4人回复
考博已经有5人回复
26年申博自荐-计算机视觉已经有4人回复
药化及相关博士的申请已经有3人回复

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

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

宇天旅梦

铜虫 (小有名气)

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

引用回帖:

15楼: Originally posted by phy_yo at 2013-08-25 20:24:55
y(2)=0 是我打错了应该是y(1)=0 不过在这里不影响这里你即便不写默认的x(1),y(1)都是0 改用其他数据的的话要记得定义数组的首项
format的用法我就不细说了这里用不到你回头自己看书吧
dat文件可以用记事 ...

这一部分明白了，可是自己写了一个和这大同小异的程序，出错了。

部分程序如下：
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))

n=(KY*50-i*SQRT(KX*KX+KY*KY))/
   $(SQRT(KX*KX+KY*KY+KY*KY*50*50)-50*50/(6.582122*10**28))

   w=(KY*50+i*SQRT(1+KY*KY))/
   $(SQRT(KX*KX+KY*KY+KY*KY*50*50))

   z=(KY*50-i*SQRT(1+KY*KY))/
   $(SQRT(KX*KX+KY*KY+KY*KY*50*50))
f=((m-n)*(w-z)*EXP(-i*SQRT(1+KY*KY))/
   $((m-w)*(n-z)*e**(2*i*SQRT(KX*KX+KY*KY)-(m-z)(n-z))
其中关系是 m(KX,KY) n(KX,KY)    w(KX,KY)  ,  z(KX,KY)

f(m(KX,KY), n(KX,KY), w(KX,KY), z(KX,KY))

给出的错误是：
--------------------Configuration: Text1 - Win32 Debug--------------------
Compiling Fortran...
C:\Program Files\DevStudio\MyProjects\Text1.for
C:\Program Files\DevStudio\MyProjects\Text1.for(28) : Error: Syntax error, found '(' when expecting one of: * :: , <END-OF-STATEMENT> ; : ) (/ [ + - . ] /) ** / > // .LT. < ...
   $(((m-w)*(n-z)*exp(2*i*SQRT(KX*KX+KY*KY))-(m-z)(n-z))
----------------------------------------------------^
C:\Program Files\DevStudio\MyProjects\Text1.for(28) : Error: Syntax error, found END-OF-STATEMENT when expecting one of: , )
   $(((m-w)*(n-z)*exp(2*i*SQRT(KX*KX+KY*KY))-(m-z)(n-z))
----------------------------------------------------------^
C:\Program Files\DevStudio\MyProjects\Text1.for(15) : Error: The shapes of the array expressions do not conform. [M]
m=(KY*50+i*SQRT(KX*KX+KY*KY))/
--------^
C:\Program Files\DevStudio\MyProjects\Text1.for(18) : Error: The shapes of the array expressions do not conform. [N]
n=(KY*50-i*SQRT(KX*KX+KY*KY))/
--------^
C:\Program Files\DevStudio\MyProjects\Text1.for(21) : Error: The shapes of the array expressions do not conform. [W]
   w=(KY*50+i*SQRT(1+KY*KY))/
------^
C:\Program Files\DevStudio\MyProjects\Text1.for(24) : Error: The shapes of the array expressions do not conform. [Z]
   z=(KY*50-i*SQRT(1+KY*KY))/
------^
Error executing df.exe.

Text1.exe - 6 error(s), 0 warning(s)

赞一下

回复此楼

16楼2013-08-25 21:45:24

已阅回复此楼关注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 是虚数