function[koren]=bisekce(fce,a,b,iterace);
%resi rovnici metodou bisekce
pravej=b;
levej=a;
f=inline(fce);

for i=1:iterace,
   stred=(levej+pravej)/2;
   hodnota=f(stred);
   if hodnota*f(pravej)>0,
       pravej=stred;
   else
       levej=stred;
       
   end
end
 koren=(levej+pravej)/2;
end