Fortran Coder

查看: 704|回复: 4
打印 上一主题 下一主题

[求助] 这个product的使用报错,不知道如何更改

[复制链接]

22

帖子

11

主题

0

精华

入门

F 币
97 元
贡献
54 点
跳转到指定楼层
楼主
发表于 2023-12-15 09:50:00 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式
[Fortran] 纯文本查看 复制代码
Subroutine Wii_()
Implicit none
real*8::x(11),y(11),pi=3.141592653589793,L
Integer::N,i,j
L=1
N=11
Do i=1,N,1
    x(i)=(L/2)*(1-cos(((i-1)/(N-1))*pi))
Do j=1,N,1
y(j)=(L/2)*(1-cos(((j-1)/(N-1))*pi))
  END Do
  END Do
end Subroutine Wii_
    
    program main_mik
    Implicit real*8(a-z)
    real*8::G(1),W(1),x(11),y(11)
    Integer::k
    k=1
call  Wii_()
 G(k)=product(x(i)-y(j),y(j)/=x(i))
 W(k)=1/G(k)
  write(*,*)"W(k)=W(k)"
  end program main_mik

微信截图_20231215094920.png (35.11 KB, 下载次数: 67)

这是问题

这是问题
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

835

帖子

2

主题

0

精华

大宗师

F 币
3926 元
贡献
2334 点
沙发
发表于 2023-12-15 11:20:35 | 只看该作者
product(x(i)-y,y/=x(i))

724

帖子

4

主题

0

精华

大师

农村外出务工人员

F 币
653 元
贡献
335 点

新人勋章爱心勋章水王勋章元老勋章热心勋章

板凳
发表于 2023-12-15 11:58:40 | 只看该作者
[Fortran] 纯文本查看 复制代码
Subroutine Wii_(x,y)
  Implicit none
  real*8::x(11),y(11),pi=acos(-1.0),L
  Integer::N,i,j
  L=1
  N=11
  Do i=1,N
    x(i)=(L/2)*(1-cos(((i-1.0)/(N-1.0))*pi))
  END Do
  Do j=1,N
    y(j)=(L/2)*(1-sin(((j-1.0)/(N-1.0))*pi))
  END Do
end Subroutine Wii_

program main_mik
  Implicit real*8(a-z)
  real*8::G(1),W(1),x(11),y(11)
  Integer::k
  k=1
  call  Wii_(x,y)
  G(k)=product(x-y,x/=y)
  W(k)=1/G(k)
  write(*,*)"W(k)=",W(k)
end program main_mik

22

帖子

11

主题

0

精华

入门

F 币
97 元
贡献
54 点
地板
 楼主| 发表于 2023-12-15 15:39:47 | 只看该作者
谢谢大佬,已经解决了

22

帖子

11

主题

0

精华

入门

F 币
97 元
贡献
54 点
5#
 楼主| 发表于 2023-12-15 15:40:28 | 只看该作者
li913 发表于 2023-12-15 11:20
product(x(i)-y,y/=x(i))

谢谢您的帮助,已解决
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-11-23 05:47

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表