programa matlab Bicessao

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

%Este programa aplica o método da Bissecção para solução
%de equações não-lineares
close all
clear all
clc
%Dados de entrada
f=@(x)(exp(-0.005*x))*cos((sqrt(2000-x^2/100))*0.05)-0.01;
a=-319;b=-318;
tol=1e-5;
kmax=30;
%Estimativa Inicial
k=1;
x=(a+b)/2;
fx=f(x);
fa=f(a);
fb=f(b);
ER=1e20;
fprintf('k xk ER |f(x)|\n')
fprintf('%d %15.12f %15.12e %15.12e\n',k,x,ER,abs(f(x)))
%Processo Iterativo
while ER > tol
%while abs(f(x)) > tol
 if fa*fx <= 0
 b=x;
 fb=fx;
 else
 a=x;
 fa=fx;
 end
 k=k+1;
 aux=x;
 x=(a+b)/2;
 fx=f(x);
 ER=abs((aux-x)/x);
 fprintf('%d %15.12f %15.12e %15.12e\n',k,x,ER,abs(f(x)))
end