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