function x=nullstelle(f,a,b,tol) m=(b+a)/2; if f(a)>= 0 | f(b) <=0 disp('f verletzt Bedingung') end while abs(m(end)-a) >= tol if f(m(end)) == 0 x=m(end); break elseif f(m(end))>0 b=m(end); else a=m(end); end m(end+1)=(a+b)/2; end x=m(end); end