Simpson's Rule
![]()
|
程式碼(最少47
Steps)
|
|||||||||
| 1. | f(K1) | 2. | KIN3 | 3. | K+4 | 4. | 1 | 5. | K+5 |
| 6. | 2 | 7. | + | 8. | 1 | 9. | +/- | 10. | XY |
| 11. | KOUT5 | 12. | - | 13. | KOUT5 | 14. | X2 | 15. | 1/X |
| 16. | FIX0 | 17. | RND | 18. | = | 19. | × | 20. | KOUT3 |
| 21. | = | 22. | K+4 | 23. | KOUT2 | 24. | - | 25. | KOUT1 |
| 26. | = | 27. | ÷ | 28. | MR | 29. | = | 30. | K+1 |
| 31. | K+2 | 32. | X |
33. | KOUT5 | 34. | X<M | 35. | MR |
| 36. | X |
37. | K-4 | 38. | 3 | 39. | K÷4 | 40. | 2 |
| 41. | K÷3 | 42. | X |
43. | K×4 | 44. | KOUT3 | 45. | X! |
| 46. | KOUT4 | 47. | NRM | ||||||
f(K1) 的意思是以K1作為變數 x 的值,輸入f(x)的程式碼。
執行程式前須先按 Shift AC ,把所有K寄存器set zero,然後把 a,b 分別存入K1和K2中,而 subinterval數則存入M中。
範例
分別用Simpson's 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. | 1 |
| 11. | K+5 | 12. | 2 | 13. | + | 14. | 1 | 15. | +/- |
| 16. | XY | 17. | KOUT5 | 18. | - | 19. | KOUT5 | 20. | X2 |
| 21. | 1/X | 22. | FIX0 | 23. | RND | 24. | = | 25. | × |
| 26. | KOUT3 | 27. | = | 28. | K+4 | 29. | KOUT2 | 30. | - |
| 31. | KOUT1 | 32. | = | 33. | ÷ | 34. | MR | 35. | = |
| 36. | K+1 | 37. | K+2 | 38. | X |
39. | KOUT5 | 40. | X<M |
| 41. | MR | 42. | X |
43. | K-4 | 44. | 3 | 45. | K÷4 |
| 46. | 2 | 47. | K÷3 | 48. | X |
49. | K×4 | 50. | KOUT3 |
| 51. | X! | 52. | KOUT4 | 53. | NRM |
對其他函數 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 |
得到下列結果:
25
(x2 + x + 4) dx = 61.5