bisseção, algoritmo

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function r=biss(a,b,F,E,n)
 x=sym('x');
 format long
 x1=(a+b)/2;
 for k=1:n
 d=(subs(F,x,a)*(subs(F,x,x1))); 
 if (d<0)
 c=(a+x1)/2;
 x1=c;
 if (abs(a-x1)<E)
 r=double(x1)
 fprintf('\nraiz = %d\n',r)
 return
 end
 if(abs(subs(F,x,x1)))<E
 r=double(x1)
 fprintf('\nraiz = %d\n',r)
 return
 end
 end
 if (d>=0)
 c=(x1+b)/2;
 x1=c;
 if (abs(b-x1))<E
 r=double(x1)
 fprintf('\nraiz = %d\n',r)
 return
 end
 if (abs(subs(F,x,x1)))<E
 r=double(x1)
 fprintf('\nraiz = %d\n',r)
 return
 end
 
 end
 fprintf('\nraiz = %d\n',x1)
 end
end