Rešitev enačbe

Besedilo naloge je:

Z metodo bisekcije poiščite na šest decimalk točno vse tri rešitve enačbe f(x) = 0, kjer je funkcija podana s predpisom

f(x) = 3.39 sin (x − 0.57) ; 0 < x < Pi (stroga enačaja)

S sledečo metodo mi je uspelo izračunati dve drugi funkciji, ki sta bili tudi podani, pri tej pa se reč v zanki zacikla:

format long
f=inline('3.39*sin(x-0.57)','x')
a=-0
b=-pi
c=(a+b)/2

while abs(f(c))>1e-20
if f(a)*f(c)<0
b=c;

else
a=c
end
c=(a+b)/2
end

c

Kak namig?

Računam v Octave/Matlab.

Sicer pa ne vem, kako Mathlab deluje - ali lahko narediš cca 100 iteracij, pa si izpišeš a, b in f(a), f(b)? Ne vem, kako dober je oni long; npr Excel se "zacikla", ko najde ničlo pri -2,571592654 - ker bolj natančno pač ne zna (zanjga je na tem mestu a=b in od tu dalje ne gre, f(a) pa je cca 10^-16).