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 |
| 36. | KOUT5 | 37. | X<M | 38. | NRM | 39. | MR | 40. | X |
| 41. | K-4 | 42. | X |
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,然後把 a,b 分別存入K1和K2中,而 subinterval數則存入M中。
範例
分別用Simpson's Rule 和 Trapezoidal
Rule,當 f(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 |
42. | KOUT5 | 43. | X<M | 44. | NRM | 45. | MR |
| 46. | X |
47. | K-4 | 48. | X |
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=2,b=5 分別存入K1和K2中,subinterval數=2 存入M中:
|
2 Kin 1 5 Kin 2 2 SHIFT MR |
執行程式:
|
P1 RUN |
得到下列結果:
|
Trapezoidal Rule
|
62.625
|
|
Simpson's Rul
|
61.5
|
注意:subinterval數必須要是雙數,否則只能得Trapezoidal Rule的結果,當按RUN 便 E。