24小时热门版块排行榜    

查看: 11413  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0703化学调剂,求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 8/400 2026-03-16 13:42 by Dream007008
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 311求调剂 +5 牛乳糖的卡卡 2026-03-10 5/250 2026-03-14 00:05 by JourneyLucky
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 333求调剂 +3 球球古力 2026-03-11 3/150 2026-03-13 21:27 by JourneyLucky
[考研] 302求调剂 +6 负心者当诛 2026-03-11 6/300 2026-03-13 16:11 by JourneyLucky
[考研] 085600材料与化工 309分请求调剂 +7 dtdxzxx 2026-03-12 8/400 2026-03-13 14:43 by jxchenghu
[考研] 0703一志愿211 285分求调剂 +4 ly3471z 2026-03-13 4/200 2026-03-13 13:00 by JourneyLucky
[考研] 304求调剂(085602一志愿985) +12 化工人999 2026-03-09 12/600 2026-03-13 12:02 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 321求调剂(食品/专硕) +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
信息提示
请填处理意见