Scripts de Cálculo Numérico

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

%ex1 aula 1808
function[raiz] = bissecao(a,b,prec)
F=input('digite a eq. nao linear:');
fa=F(a);
fb=F(b);
if fa*fb >= 0
 disp(['nao existe raiz no intervalo',num2str(a) ',' num2str(b)]);
 raiz='inexistente';
 return
end
erro=10000;
c=(a+b)/2;
cv=c;
while erro >= prec
 fc=F(c);
 if fa*fc < 0
 b=c;
 fb=F(b);
 elseif fb*fc < 0
 a=c;
 fa=F(a);
 else 
 raiz=c;
 return
 end
c=(a+b)/2;
erro=abs((c-cv)/c)*100;
cv=c;
end
raiz=c;