±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 385·Ö ÉúÎïѧ£¨071000£©Çóµ÷¼Á +12 qf626 2026-04-01 12/600 2026-04-06 17:36 by 52305043001
[¿¼ÑÐ] 324Çóµ÷¼Á +3 k¿ÉÀÖ 2026-04-05 4/200 2026-04-06 09:54 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 377Çóµ÷¼Á +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[¿¼ÑÐ] 304Çóµ÷¼Á +3 luoye0105 2026-04-05 3/150 2026-04-05 18:16 by ÍÁľ˶ʿÕÐÉú
[¿¼ÑÐ] 296Çóµ÷¼Á +3 Íô£¡£¿£¡ 2026-04-05 5/250 2026-04-05 17:38 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 271·ÖÇóµ÷¼ÁѧУ +12 zph158488£¡ 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[¿¼ÑÐ] »·¾³285·Ö£¬¹ýÁù¼¶£¬Çóµ÷¼Á +10 xhr12 2026-04-02 10/500 2026-04-04 21:53 by bn53987
[¿¼ÑÐ] 081200-11408-276ѧ˶Çóµ÷¼Á +5 ´Þwj 2026-03-31 5/250 2026-04-04 19:45 by 1753564080
[¿¼ÑÐ] Ò»Ö¾Ô¸0817»¯Ñ§¹¤³ÌÓë¼¼Êõ£¬Çóµ÷¼Á +24 ÎÒ²»ÊÇÖ»Òò 2026-04-02 28/1400 2026-04-04 15:15 by dongzh2009
[¿¼ÑÐ] 26µ÷¼Á 086003 +6 ʧ»îµÄϸ°û 2026-04-04 6/300 2026-04-04 09:50 by zhangdingwa
[¿¼ÑÐ] 305Çóµ÷¼Á +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
[¿¼ÑÐ] 288Çóµ÷¼Á Ò»Ö¾Ô¸¹þ¹¤´ó ²ÄÁÏÓ뻯¹¤ +39 ÂåÉñ¸ç¸ç 2026-03-31 41/2050 2026-04-03 21:51 by qlm5820
[¿¼ÑÐ] 334Çóµ÷¼Á +9 Trying] 2026-03-31 9/450 2026-04-03 15:18 by ×ÁçíØ¼
[¿¼ÑÐ] 303Çóµ÷¼Á +3 һɫÇåÓð 2026-04-02 4/200 2026-04-03 10:22 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 330Çóµ÷¼Á +3 °×ÉñÎØºôºô 2026-04-02 3/150 2026-04-03 10:15 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 325·Ö»¯Ñ§µ÷¼Á +5 15771691647 2026-04-02 5/250 2026-04-03 09:58 by ChemPharm
[¿¼ÑÐ] 312 »¯¹¤»òÖÆÒ©µ÷¼Á +8 ССī123 2026-04-02 9/450 2026-04-03 09:12 by zhouxiaoyu
[¿¼ÑÐ] 366Çóµ÷¼ÁÒ»Ö¾Ô¸¶«±±´óѧ +8 ÔËÆøÀ´µÃÈôÓÐËÆÎ 2026-04-02 8/400 2026-04-02 21:39 by dongzh2009
[¿¼ÑÐ] 302Çóµ÷¼ÁÒ»Ö¾Ô¸±±º½070300£¬±¾¿ÆÖ£´ó»¯Ñ§ +8 Ê¥ÈÕ¶úÂüÌõ 2026-04-01 11/550 2026-04-02 07:40 by chemdavid
[¿¼ÑÐ] Ó¢Ò»ÊýÒ»×Ü·Ö334Çóµ÷¼Á +4 ³ÂÑôÀ¤ 2026-03-31 4/200 2026-03-31 14:22 by ¼Çʱ¾2026
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û