preserva=100;

prenta=200;

for z=1:30

for i=1:30

    camaso(i)=binornd(15,0.4);

end

    for j=1:15

        costo_reserva=j*preserva;

        for i=1:30

            costo_renta=max(camaso(i)-j,0)*prenta;

            costo_opcion(i,j)=costo_renta+costo_reserva;

        end

    end

    for i=1:15

       sumas(i)=sum(costo_opcion(:,i));

   end

        muestra(z,:)=sumas;

    end

    [muc,sigc,cimc,cisic]=normfit(muestra)

    plot((1:15),muc,(1:15),cimc(1,:),(1:15),cimc(2,:))