Phương pháp chia đôi tìm nghiệp gần đúng của phương trình f(x) = 0 với khoảng cách ly nghiệm [a,b]
Code:
(*Phuong pháp chia dôi tim nghiem gan dung cua phuong trinh :*)
Clear[a, b, x, e, c]
PPChiaDoi[ss_] :=
Module[{a, b, x, e, c},
a = Input["Nhap vao a= "];
b = Input["Nhap vao b ="];
f[x_] = Input["Nhap vao phuong trinh f(x) = 0"];
e = b - a;
While[e > ss,
(c = (a + b)/2;
If[f[c]*f[a] < 0, b = c, a = c];
e = b - a)];
Return[N[a]]
]
(*Ví du : Phương trình x^3 - x - 1 = 0, khoảng phân ly nghiệm là[1; 2] và sai số 10^-3.
Đáp số : x ≈ 1.324219*)
Print["Ket qua:", PPChiaDoi[0.003]]