Pa napisi funkciju koja ce da ti radi polovljenje intervala. Nemoj da ocekujes da je NSolve svemoguca
Code:
PolovljenjeIntervala[f0_, {x_, a0_, b0_}] := Module[
{f, a, b, c},
a = N[a0];
b = N[b0];
f = Function[x, f0];
c = a + 0.5*(b - a);
While[a < c < b,
If[Sign[f[c]] == Sign[f[a]],
a = c, b = c];
c = a + 0.5*(b - a)];
c];
Eto, ovo ce ti resiti problem.
Pozivas je na sledeci nacin:
Code:
PolovljenjeIntervala[Exp[x] - x - 1, {x, -1, 1}]
i dobijes
Code:
-1.05367*10^-8
a to je 0 u sustini
Lupio sam ove granice, namerno nisam stavio 0, nego oko nje. Znaci ovo a i b su pocetne tacke. Znas vec i sam, da ti ne pricam mnogo. Naravno, funkciju mozes nazvati kako god hoces.
If you can't explain it simply, you don't understand it well enough. A. Einstein