Simpson's Rule and Trapezoidal Rule

程式碼(最少56 Steps)
1. f(K1) 2. KIN3 3. K+4 4. K+4 5. ×
6. 1 7. K+5 8. KOUT5 9. X2 10. 1/X
11. FIX0 12. RND 13. = 14. K-4 15. 1
16. + 17. 1 18. +/- 19. XY 20. KOUT5
21. = 22. × 23. KOUT3 24. = 25. K+6
26. KOUT2 27. - 28. KOUT1 29. = 30. ÷
31. MR 32. = 33. K+1 34. K+2 35. X<>Y
36. KOUT5 37. X<M 38. NRM 39. MR 40. X<>K3
41. K-4 42. X<>Y 43. K×4 44. K×6 45. KOUT4
46. K+6 47. ÷ 48. 2 49. K÷3 50. =
51. HLT 52. KOUT3 53. X! 54. 3 55. K÷6
56. KOUT6                

f(K1) 的意思是以K1作為變數 x 的值,輸入f(x)的程式碼。

執行程式前須先按 Shift AC ,把所有K寄存器set zero,然後把 ab 分別存入K1K2中,而 subinterval數則存入M

範例

分別用Simpson's Rule Trapezoidal Rulef(x) = x2 + x + 4,以 2 subintervals,計算25 f(x) dx。

f(K1) 的程式碼如下:

1. Kout1 2. X2 3. + 4. Kout1 5. +
6. 4 7. =            

整體的程式碼:

1. Kout1 2. X2 3. + 4. Kout1 5. +
6. 4 7. = 8. KIN3 9. K+4 10. K+4
11. × 12. 1 13. K+5 14. KOUT5 15. X2
16. 1/X 17. FIX0 18. RND 19. = 20. K-4
21. 1 22. + 23. 1 24. +/- 25. XY
26. KOUT5 27. = 28. × 29. KOUT3 30. =
31. K+6 32. KOUT2 33. - 34. KOUT1 35. =
36. ÷ 37. MR 38. = 39. K+1 40. K+2
41. X<>Y 42. KOUT5 43. X<M 44. NRM 45. MR
46. X<>K3 47. K-4 48. X<>Y 49. K×4 50. K×6
51. KOUT4 52. K+6 53. ÷ 54. 2 55. K÷3
56. = 57. HLT 58. KOUT3 59. X! 60. 3
61. K÷6 62. KOUT6            

對其他函數 f(x),可在 EDIT 狀態下,用SHIFT step1 step7 刪除,再輸入新函數 f(x)的程式碼,這樣就可把這程式應用到任何函數中。這亦是 EDIT 狀態的一大好處。

把所有K寄存器set zero

Shift AC

a=2b=5 分別存入K1K2中,subinterval=2 存入M

2 Kin 1 5 Kin 2 2 SHIFT MR

執行程式:

P1 RUN

得到下列結果:

25 (x2 + x + 4) dx
Trapezoidal Rule
62.625
Simpson's Rul
61.5

注意:subinterval數必須要是雙數,否則只能得Trapezoidal Rule的結果,當按RUN 便 E