Trabalho de cálculo numérico - método bissecção no Matlab

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

clear all
clc
fprintf('\n*------------*------------*------------*------------*------------*');
fprintf('\n Claudine Moll -- S41 -- 2017.2');
fprintf('\n Algoritmo para encontrar zeros de funcoes reais');
fprintf('\n Método da bisseccao');
fprintf('\n*------------*------------*------------*------------*------------*');
fprintf('\nFuncao f(x)=x^3+4x^2-10 ');
%Parametros
a=1;
b=2;
epsilon = 0.000001;
k=0;
m=(a+b)/2;
erro=(b-a)/2;
fa=(a*a*a+4*a*a-10);
fm=(m*m*m+4*m*m-10);
while erro > epsilon
 if fa*fm > 0
 a=m;
 else
 b=m;
 end
 
 k=k+1;
 m=(a+b)/2;
 erro=(b-a)/2;
 fa=(a*a*a+4*a*a-10);
 fm=(m*m*m+4*m*m-10);
end
fprintf('\nForam feitas %f iteracoes e o zero da raiz é %f', k+1,m);
fprintf('\n*------------*------------*------------*------------*------------*');