Phương pháp chia đôi tìm nghiệm gần đúng của phương trình f(x) = 0 với khoảng cách li nghiệm [a,b] sai số gần đúng Epxilon.

Mã:
Clear[f, a, b, c, e, Epxilon];
f[x_] = Input["Nhap f(x)"];
a = Input["Nhap a"];
b = Input["Nhap b"];
Epxilon = Input["Nhap vao sai so gan dung Epxilon"]; 
Print["Giai phuong trinh"];
Print[f[x], " = 0"];
Print["Khoang cach li nghiem [", a, ",", b, "]"];
Print["Sai so gan dung Epxilon = ", Epxilon];

(*Phuong phap chia doi*)
	If[f[a] f[b] <= 0,
    		e = Abs[b - a];
    		While[e > Epxilon,
      			c = (a + b)/2;
      			If[f[a] f[c] < 0, b = c , a = c];
      			e = Abs[b - a];
      		];
    		Print["Nghiem x = ", N[a]];
    	,
    		Print["[", a, ",", b, "] khong phai la khoang cach li nghiem"];
    	];
   	Plot[f[x], {x, 5 a , 5 b}];