Trapezoidal Rule
![]()
|
程式碼(最少34
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. | KOUT2 |
| 16. | - | 17. | KOUT1 | 18. | = | 19. | ÷ | 20. | MR |
| 21. | = | 22. | K+1 | 23. | K+2 | 24. | X |
25. | KOUT5 |
| 26. | X<M | 27. | NRM | 28. | KOUT3 | 29. | K-4 | 30. | X |
| 31. | K×4 | 32. | 2 | 33. | K÷4 | 34. | KOUT4 | ||
f(K1) 的意思是以K1作為變數 x 的值,輸入f(x)的程式碼。
執行程式前須先按 Shift AC ,把所有K寄存器set zero,然後把 a,b 分別存入K1和K2中,而 subinterval數則存入M中。
範例
用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. | KOUT2 | 22. | - | 23. | KOUT1 | 24. | = | 25. | ÷ |
| 26. | MR | 27. | = | 28. | K+1 | 29. | K+2 | 30. | X |
| 31. | KOUT5 | 32. | X<M | 33. | NRM | 34. | KOUT3 | 35. | K-4 |
| 36. | X |
37. | K×4 | 38. | 2 | 39. | K÷4 | 40. | KOUT4 |
對其他函數 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 約等於 62.625。