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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【交流】Fortran语言答疑专帖
1586/16返回列表首页上一页456789下一页末页
查看: 11413  |  回复: 157
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 程序强帖 ，点击这里进行查看

lirenshun

铁虫 (初入文坛)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
jjdg(金币+1,VIP+0):欢迎新虫！ 11-28 12:09
你好，我在编程时错误显示D:\PROJECTS\hhh\hhh.f90(54) : Error: The attributes of this name conflict with those made accessible by a USE statement.   [P]
其中p是我定义的一个函数.谢谢
赞一下(7人) 回复此楼
51楼2009-11-17 21:12:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
引用回帖:
Originally posted by lirenshun at 2009-11-17 21:12:
你好，我在编程时错误显示D:\PROJECTS\hhh\hhh.f90(54) : Error: The attributes of this name conflict with those made accessible by a USE statement.   [P]
其中p是我定义的一个函数.谢谢

能不能把其中一部分贴出来看看？
赞一下 回复此楼
52楼2009-11-18 08:53:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lirenshun

铁虫 (初入文坛)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
程序如下，是计算辛普森积分的一个程序

module INTEGRAL
  implicit none
  real, parameter :: t=3000
contains
! 产生数列
   
  subroutine GenerateData(datas, width, f)
    real datas(, width
        real, external :: f
        real r
        integer i,n
        n = size(datas,1)
        width = t/(n-1)
        r = 0
        do i=1,n
          datas(i) = f(r)
          r = r+width
        end do
  end subroutine
   real function f(t)
     implicit none
     real,parameter::x1=1e-7
     real,parameter::d=1e-10
     real,parameter::pi=3.1415926
     integer::t
     f=(x1/(60*t*((4*pi*d*t)**0.5)))*exp((-x1)/(4*d*t))
    end function


real function Simpson_Integral(datas, width)
  IMPLICIT NONE
  real datas(, width
  real sum
  integer i,n
  n = size(datas,1)
  if ( mod(n,2)==0 ) then
    write(*,*) "要有奇数条数据"
        stop
  end if
  sum = datas(1) + datas(n) ! 先算出头尾的和
  do i=2,n-1
    if ( mod(i,2)==0 ) then
      sum = sum + 4*datas(i) ! 把4*f(x)的部分累加起来
        else
      sum = sum + 2*datas(i) ! 把2*f(x)的部分累加起来
        end if
  end do
  Simpson_Integral = sum * width/3.0        ! SUM再乘上H/3 就好了
  return
end function
end module
! SIMPSON法积分范例
program main
use integral
implicit none
integer,parameter::n=1001
real,external::f

real datas(n),width
call generatedata(datas,width,f)
write(*,*)simpson_integral(datas,width)
stop
end
赞一下(1人) 回复此楼
53楼2009-11-18 09:17:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★
余泽成(金币+5,VIP+0):专家辛苦了！ 11-18 10:48
senlia(金币+3,VIP+0):辛苦了！专家同志，能不能麻烦您顺带谢谢是第几行有错误[注释一下就行]，方便提问者[毕竟人家也是有疑惑才问嘛，您的思路大概比人家清楚些] 11-18 12:18
senlia(金币+0,VIP+0):顺带写写[不好意思 上面打错两个字] 11-18 12:18
module INTEGRAL
  implicit none
  real, parameter :: t=3000
contains
! 产生数列
   
  subroutine GenerateData(datas, width, f)
    real datas( :   ), width
        real, external :: f
!_ 原来你声明的是integer      
       real r
        integer i,n
        n = size(datas,1)
        width = t/(n-1)
        r = 0
        do i=1,n
               
          datas(i) = f(r)
          r = r+width
       
        end do
               
  end subroutine

   real function f(t)
     implicit none
     real,parameter::x1=1e-7
     real,parameter::d=1e-10
     real,parameter::pi=3.1415926
     real::t
！——请确认你的函数形式正确性
     f=(x1/(60*t*((4*pi*d*t)**0.5)))*exp((-x1)/(4*d*t))
  
        end function


real function Simpson_Integral(datas, width)
  IMPLICIT NONE
  real datas(:    ), width
  real sum
  integer i,n
  n = size(datas,1)
  if ( mod(n,2)==0 ) then
    write(*,*) "要有奇数条数据"
        stop
  end if

  sum = datas(1) + datas(n) ! 先算出头尾的和
  do i=2,n-1
    if ( mod(i,2)==0 ) then
      sum = sum + 4*datas(i) ! 把4*f(x)的部分累加起来
        else
      sum = sum + 2*datas(i) ! 把2*f(x)的部分累加起来
        end if
  end do

  Simpson_Integral = sum * width/3.0        ! SUM再乘上H/3 就好了
  return
end function
end module


! SIMPSON法积分范例
program main
use integral
implicit none
integer,parameter::n=1001
!——这里我去掉了下面的话
!Li real,external::ff

real datas(n),width


call generatedata(datas,width,f)

write(*,*)simpson_integral(datas,width)
stop
end


有些自变量类型不对，还有你的函数也有问题，对比旧的看一下便知道了。。

[ Last edited by maomao1210 on 2009-11-18 at 13:22 ]
赞一下(22人) 回复此楼
54楼2009-11-18 10:45:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liumin2532

金虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
你好，我以前本科时学过Fortran77，用的是Microsoft Fortranpower station 4.0 现在用的是compaq visual Fortran 6.5，忘记了以前是怎么输入并执行的了，您能不能用图示或者详细说明一下啊，谢谢！
赞一下(1人) 回复此楼
55楼2009-11-24 10:51:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
nono2009(金币+0,VIP+0):建议通过PM或“引用回复该帖”，以便别人及时收到你的message. 11-24 12:10
wangen994(金币+50,VIP+0):代woshilshi奖励500金币 11-27 09:29
wangen994(金币+50,VIP+0):代woshilshi奖励500金币 11-27 09:29
wangen994(金币+50,VIP+0):代woshilshi奖励500金币 11-27 09:29
wangen994(金币+50,VIP+0):代woshilshi奖励500金币 11-27 09:29
wangen994(金币+50,VIP+0):代woshilshi奖励500金币 11-27 09:30
wangen994(金币+50,VIP+0):代woshilshi奖励500金币 11-27 09:30
wangen994(金币+50,VIP+0):代woshilshi奖励500金币 11-27 09:30
wangen994(金币+50,VIP+0):代woshilshi奖励500金币 11-27 09:31
wangen994(金币+50,VIP+0):代woshilshi奖励500金币 11-27 09:31
wangen994(金币+50,VIP+0):代woshilshi奖励500金币 11-27 09:31
手头上没有入门的fortran书籍吗？呵呵，书上都应该有的吧
赞一下(110人) 回复此楼
56楼2009-11-24 12:03:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

波不动

木虫 (正式写手)

Wave No Move
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+0,VIP+0):建议通过PM或“引用回复该帖”，以便别人及时收到你的message. 11-24 12:10
nono2009(金币+1,VIP+0):辛苦了！ 11-24 12:11
对啊。没有的话，可以下个电子版的书看看，平时查起来也方便。这个链接里的书就不错。
http://muchong.com/bbs/viewthread.php?tid=1366167&fpage=1
赞一下(21人) 回复此楼
端好自己的碗，吃好自己的饭。
57楼2009-11-24 12:06:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshilsh

荣誉版主 (职业作家)

优秀版主
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
小木虫(金币+0.5):给个红包，谢谢回帖交流
wangen994(金币-49,VIP+0):程序软件版为你减免50% 11-27 09:24
wangen994(金币-49,VIP+0):程序软件版为你减免50% 11-27 09:25
wangen994(金币-49,VIP+0):程序软件版为你减免50% 11-27 09:25
wangen994(金币-49,VIP+0):程序软件版为你减免50% 11-27 09:25
wangen994(金币-49,VIP+0):程序软件版为你减免50% 11-27 09:25
wangen994(金币-5,VIP+0):程序软件版为你减免50% 11-27 09:25
maomao帮助我进行了fortran解答，请从我账号扣除500金币给她，谢谢了！哈哈！
赞一下(61人) 回复此楼
[center][url=http://www.91cool.net/][img]http://id.91cool.net/sign/?name=小木虫印&say=各位版主辛苦了！[/img][/url][/center]
58楼2009-11-27 08:33:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangen994

荣誉版主 (著名写手)

将军
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
woshilsh(金币+2,VIP+0):3Q,好好忙你的，复习要紧！ 11-27 09:32
引用回帖:
Originally posted by woshilsh at 2009/11/27 08:33:
maomao帮助我进行了fortran解答，请从我账号扣除500金币给她，谢谢了！哈哈！

这个我来扣除吧，呵呵
赞一下(21人) 回复此楼
59楼2009-11-27 09:17:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

maomao1210

金虫 (正式写手)
nono2009(金币+0,VIP+0):快要领工资，这个要算上。呵呵 11-27 11:40
引用回帖:
Originally posted by wangen994 at 2009-11-27 09:17:

这个我来扣除吧，呵呵

惭愧，惭愧。。。。。
回复此楼
60楼2009-11-27 11:01:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 maomao1210 的主题更新
1586/16返回列表首页上一页456789下一页末页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本人考085602 化学工程 专硕 +11 不知道叫什么！ 2026-03-15 13/650 2026-03-16 15:10 by 了了了了。。
[考研] 材料与化工求调剂 +3 为学666 2026-03-16 3/150 2026-03-16 15:09 by 加号+
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考博] 欢迎申博同学联系 +3 天道酬勤2026686 2026-03-10 7/350 2026-03-15 19:03 by 天道酬勤2026686
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 328求调剂 +3 5201314Lsy！ 2026-03-13 6/300 2026-03-14 15:31 by hyswxzs
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 一志愿天津大学，英一数二305分求调剂，四六级已过 +8 小小番的茄 2026-03-09 8/400 2026-03-14 01:53 by JourneyLucky
[考研] 云南财经大学信息学院计算机学硕专硕学位点 +3 zjptai 2026-03-10 5/250 2026-03-14 01:23 by 飞行琦
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] 307求调剂 +7 超级伊昂大王 2026-03-10 7/350 2026-03-14 00:49 by JourneyLucky
[考研] 271求调剂 +10 生如夏花… 2026-03-11 10/500 2026-03-14 00:35 by 卖报员小雨
[考研] 材料工程,326分，求调剂 +6 KRSLSR 2026-03-10 6/300 2026-03-13 23:47 by JourneyLucky
[考研] 0805，333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 285 求调剂 资源与环境 一志愿北京化工大学 +3 未名考生 2026-03-10 3/150 2026-03-13 23:04 by JourneyLucky
[考研] 304求调剂（085602一志愿985） +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[基金申请] 提交后的基金本子，已让学校撤回了，可否换口子提交 +3 dut_pfx 2026-03-10 3/150 2026-03-11 08:38 by kudofaye
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
关闭
请填处理意见
关闭