ÇÁ·Î±×·¡¹Ö
ù ÆäÀÌÁö ·£´ý ±Û ȸ¿ø°¡ÀÔ ·Î±×ÀÎ
ºñ°ø°³ ¼Õ´Ô ¡¦ 2008-06-16 00:34:14
URL https://te31.com/rgr/view.php?id=study&no=584 ¸ð¹ÙÀÏ È­¸é
ÇÁ·Î±×·¥ °í¼ö´Ôµé ºÁÁÖ¼¼¿ä..¤Ð

½É½¼ ¹ýÄ¢(simpson's rule)



f(x)=6+3cosx  ±¸°£Àº  a = 0 ,   b = ¥ð/2



@ °¡»óÄÚµå @

(a) ´ÜÀϱ¸°£ Simpson 1/3 °ø½Ä

       function simp13 (h, f0, f1, f2)

          simp13 = 2*h* (f0+4*f1+f2) / 6

       end simp13



(b) ´ÜÀϱ¸°£ Simpson 3/8 °ø½Ä

       function simp38 (h, f0, f1, f2, f3)

          simp38 = 3*h* (f0+3*(f1+f2)+f3) / 8

       end simp38



(c) º¹ÇÕ±¸°£ Simpson 1/3 °ø½Ä

      function simp13m ( h, n, f)

          sum = f(0)

          dofor i = 1, n-2, 2

              sum = sum + 4 * fi + 2 * fi+1

          end do

          sum = sum + 4 * fn-1 + fn

          simp13m = h * sum / 3

       end simp13m



(d) ±¸°£¼ö°¡ ¦¼ö ¹× Ȧ¼öÀÎ °æ¿ì¸¦ À§ÇÑ º¹ÇÕ±¸°£ Simpson 3/8°ø½Ä

       function simpint(a, b, n, f)

           h = (b - a) / n

           if n = 1 then

               sum = trap(h, fn-1, fn)

           else

               m = n

               odd = n / 2 - int(n / 2)

               if odd > 0 and n > 1 then

                   sum = sum+simp38(h, fn-3, fn-2, fn-1, fn)

                   m = n-3

               end if

               if m > 1 then

                   sum = sum + simp13m(h, m , f)

               end if

           end if

           simpint = sum

       end simpint



À§ÀÇ °¡»óÄڵ带 Åä´ë·Î Çؼ­ ÇÁ·Î±×·¥À» Â¥ÁÖ¼¼¿ä. º¹ÀâÇÑ°Å Àý´ë ÇÊ¿ä¾ø¾î¿ä.
±×³É ÀԷºθ¸ ³Ö¾î¼­ ÀûºÐ°ª¸¸ ¾Ë¼ö ÀÖ°Ô¸¸ Â¥ÁÖ¼¼¿ä. ÃÖ´ëÇÑ °£°á°£´ÜÇÏ°Ô ÇØÁÖ¼¼¿ä¤Ð

Àâ´ã | 2237¸íÀÌ Àоú¾î¿ä. 3.14.253.152

0
1 ºñ°ø°³ ¼Õ´Ô ¡¦ 2008-06-18 09:18:20
;;ÄÚ¸¦ ÆÄ´Þ¶ó°í ÇÏ¸é ¾î¶±ÇØ¿ä
´ñ±ÛÀ» ÀÛ¼ºÇÏ½Ç ¼ö ¾ø½À´Ï´Ù.
(±ÇÇÑÀÌ ¾ø´Â ȸ¿ø·¹º§)
¸ñ·ÏÀ¸·Î
¡â ÀÌÀü±Û: [¼ö´É] ÀÌ°úÀε¥ ¼ö¸®¿µ¿ª ¸ðÀÇ°í»ç ¹®Á¦Áý ÃßõÁ»¿ä [8]
¡ä ´ÙÀ½±Û: [¼ö´É] ºÎ»ê¿¡ Àç¼öÇпø ÁÁÀº°÷ ¾Æ´Â ÃÊ°í¼ö¸¸ [2]
ÀÌ¿ë¾à°ü | ±¤°í/Á¦ÈÞ | °³ÀÎÁ¤º¸Ãë±Þ¹æħ | ¹®ÀÇ/½Å°í | ¸ð¹ÙÀÏ TE31 | ¼­¹ö ºÎÇÏ : 15%
½Ç½Ã°£ Issue Ä¿¹Â´ÏƼ TE31 [¾ËÁö·Õ] ¨Ï 2002-2024
TOP arrow_upward