pmi=0.5;
pmizo=0.8;
pmizf=0.4;

pmiz1= 0.8*0.8*0.4*0.4;
pmiz2=pmiz1;
pmiz4=pmiz1;
pmiz5=pmiz1;
pmiz3= 0.8^3*0.4;
pmizs=[pmiz1,pmiz2,pmiz3,pmiz4,pmiz5];

pmizl= pmizs(1)
for i = 2:length(pmizs)
    i
    pmizc= 1/(1+((1-pmizs(i))/(pmizs(i))*(1-pmizl)/(pmizl)* pmi/(1-pmi)));
    pmizl= pmizc
end

pmizl2