Máy tính tích phân Simpson là gì?
Công cụ này tính gần đúng giá trị tích phân xác định của một hàm số \(f(x)\) trên khoảng \([a, b]\) bằng quy tắc Simpson kép (composite Simpson). Đây là một phương pháp toán học thuần túy nên áp dụng được ở mọi nơi, không phụ thuộc vào quốc gia hay quy ước nào. Máy tính sẽ liên tục tinh chỉnh kết quả bằng cách tăng gấp đôi số khoảng chia nhỏ (2, 4, 8, 16, ... cho đến giá trị tối đa \(N\) mà bạn chọn), giúp bạn theo dõi quá trình giá trị hội tụ dần về kết quả chính xác.
Cách sử dụng
Nhập hàm số theo biến \(x\) (ví dụ 4/(1+x^2) hoặc sin(x)), thiết lập cận dưới \(a\) và cận trên \(b\), rồi chọn số khoảng chia tối đa \(N\). Bạn cũng có thể tùy chọn số chữ số có nghĩa cần hiển thị. Cả phần cận lẫn hàm số đều chấp nhận các hằng số như pi và e, cùng nhiều hàm toán học như sin, cos, tan, exp, ln, log10, sqrt, abs và nhiều hàm khác.
Giải thích công thức
Với số khoảng chia \(n\) là số chẵn trên đoạn \([a, b]\), bước nhảy được tính bằng \(h = (b - a) / n\). Tại các điểm nút \(x_i = a + i\cdot h\), quy tắc Simpson gán trọng số 1 cho hai điểm đầu mút, trọng số 4 cho các điểm trong có chỉ số lẻ, và trọng số 2 cho các điểm trong có chỉ số chẵn, rồi nhân tổng với \(h/3\). Phương pháp này có sai số bậc \(O(h^4)\) và cho kết quả chính xác tuyệt đối với mọi đa thức bậc ba trở xuống.
$$\int_{a}^{b} f(x)\,dx \approx \frac{h}{3}\left[ f(x_0) + 4\sum_{i\,\text{odd}} f(x_i) + 2\sum_{i\,\text{even}} f(x_i) + f(x_N) \right]$$ $$\text{where}\quad \left\{ \begin{aligned} h &= \dfrac{b - a}{N} \\ x_i &= a + i\,h \end{aligned} \right.$$
Ví dụ minh họa
Xét \(f(x) = 4/(1+x^2)\), với \(a = 0\), \(b = 1\), \(n = 4\). Khi đó \(h = 0{,}25\) và giá trị tại các nút lần lượt là \(4.000000\), \(3.764706\), \(3.200000\), \(2.560000\), \(2.000000\). Áp dụng quy tắc Simpson:
$$S = \frac{0{,}25}{3}\cdot\left[4 + 4\cdot(3.764706 + 2.560000) + 2\cdot 3.200000 + 2\right] = 3{,}141569.$$Khi tăng lên \(N = 64\), kết quả hội tụ về \(\pi \approx 3{,}14159265358979\).
Câu hỏi thường gặp
Tại sao n bắt buộc phải là số chẵn? Quy tắc Simpson ghép từng cặp khoảng liền kề để dựng một parabol đi qua ba điểm, nên tổng số khoảng phải là số chẵn. Các lũy thừa của 2 dùng trong máy tính này luôn là số chẵn.
Nếu b nhỏ hơn a thì sao? Kết quả khi đó chỉ đơn giản là số đối của tích phân trên \([b, a]\); công thức tự xử lý đúng trường hợp bước nhảy âm.
Còn các điểm kỳ dị (singularity) thì sao? Nếu \(f(x)\) không xác định tại một điểm nút (chia cho 0, lấy log của số không dương, hoặc căn bậc hai của số âm), kết quả sẽ không còn đáng tin cậy; lúc này máy tính sẽ báo lỗi thay vì đưa ra một con số sai lệch gây hiểu nhầm.
