p12=0.05; p21=p12;
p23=0.07; p32=p23;
p34=0.12; p43=p34;
p14=1e-3; p41=p14;
p45=0.2; p54=p45;
p56=0.08; p65=p56;
p16=0.03; p61=p16;


%Paths 1->...->3

P1= 1-(1-p12)*(1-p23)
P2=1-(1-p14)*(1-p43)
P3=1-(1-p16)*(1-p65)*(1-p54)*(1-p43)

%Paths 1-5

P4=1-(1-p12)*(1-p23)*(1-p34)*(1-p45)
     P5=1-(1-p14)*(1-p45)
     P6=1-(1-p16)*(1-p65)

%Paths 3-5
P7=1-(1-p34)*(1-p45)
     P8=1-(1-p32)*(1-p21)*(1-p14)*(1-p45)
     P9=1-(1-p32)*(1-p21)*(1-p16)*(1-p65)
     P10=1-(1-p34)*(1-p41)*(1-p16)*(1-p65)

     disp([P1,P2,P3])
     disp([P4,P5,P6])
     disp([P8,P9,P10])