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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » Fortran » 【交流】Fortran语言答疑专帖
1586/16返回列表首页上一页456789下一页末页
查看: 10722  |  回复: 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下一页末页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭