Matlab-bisekcija

[johnko]

A mi lahko en razloži zakaj mi tole ne dela? Sam ne najdem napake...


f=@(x)(x^3-6*x^2-72*x-27);
a=input('Podaj spodnjo mejo intervala: ');
b=input('Podaj zgornjo mejo itervala: ');
if f(a)*f(b)>0
disp('Na tem intervalu ni nicel, vstavi drug interval');
end
eps=input('podaj najmanjso razliko intervala epsilon:');
c=(a+b)/2;
int=(a-b);
stevec=0;
while abs(int)<eps
if f(c)==0
break
end
if f(a)*f(c)<0
b=c;
else
a=c;
end
stevec=stevec+1;
end
fprintf('nicla=%12.f \n',c)
fprintf('število korakov=%12.f \n',stevec)

[Aniviller]

Daj povej kako si ga klical (kaj si dal za a in b), in na kakšen način ti ne dela... kaj je narobe?