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

24小时热门版块排行榜     石溪大学接受考研调剂申请>

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【交流】Fortran语言答疑专帖
【调剂】北京石油化工学院2024年16个专业接受调剂
1586/16返回列表首页上一页456789下一页末页
查看: 8798  |  回复: 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下一页末页
普通表情 龙 兔 虎 猫 高级回复(可上传附件) SCI英文润色翻译-EnPapers美团队-先修改后付款
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 前同事拿我以前未中的国自然标书申请今年的国自然，现在基金委发来调查函 +8 bjdxyxy 2024-04-19 10/500 2024-04-19 22:50 by 凌绝顶
[考研] 想问一下有没有需要文章但是没时间带学生的研究生导师。 +8 lekinna 2024-04-19 17/850 2024-04-19 22:06 by ld_cauchy
[考研] 298求调剂 +9 孙大大@ 2024-04-17 9/450 2024-04-19 21:50 by 刘国宁
[论文投稿] 一审一个审稿人，小修，会怎么样呀？ +9 林师妹 2024-04-18 9/450 2024-04-19 20:00 by nono2009
[论文投稿] 最近遇到这样一个问题 2+4 asd123gfa689 2024-04-18 10/500 2024-04-19 18:57 by asd123gfa689
[找工作] 事业单位还是大学好？ +17 青萍之沫 2024-04-16 18/900 2024-04-19 17:45 by charles-c
[找工作] 杭州国企和浙江高校如何选择？ +15 restart2024 2024-04-15 21/1050 2024-04-19 16:05 by l419110028
[基金申请] 院士建议：对勇于提出解决卡脖子问题新解决方案的学者加大资助力度 +7 zju2000 2024-04-16 10/500 2024-04-19 11:10 by dwuab
[有机交流] 紫外光谱 50+3 54胡 2024-04-17 3/150 2024-04-19 10:48 by Nanamiwww
[论文投稿] 求助 100+3 guhangyu 2024-04-15 3/150 2024-04-19 08:48 by bobvan
[博后之家] 博后换方向可行吗？ +3 越越不暴躁 2024-04-15 3/150 2024-04-18 10:58 by ciompman
[考研] 工学322，英语六级求调剂 +3 domax 2024-04-17 10/500 2024-04-17 21:25 by 东方花园小万
[考研] 322求调剂 +7 本己上岸 2024-04-16 7/350 2024-04-17 11:49 by duanxz
[考研] 332求调剂 +3 木叶下1999 2024-04-16 5/250 2024-04-17 00:56 by angeliar
[考研] 296求调剂 +3 Cclocomotive 2024-04-16 4/200 2024-04-16 10:04 by 19862091
[考研] 求调剂 +4 桃岸雪 2024-04-15 5/250 2024-04-15 18:49 by mthwyj
[论文投稿] with efitor 越久是不是越容易拒稿。我的已经一个多月了 +5 lizhengke06 2024-04-14 5/250 2024-04-15 18:33 by jonewore
[考研] 334求调剂 +4 学药救人 2024-04-14 4/200 2024-04-15 15:05 by hunanzang
[考研] 化学、材料类最后调剂机会！！！ +3 加油努力就好 2024-04-14 10/500 2024-04-15 09:05 by 任pen
[考研] 086000生医没考数学304分，一志愿西农求调剂 +3 王有田有学上 2024-04-13 5/250 2024-04-13 16:18 by coco1981
信息提示
关闭
请填处理意见
关闭