function svar=intervall(a,b,fel) 
%svar=intervall(a,b,fel)
% intervallhalvering    [a b] startintervall
% funktionen vars nollställen bestäms heter fex
while b-a>fel
m=(a+b)/2;
if fex(a)*fex(m)<0
b=m;
else a=m;
end
[a b]
end
svar=[a b];
