±±¾©Ê¯ÓÍ»¯¹¤Ñ§Ôº2026ÄêÑо¿ÉúÕÐÉú½ÓÊÕµ÷¼Á¹«¸æ
²é¿´: 1175  |  »Ø¸´: 7
¡¾½±Àø¡¿ ±¾Ìû±»ÆÀ¼Û4´Î£¬×÷ÕßѧԱ7eU8hQÔö¼Ó½ð±Ò 3.5 ¸ö
µ±Ç°Ö÷ÌâÒѾ­´æµµ¡£
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

zsuzwr

ľ³æ (ÕýʽдÊÖ)


[×ÊÔ´] ¡¾·ÖÏí¡¿²úÉú¸ß˹·Ö²¼Ëæ»úÊýÔ´´úÂë(Fortran)

subroutine gauss(x1,x2,x,w,n)
      integer n
      real x1,x2,x(n),w(n)
      double precision xm,xl,z,p1,p2,p3,pp,z1
      integer m,i,j
      parameter (eps=3.d-14)
      m=(n+1)/2
      xm=0.5d0*(x2+x1)
      xl=0.5d0*(x2-x1)
      do 12 i=1,m
        z=cos(3.141592654d0*(i-.25d0)/(n+.5d0))
1       continue
          p1=1.d0
          p2=0.d0
          do 11 j=1,n
            p3=p2
            p2=p1
            p1=((2.d0*j-1.d0)*z*p2-(j-1.d0)*p3)/j
11        continue
          pp=n*(z*p1-p2)/(z*z-1.d0)
          z1=z
          z=z1-p1/pp
        if(abs(z-z1).gt.eps)go to 1
        if (abs(z).lt.eps) z=0.
        x(i)=xm-xl*z
        x(n+1-i)=xm+xl*z
        w(i)=2.d0*xl/((1.d0-z*z)*pp*pp)
        w(n+1-i)=w(i)
12    continue
      return
      end

[ Last edited by woshilsh on 2009-2-20 at 16:32 ]
»Ø¸´´ËÂ¥

» ÊÕ¼±¾ÌûµÄÌÔÌûר¼­ÍƼö

ÎÄÏײéÔÄ

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

lvzhu2007

ÈÙÓþ°æÖ÷ (ÖøÃûдÊÖ)


Ò»°ãÓþùÔÈ·Ö²¼µÄËæ»úÊýÀ´²úÉú¸ß˹·Ö²¼µÄ£¬ºÜ³£¼ûÒ²ºÜ¼òµ¥
http://www.taygeta.com/random/gaussian.html
4Â¥2009-02-21 09:22:51
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
²é¿´È«²¿ 8 ¸ö»Ø´ð

chenjinglzz

гæ (³õÈëÎÄ̳)



Ö§³ÖÏÂ
²»¹ýFORTRANÓïÑÔÎÒ²»ÊìϤ
Ï£ÍûÄÜ¿´Ã÷°×
2Â¥2008-05-05 17:35:26
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

scdxchenfang

ľ³æ (ÕýʽдÊÖ)


¡ï¡ï¡ï¡ï¡ï ÎåÐǼ¶,ÓÅÐãÍÆ¼ö

¡ï ¡ï ¡ï
woshilsh(½ð±Ò+3,VIP+0):¶àл½»Á÷£¬³£À´£¬»¥Ïà°ïÖú 2-20 16:32
Â¥Ö÷£¬ÎÒÔõôÔËÐÐÓкöà´íÎó°¡
--------------------Configuration: shuijishu - Win32 Debug--------------------
Compiling Fortran...
D:\exe\shuijishu.for
f90: Error: Illegal character in statement label field  [s]
f90: Error: Illegal character in statement label field  
f90: Error: Illegal character in statement label field  
f90: Error: Illegal character in statement label field  [r]
f90: Error: Illegal character in statement label field  [o]
f90: Error: First statement in file must not be continued
D:\exe\shuijishu.for(1) : Error: Syntax error, found END-OF-STATEMENT when expecting one of: ( % . = =>
subroutine gauss(x1,x2,x,w,n)
-----------------------------^
D:\exe\shuijishu.for(30) : Error: A RETURN statement is invalid in the main program.
      return
------^
D:\exe\shuijishu.for(3) : Error: An automatic object is invalid in a main program.   [X]
      real x1,x2,x(n),w(n)
-----------------^
D:\exe\shuijishu.for(3) : Error: An automatic object is invalid in a main program.   [W]
      real x1,x2,x(n),w(n)
----------------------^
D:\exe\shuijishu.for(2) : Error: A specification expression object must be a dummy argument, a COMMON block object, or an object accessible through host or use association   [N]
      integer n
--------------^
Error executing df.exe.

shuijishu.obj - 11 error(s), 0 warning(s)
ÎÒ¶Ôfortran²»ÊÇÌ«¶®£¬Äã¿ÉÒÔ¸øÎÒµÄÒâ¼ûÂð£¬Ð»Ð»
3Â¥2009-02-20 15:37:15
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

scdxchenfang

ľ³æ (ÕýʽдÊÖ)


Ŷ£¬Ð»Ð»Ä㣬ÎÒÏÈѧϰѧϰ
5Â¥2009-02-21 09:59:30
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
¡î ÎÞÐǼ¶ ¡ï Ò»ÐǼ¶ ¡ï¡ï¡ï ÈýÐǼ¶ ¡ï¡ï¡ï¡ï¡ï ÎåÐǼ¶
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] ¿¼Ñе÷¼Á +4 Sanmu-124 2026-03-26 4/200 2026-03-27 17:49 by kiokin
[¿¼ÑÐ] 328Çóµ÷¼Á +6 àŵεĻù±¾¶¼ 2026-03-27 6/300 2026-03-27 16:49 by sanrepian
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤085600£¬×Ü·Ö304£¬±¾¿ÆÓÐÁ½Æªsci²ÎÓ룬Çóµ÷¼Á +10 ÐÒÔ˵Ľ´½´ 2026-03-22 12/600 2026-03-27 16:08 by muchong357
[¿¼ÑÐ] 322Çóµ÷¼Á +4 ¾É… 2026-03-24 4/200 2026-03-27 15:38 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤µ÷¼Á +10 A-¶ßÀ²ZÃÎ 2026-03-23 16/800 2026-03-27 15:13 by caszguilin
[¿¼ÑÐ] 085600£¬²ÄÁÏÓ뻯¹¤321·Ö£¬Çóµ÷¼Á +9 ´ó²öС×Ó 2026-03-27 9/450 2026-03-27 14:30 by mmm just
[¿¼ÑÐ] Ò»Ö¾Ô¸211£¬335·Ö£¬0856£¬Çóµ÷¼ÁԺУºÍµ¼Ê¦ +4 Çã____Ïô 2026-03-27 5/250 2026-03-27 11:52 by zhshch
[¿¼ÑÐ] 0703»¯Ñ§/290Çóµ÷¼Á/±¾¿Æ¾­Àú·á¸»/¹¤¿ÆÒ²¿É +7 µ¤ÇàÄÌ¸Ç 2026-03-26 8/400 2026-03-27 11:24 by µ¤ÇàÄ̸Ç
[¿¼ÑÐ] 312Çóµ÷¼Á +9 Éϰ¶°ÉZJY 2026-03-22 13/650 2026-03-27 11:24 by sanrepian
[¿¼ÑÐ] 276Çóµ÷¼Á¡£ÓаëÄêµç³ØºÍ°ëÄê¸ß·Ö×Óʵϰ¾­Àú +10 ²ÄÁÏѧ257Çóµ÷¼Á 2026-03-23 11/550 2026-03-27 10:13 by YCIT- LHL
[¿¼ÑÐ] Ò»Ö¾Ô¸±±»¯Çóµ÷¼Á +3 Jsman 2026-03-22 3/150 2026-03-26 21:06 by ajpv·çÀ×
[¿¼ÑÐ] 340Çóµ÷¼Á +3 Amber00 2026-03-26 3/150 2026-03-26 18:57 by ²»³Ôô~µÄ؈
[¿¼ÑÐ] Ò»Ö¾Ô¸ÉϺ£½»´óÉúÎïÓëҽҩר˶324·Ö£¬Çóµ÷¼Á +6 jiajunX 2026-03-22 6/300 2026-03-25 23:05 by licg0208
[¿¼ÑÐ] 318Çóµ÷¼Á +3 plumÀî×Ó 2026-03-23 3/150 2026-03-25 09:42 by ÎíÉ¢ºóÏàÓölc
[¿¼ÑÐ] Ò»Ö¾Ô¸ÎäÀí085500»úеרҵ×Ü·Ö300Çóµ÷¼Á +3 an10101 2026-03-24 7/350 2026-03-25 00:00 by ɽ¹í0-
[¿¼ÑÐ] 292Çóµ÷¼Á +4 ¶ì¶ì¶ì¶î¶î¶î¶î¶ 2026-03-24 4/200 2026-03-24 16:41 by peike
[¿¼ÑÐ] 341Çóµ÷¼Á(Ò»Ö¾Ô¸ºþÄÏ´óѧ070300) +5 ·¬ÇÑÍ·--- 2026-03-22 6/300 2026-03-23 23:45 by Txy@872106
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖØÇì´óѧ085700×ÊÔ´Óë»·¾³£¬×Ü·Ö308Çóµ÷¼Á +7 īīĮ 2026-03-23 8/400 2026-03-23 20:36 by Creta
[¿¼ÑÐ] Çóµ÷¼ÁԺУÐÅÏ¢ +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by ÎÞи¿É»÷111
[¿¼ÑÐ] ÇóÖú +5 ÃÎÀïµÄÎÞÑÔ 2026-03-21 6/300 2026-03-21 17:51 by ѧԱ8dgXkO
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û