Reporte

INSTITUTO TECNOLOGICO DE POCHUTLA

ING. EN SISTEMAS COMPUTACIONALES

TRABAJO:
USO DE HERRAMIENTAS COMPUTACIONALES

4 SEMESTRE
INTEGRANTES DE EQUIPO:
FLAVIO PEDRO SANTIAGO
IRVING CASTILLO PEREZ
BERTIN RAMIREZ VALENCIA
LUCINA JIJON MARTINEZ
OLAY EDWIN MARTINEZ MARTINEZ

MATERIA:
METODOS NUMERICOS

UNIDAD 2


SAN PEDRO POCHUTLA, MARZO 2012

CODIGO DEMETODO DE BISECCION
%****************************************************************
clear; %borra las variables en el espacio de trabajo
clc; %**limpia comandos en la ventana
%****************************************************************
%** Se ingresan los datos de la ecuación que deseamos **
%** conocer sus raíces, entre que intervalosy con que **
%** porcentaje de error. **
%****************************************************************
fprintf('\cálculo de la raíz de una ecuación por método de Bisección\n\n'); %**muestra en pantalla el texto solicitado con ' '
y=input('Dame la función :','s'); %** y guardara la función solicitada
xl=input('Dame el intervalo inferior : '); %**xl o intervalo inicial guardara el valor solicitado
xu=input('Dame el intervalo superior : '); %**xu o intervalo ultimo almacenara el valor solcitado
e=input('Dame el porciento del error : '); %**e o porcentaje de error almacenara el error maximo al que se quiere calcular la iteraciones
xi=0; %**xi seiguala a cero
ea=100; %**ea o error absoluto igual a cien
c=1; %**contador igual a 1
%****************************************************************
%** Se hace el calculo para determinar si en los intervalos **
%** para los que se va a calcular esta comprendida la raiz **
%** de la ecuación**

%****************************************************************
x=xl; %**x es igual a intervalo inicial
a=eval(y);
x=xu; %**x es igual a intervalo ultimo
b=eval(y);
cc=a*b; %**cc es igual a a por b
if cc>0 %** si cc es mayor que cero, hacer...
fprintf('\n\nLos intervalos que haingresado para hacer el calculo de la raiz\n\n'); %**imprime en pantalla el texto solicitado
fprintf('\nde la ecuacion no son los adecuados por no existir una raíz entre \n\n'); %**imprime en pantalla el texto solicitado
fprintf('\nestos, Execute de nuevo el programa, por favor.\n\n'); %**imprime en pantalla el texto solicitado
break %**rompe el bucle while
end %**termina la sentencia%****************************************************************
While ea>e %** mientras error absoluto sea menor que error, hacer
xr=(xl+xu)/2; %** raiz nueva es igual a (intervalo inicial + intervalo ultimo) sobre dos
x=xl; %** x es igual a intervalo inicial
yl=eval(y);
x=xr; %** x es igual a raiz nueva
yr=eval(y);
ea=abs((xr-xi)/xr)*100; %error absoluto es igual a valor absoluto de((raiz nueva-intervalo inicial)/raiz nueva) * 100
z=yl*yr;
if z<0 %** si z es menor que cero, hacer...
xu=xr; %**intervalo ultimo es igual a raiz nueva
end %**finaliza condicion
if z>0 %**si z es mayor que cero, hacer...
xl=xr; %**intervalo inicial es igual a raiz nueva
end %**finaliza condicion
if z==0 %**si z es igual a cero, hacer...
fprintf('\n\n\n\nLa raíz exacta es: %d',xr)%** imprime pantalla el texto solicitado y muestra la raiz nueva
fprintf('\n\nNumero de iteraciones: %d',c); %** imprime pantalla el texto solicitado y muestra el contador o numero de iteraciones
break %**finaliza el bucle while
end %** finaliza el bucle en ejecucion
xi=xr; %** intervalo inicial es igual a raiz nueva
c=c+1; %** contador es igual a contador mas uno
end %** finaliza...