Práctica 4

Contents

Ejercicio 1

Función usada:

x = [1 42 54 66 13 52 7 43 44 1098];

[xpar, ximpar] = SepararParImpar(x);

disp(['Lista: ' num2str(x)])
disp(['Los pares: ' num2str(xpar)])
disp(['Los impares: ' num2str(ximpar)])
Lista: 1    42    54    66    13    52     7    43    44  1098
Los pares: 42    54    66    52    44  1098
Los impares: 1  13   7  43

Ejercicio 2

Funciones usadas:

N = 1000;
[M, xy, id] = PuntosDentroDelCirculo(N);
EstimPi = 4*M/N;

disp([' pi estimado = ' num2str(EstimPi) '; El error relativo = ' ...
    num2str(100*(EstimPi-pi)/pi) '%'])

figure
plot(xy(1,~id),xy(2,~id),'k.','MarkerSize',7);
hold on
plot(xy(1,id),xy(2,id),'r.','MarkerSize',7)
axis([-1.2 1.2 -1.2 1.2])
 pi estimado = 3.036; El error relativo = -3.3611%