²é¿´: 648  |  »Ø¸´: 1

lernen8

гæ (³õÈëÎÄ̳)

[ÇóÖú] Çó°ïÖú£¡¸ÄFunctionÖеĿɱäϵÊý½øÐÐÑ­»·¼ÆËãµÄÎÊÌâ

ûÓнð±ÒÁ˰¡£¡£¡£¡Ï£Íû·¢ÑïÀ׷澫Éñ°¡£¡Çó´óÉñ°ïæ½â¾öÎÊÌâ°¡£¡
ÎÒÏëÓÃFortran Çó½â·ÇÏßÐÔ·½³Ì×飬³ÌÐòÊÇMonte Carlo·¨
³ÌÐòÓÉÖ÷³ÌÐò£¬Function, SubroutineÈý²¿·Ö×é³É¡£
Function¸ø³öÁË·ÇÏßÐÔ·½³Ì×飬ͨ¹ýsubroutine½øÐмÆËã·ÇÏßÐÔ·½³Ì×飬Ö÷³ÌÐò¸ø³öÁ˳õʼ¼ÆËãÖµX1,X2,X3¡£

ÎÊÌâÊÇ£¬ Function ÖÐÓиöϵÊýMÊDZ仯µÄ£¬
ÔÚsubroutine£¨»òÕßÖ÷³ÌÐò£©ÖУ¬ÏÈÉ趨Á˳õʼֵM0£¬Ã¿´Î¼ÆËã½áÊøÁËXºó£¬ÓÉX¼ÆËãеÄM=f(x).
Èç¹ûÐÂMc ºÍ¾ÉM²î±ðÔÚ1.0E-05ÄھͽáÊø¼ÆË㣬
·ñÔòÖØÐ¸³ÓèFunction ϵÊýM, ÖØÐµ÷ÓÃ×Ó³ÌÐòÑ­»·¼ÆËãX.

     DIMENSION X(N), Y( N)
     REAL NRND1
     External F
      Double precision R, M, M0, Mc
      Double precision R=1.0E-4; d=30.0E-3; k=383.0;
     * h=40000.0; T=303.0
        A= B
        K= 1
        R= 1.0E0
        M0=450.0
40      M=M0
       Z= F(X, N, M)
  10   IF (A. GE. EPS) THEN
      L= L+ 1
      DO 20 I= 1, N
  20  Y(I) = -A+ 2.0* A* NRND1(R) + X(I)
      Z1= F(Y, N, M)
      K= K+ 1
      IF ( Z1. GE. Z) THEN
      IF (K .GT. M) THEN
       K= 1
       A= A/2.0
      ENDIF
      GOTO 10
      ELSE
       K= 1
      DO 30 I= 1, N
  30   X(I) = Y( I)
       Z= Z1
       IF ( Z. GE. EPS) GOTO 10
      ENDIF
       ENDIF
      
      Mc=T+(1000.0*X(1)-M)*(d/k+1/h)/R
      IF(ABS(M-Mc) .GE. EPS) THEN
      M0=(M+Mc)/2
      GOTO 40   
      
      ENDIF
      END

ÉÏÃæ³ÌÐòÊǶÔsubroutineµÄÐ޸ġ£µ¹ÊýµÚ6ÐÐÖ±½Óµ÷ÓÃX1¶ÔеÄMc¼ÆË㣬²¢¸³Öµ¸øM0£¬È»ºó·µ»ØÐÐ40¡£
µ«ÊÇÃ²ËÆ¸ù±¾¾ÍûÓаÑеÄM´«µÝ¸ø·½³ÌF(X,N,M)½øÐÐÑ­»·¼ÆËã¡£×îºóÖ»ÓеÚÒ»´ÎµÄ½á¹û¡£

ÄѵÀ¼ÆËãµÃµ½µÄX1²»ÄÜÖ±½Óµ÷Óã¬ÓÃÓÚ·½³Ì¼ÆËãô£¿(³¢ÊÔ¹ý°ÑX1Ö±½Óµ÷³öÖ÷³ÌÐò¼ÆËãµ¹ÊýµÚ6Ðеķ½³Ì£¬·½³Ì½á¹û¸ù±¾¾Í²»¶Ô)
´Ë×Ó³ÌÐòÓзñ´íÎóÄØ£¿
лл£¡
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:

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

ѧ³ædj

½ð³æ (ÕýʽдÊÖ)

¿´Á˰ëÌìû¿´¶®~Â¥Ö÷×îºÃ°Ñ´úÂë´ò°ü£¬ÎÊÌâ˵Çå³þ~
ÇîÔò¶ÀÉÆÆäÊ£¬´ïÔò¼æ¼ÃÌìÏÂ
2Â¥2014-05-26 16:40:35
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ lernen8 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 275Çóµ÷¼Á +3 Micky11223 2026-03-25 5/250 2026-03-25 22:32 by barlinike
[¿¼ÑÐ] Ò»Ö¾Ô¸ºÓ¹¤´ó 081700 276Çóµ÷¼Á +3 µØÇòÈÆ×ÅÌ«Ñôת 2026-03-23 3/150 2026-03-25 19:10 by ÎíÉ¢ºóÏàÓölc
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤ 322Çóµ÷¼Á +6 È»11 2026-03-19 6/300 2026-03-25 18:37 by haxia
[¿¼ÑÐ] ¿¼ÑÐÒ»Ö¾Ô¸ËÕÖÝ´óѧ³õʼ315£¨Ó¢Ò»£©Çóµ÷¼Á +3 sbdksD 2026-03-24 4/200 2026-03-25 18:16 by xcjcqu
[¿¼ÑÐ] ÍøÂç¿Õ¼ä°²È«0839Õе÷¼Á +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
[¿¼ÑÐ] 271Çóµ÷¼Á +4 ÉúÈçÏÄ»¨¡­ 2026-03-22 4/200 2026-03-25 11:25 by userper
[¿¼ÑÐ] 278Çóµ÷¼Á +4 ÎÒ¿ÉÒÔÉϰ¶µÄ¶Ô 2026-03-19 4/200 2026-03-25 11:01 by userper
[¿¼ÑÐ] ÉúÎïѧѧ˶Çóµ÷¼Á +7 СÑò˯×ÅÁË? 2026-03-23 10/500 2026-03-25 02:24 by Çå·ç·÷Ñï¡£ m
[¿¼ÑÐ] ²ÄÁÏѧÇóµ÷¼Á +6 Stella_Yao 2026-03-20 6/300 2026-03-25 00:37 by baoball
[¿¼ÑÐ] 070300»¯Ñ§Çóµ÷¼Á +9 Ô·¶¹¶¹ 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[¿¼ÑÐ] 307Çóµ÷¼Á +3 ÓàÒâÇä 2026-03-21 6/300 2026-03-24 15:03 by ÓàÒâÇä
[¿¼ÑÐ] 276Çóµ÷¼Á¡£ÓаëÄêµç³ØºÍ°ëÄê¸ß·Ö×Óʵϰ¾­Àú +9 ²ÄÁÏѧ257Çóµ÷¼Á 2026-03-23 10/500 2026-03-24 07:36 by wangy0907
[¿¼ÑÐ] 263Çóµ÷¼Á +6 yqdszhdap£­ 2026-03-22 9/450 2026-03-23 12:57 by yqdszhdap£­
[¿¼ÑÐ] 291Çóµ÷¼Á +5 ‹üÈA 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[¿¼ÑÐ] 293Çóµ÷¼Á +3 ÌÎÌÎWjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[¿¼ÑÐ] Çóµ÷¼Á +3 13341 2026-03-20 3/150 2026-03-21 18:28 by ѧԱ8dgXkO
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖØÇì´óѧ085700×ÊÔ´Óë»·¾³×Ü·Ö308Çóµ÷¼Á +7 īīĮ 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[¿¼ÑÐ] 330Çóµ÷¼Á0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by ²«»÷518
[¿¼ÑÐ] 085601µ÷¼Á 358·Ö +3 zzzzggh 2026-03-20 4/200 2026-03-21 10:21 by luoyongfeng
[¿¼ÑÐ] Çóµ÷¼ÁÒ»Ö¾Ô¸ÄϾ©º½¿Õº½Ìì´óѧ289·Ö +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û