PROGRAMA CON LAS SALIDAS

[pic 1][pic 2][pic 3][pic 4][pic 5][pic 6]

PROGRAMA CON LAS SALIDAS

%PRACTICA 2_CCALDERON

clc;

clear all;

syms zc;

in1=0;

x=0;

x2=0;

x3=0;

zc2=5;

sum=0;

muest=input('Digite su muestra: ');

disp('     ');

vp=input('Digite valor de prueba: ');

disp('     ');

alfa=input('Digite valor de significancia en forma porcentual: ');

a=alfa/100;

disp('    ');

disp('Conoce el valor de la desviación estandar,ingrese el numero correspindiente a su respuesta:');

disp('1.No     2.Si');

opc=input('  ');

if opc==1

dv=input('Digite valor de desviación estandar: ');

h=(0.2)/400;

for i=1:1:399

    in11=(1/((-0.2+(i*h))^2))*exp((-1)/(2*((-0.2+(i*h))^2)));

    in1=in1+in11;

end

in1=in1*h;

in1=in1+((h/2)*((1/((-0.2)^2))*exp((-1)/(2*((-0.2)^2)))));

vr=((a*(sqrt(2*pi)))/2)-in1;

error=10;

b=-2.5;o

b1=-5;

a1=-5;

b2=0;

while error>0.1

h1=(b-a1)/100000;

for i=1:1:99999

    x2=exp(((-1)*((-5+(h1*i))^2))/(2));

    x3=x3+x2;

end

x3=x3*h1;

x3=x3+((h1/2)*exp(((-1)*(5^2))/(2)))+((h1/2)*(exp((-1*(b^2))/(2))));%suma f(x0)+f(xi)+f(zc)

error=(abs((vr-x3))/vr)*100;

if error>0.1

if vr>x3

b2=b;

b=((abs((b1-b)))/2)+b;

b1=b2;

end

if vr

b2=b;

b=-((abs((b1-b)))/2)+b;

b1=b2;

end

end

end

n=numel(muest);

zc2=b;

for i=1:1:n

    sum1=muest(i);

    sum=sum+sum1;

end

media=sum/n;

min=media-abs((zc2*(dv/sqrt(n))));

max=media+abs((zc2*(dv/sqrt(n))));

graf1=[];

t=[];

t1=[];

t2=[];

gmin=[];

interv=[];

val=[];

for i=-4:0.1:4

    ev=(exp((-1*(i^2))/(2)))/(sqrt(2*pi));

    graf1=[ev graf1];

end

for i=-4:0.1:4

    t=[t i];

end

n3=(4-abs(zc2))/20;

for i=-4:n3:-abs(zc2)

    ev1=(exp((-1*(i^2))/(2)))/(sqrt(2*pi));

    gmin=[gmin ev1];

end

for i=-4:n3:-abs(zc2)

    t1=[t1 i];

end

for i=abs(zc2):n3:4

    ev2=(exp((-1*(i^2))/(2)))/(sqrt(2*pi));

    gmax=[gmax ev2];

end

for i=abs(zc2):n3:4

    t2=[t2 i];

end

n2=(max-min)/20;

for i=min:n2:max

    interv=[interv i];

end

for i=min:n2:max

    val=[1 val];

end

estadistico=(media-vp)/(dv/sqrt(n));

intervalo1=subplot(2,1,1),plot(t,graf1);

title('Función Densidad: ');

hold on

subplot(2,1,1),area(t1,gmin)

hold on

subplot(2,1,1),area(t2,gmax)

text(-abs(zc2),-0.02,num2str(-abs(zc2)));

text(abs(zc2),-0.02,num2str(abs(zc2)));

text(4.2,0,'z');

intervalo1.YAxis.Visible='off';

intervalo1.XAxis.Visible='off';

intervalo=subplot(2,1,2),plot(interv,val),line([min min],[0.8 1.2]),line([max max],[0.8 1.2]);

title('Intervalo de Confiaza');

intervalo.YAxis.Visible='off';

intervalo.XAxis.Visible='off';

text(max,0.78,num2str(max));

text(min,0.78,num2str(min));

hold on

subplot(2,1,2),plot(estadistico,1,'c*')

end

if opc==2

    n=numel(muest);

    v=n+1;

    h=1/400;

    k1=(v+1)/2;

    k2=v/2;

    sk1=0;

    sk2=0;

for i=1:1:399

    sumk1=((h*i)^(k1-1)*exp(-h*i));

...