8/18/2019 metodos labo

1/4

Cls; clear all;x=0:0.01:pi; %Creamos un vector con valores entre 0 y pi.y=sin(3.*x)+cos(x);subplot(1!1) % "ivi#imos la pantalla en #os secciones y $racamos en la primeraseccion

plot(xy)title(&'(x) = seno(3x) + cos(x)&)$ri#x=3:0.01:3; % "enimos el ran$o #e valores para $racar la se$un#a 'unciony=exp(x) x !;subplot(1!!)plot(xy)title (&'(x) = exp(x) x !&)$ri# %raamos la cua#ricula

clcclear allc=0.1:.01:1;'=,.-1.*0./c .*(1 exp(c.*10/0))0;plot(c')$ri#title (&'=,.-1.*0./c .*(1 exp(c.*10/0))0&)

 = cos(x)*cos2(x)+1 en el intervalo 0!04.5l si$uiente $uion escrito en 6at7ab permite obtener la anterior $ura:x=0:0.01:!0;y = cos(x).*cos2(x)+1;plot(xy)$ri#subplot(1!1)plot(xy)$ri#subplot(1!!)plot(xy)

89(0!0!0!04)

8/18/2019 metodos labo

2/4

% ?ro$rama para calcular raices #e ecuaciones basa#o en el meto#o #ebiseccion%'unction biseccionclc;clear;

a=input(&@alor #e la cota in'erior: &);b=input(&@alor #e la cota superior: &);

 A7=input(&olerancia #e error : &);maxBiter=input(&umero #e iteracciones: &);iter = 0;'print'(&Dn65A"A "5 E5CCA...Dn5F78"ADn&)'print'(&ter a 9r b '(a) '(xr) 5rrorDn&);G2ile 1 %nicio un cliclo innito #el cual sal$o al 2allar la rai o exce#ermaxBiterxr=(b+a)/!; %8proximo el valor #e la posible rai como el punto me#io #elintervalo.'a = 'uncion(a); %5valuamos la 'uncion en la cota in'erior

'b = 'uncion(b); %evaluamos la 'uncion en la cota superior'xr= 'uncion(xr); %5valuamos la 'uncion en la posible raierror = abs((b a)/a);

% mpresion #e resulta#os en ca#a iteracion'print'(&%.0'%10.'%10.'%10.'%10.'%10.'%1!.'Dn&iteraxrb'a'xrerror)i' error H= A7 %i supero el error acepta#o #eten$o el proceso'print'(&?roceso conclui#o exitosamente con el nivel #e error H=%10.'DnDn&A7)breaI;en#

i' abs('xr)H=A7 %e 2allo la aproximacion a la rai y sal$o #el ciclo innito'print'(&Dn?roceso conclui#o exitosamente con el nivel #e error H=%1!.eDnDn&A7)breaI;en#i' 'a*'xr H 0b=xr; %7a rai esta en el primer intervaloelsea=xr; %7a rai esta en el se$un#o intervaloen#i' (iter J maxBiter) %@erico si se exce#e el numero #e iteracciones'print'(&Dnumero #e iteracciones exce#i#o...DnDn&)

breaI;en#iter=iter+1; %ncremento el numero #e iteraccionesen#'print'(&ai aproxima#a: %1!.'&xr);'print'(& teraciones : %.0'Dn&iter);% 6o#elo $racox=Cotan':0.01:Cotaup;'='uncion(x);

8/18/2019 metodos labo

3/4

plot(x')$ris'unction y='uncion(x) %8 continuaciKn #ebo colocar la 'unciKn a evaluary=cos(x);% n #el pro$rama

%?ro$rama para calcular raices #e ecuaciones basa#o en el al$oritmo #e Lalsa%?osicion'unction Lalsa?osicionclc;clear all;a=input(&@alor #e la cota in'erior: &);b=input(&@alor #e la cota superior: &);

 A7=input(&olerancia #e error : &);maxBiter=input(&umero #e iteracciones: &);%6o#elo $raco

x=a:0.01:b;'='uncion(x);plot(xxx');$ri#;title (&

8/18/2019 metodos labo

4/4

en#i' (iterJmaxBiter) %@erico si se exce#e el numero #e iteracciones'print'(&Dnumero #e iteracciones exce#i#o...DnDn&)breaI;en#iter=iter+1; %ncremento el numero #e iteracciones

en#'print'(&Dnai aproxima#a: %1!.'&xr);'print'(& teraciones: %.0'Dn&iter);'unction y='uncion(x) %8ca es #on#e #obo colocar la 'uncion a evaluary=x.M3+!.*x.M!+10.*x!0;%Lin #el pro$rama