24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 11518  |  回复: 157
本帖产生 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是我定义的一个函数.谢谢
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
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 ]
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,忘记了以前是怎么输入并执行的了,您能不能用图示或者详细说明一下啊,谢谢!
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书籍吗?呵呵,书上都应该有的吧
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
端好自己的碗,吃好自己的饭。
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金币给她,谢谢了!哈哈!
[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金币给她,谢谢了!哈哈!

这个我来扣除吧,呵呵
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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 302分求调剂 一志愿安徽大学085601 +11 zyx上岸! 2026-04-04 11/550 2026-04-06 22:55 by qlm5820
[考研] 求调剂 +20 111623 2026-04-04 22/1100 2026-04-06 22:52 by chenzhimin
[考研] 华南理工0703化学,总分336求调剂 +10 手机用户 2026-04-03 11/550 2026-04-06 22:20 by michael2011
[考研] 材料调剂 +12 一样YWY 2026-04-05 13/650 2026-04-06 15:38 by lin-da
[考研] 生物学求调剂 +5 15064154688 2026-04-03 5/250 2026-04-06 11:56 by lijunpoly
[考研] 化学调剂 +17 艾志恒 2026-04-03 18/900 2026-04-06 07:10 by jj987
[考研] 294求调剂 +4 Grey_Ey 2026-04-01 5/250 2026-04-05 23:05 by Grey_Ey
[考研] 315求调剂 +5 &123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 工科求调剂 +15 11ggg 2026-04-03 15/750 2026-04-05 16:24 by zzx2138
[考研] 291求调剂 +4 迷蒙木木 2026-04-01 5/250 2026-04-04 15:59 by sihailian3
[考研] 一志愿0817化学工程与技术,求调剂 +24 我不是只因 2026-04-02 28/1400 2026-04-04 15:15 by dongzh2009
[考研] 求材料调剂,一志愿郑州大学289分 +15 硕星赴 2026-04-03 15/750 2026-04-04 01:01 by userper
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +9 哇呼哼呼哼 2026-04-02 9/450 2026-04-03 12:05 by 1753564080
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 一志愿武汉理工0856,初试334 +3 26考研材料 2026-04-02 3/150 2026-04-02 21:22 by dongzh2009
[考研] 土木304求调剂 +4 兔突突突, 2026-04-02 5/250 2026-04-02 21:16 by 兔突突突,
[考研] 调剂 +3 好好读书。 2026-04-01 6/300 2026-04-02 15:49 by liumengping
[考研] 348环境工程调剂 +3 吴彦祖24k 2026-04-01 3/150 2026-04-02 09:14 by nanaliuyun
[考研] 350求调剂 +7 阿佳~ 2026-03-31 7/350 2026-04-01 16:12 by yanflower7133
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
信息提示
请填处理意见