I have a simple code for bisection method that isn't working correctly. Can someone please tell me what is wrong.
f[x_] := x*Log[x] + x^2;
Plot[f[x], {x, 0, 1}]
i = 1;
Subscript[a, i] = 0.4;
Subscript[b, i] = 0.6;
Subscript[error, i] = Abs[f[Subscript[x, i]]];
Subscript[error, i] = 1;
While[
Subscript[error, i] >= 0.0001,
Subscript[x, i] = (Subscript[a, i] + Subscript[b, i])/2;
If[
N@f[Subscript[x, i]] > 0,
Subscript[a, i + 1] == Subscript[a, i] &&
Subscript[b, i + 1] == Subscript[x, i],
Subscript[a, i + 1] == Subscript[x, i] &&
Subscript[b, i + 1] == Subscript[b, i]
]
Print["iteration ", i , ", output", ": ", Subscript[x, i], ", ",
Subscript[a, i], ", ", Subscript[b, i]]
i++
]