function sbwinprobs=getsbwinprobs(seeds,p,q)

% seeds(m,i) comtains the indices of the playoff seeds from the two conferences.
% p(j,k)=prob that team j beats team k, if team j is at home.
% q(j,k)=prob that team j beats team k on a neutral field.

pconfchamp=zeros(6,2);  % probs of winning conf championship for the six seeds in each conf.
sbwinprobs=zeros(6,2);  % probs of winning SB for the six seeds in each conference.

for i=1:2

%        djklm is prob that the divisional round has seeds j, k, l, and m.

d1234=p(seeds(3,i),seeds(6,i))*p(seeds(4,i),seeds(5,i));
d1264=(1-p(seeds(3,i),seeds(6,i)))*p(seeds(4,i),seeds(5,i));
d1235=p(seeds(3,i),seeds(6,i))*(1-p(seeds(4,i),seeds(5,i)));
d1265=(1-p(seeds(3,i),seeds(6,i)))*(1-p(seeds(4,i),seeds(5,i)));

%        cjk is prob that conf championship has seeds j and k.

c12=p(seeds(1,i),seeds(4,i))*p(seeds(2,i),seeds(3,i))*d1234 ...
   +p(seeds(1,i),seeds(6,i))*p(seeds(2,i),seeds(4,i))*d1264 ...
   +p(seeds(1,i),seeds(5,i))*p(seeds(2,i),seeds(3,i))*d1235 ...
   +p(seeds(1,i),seeds(6,i))*p(seeds(2,i),seeds(5,i))*d1265;

c13=p(seeds(1,i),seeds(4,i))*(1-p(seeds(2,i),seeds(3,i)))*d1234 ...   
   +p(seeds(1,i),seeds(5,i))*(1-p(seeds(2,i),seeds(3,i)))*d1235;

c24=(1-p(seeds(1,i),seeds(4,i)))*p(seeds(2,i),seeds(3,i))*d1234;

c34=(1-p(seeds(1,i),seeds(4,i)))*(1-p(seeds(2,i),seeds(3,i)))*d1234;

c25=(1-p(seeds(1,i),seeds(5,i)))*p(seeds(2,i),seeds(3,i))*d1235;

c35=(1-p(seeds(1,i),seeds(5,i)))*(1-p(seeds(2,i),seeds(3,i)))*d1235;

c14=p(seeds(1,i),seeds(6,i))*(1-p(seeds(2,i),seeds(4,i)))*d1264;

c26=(1-p(seeds(1,i),seeds(6,i)))*p(seeds(2,i),seeds(4,i))*d1264 ...
   +(1-p(seeds(1,i),seeds(6,i)))*p(seeds(2,i),seeds(5,i))*d1265;

c46=(1-p(seeds(1,i),seeds(6,i)))*(1-p(seeds(2,i),seeds(4,i)))*d1264;

c15=p(seeds(1,i),seeds(6,i))*(1-p(seeds(2,i),seeds(5,i)))*d1265;

c56=(1-p(seeds(1,i),seeds(6,i)))*(1-p(seeds(2,i),seeds(5,i)))*d1265;

%         pconfchamp(j,i) is prob that seed j wins in conference i.

pconfchamp(1,i)=p(seeds(1,i),seeds(2,i))*c12 ...
              +p(seeds(1,i),seeds(3,i))*c13 ...
              +p(seeds(1,i),seeds(4,i))*c14 ...
              +p(seeds(1,i),seeds(5,i))*c15;

pconfchamp(2,i)=(1-p(seeds(1,i),seeds(2,i)))*c12 ...
              +p(seeds(2,i),seeds(4,i))*c24 ...
              +p(seeds(2,i),seeds(5,i))*c25 ...
              +p(seeds(2,i),seeds(6,i))*c26;

pconfchamp(3,i)=(1-p(seeds(1,i),seeds(3,i)))*c13 ...
              +p(seeds(3,i),seeds(4,i))*c34 ...
              +p(seeds(3,i),seeds(5,i))*c35;

pconfchamp(4,i)=(1-p(seeds(2,i),seeds(4,i)))*c24 ...
              +(1-p(seeds(3,i),seeds(4,i)))*c34 ...
              +(1-p(seeds(1,i),seeds(4,i)))*c14 ...
              +p(seeds(4,i),seeds(6,i))*c46;

pconfchamp(5,i)=(1-p(seeds(2,i),seeds(5,i)))*c25 ...
              +(1-p(seeds(3,i),seeds(5,i)))*c35 ...
              +(1-p(seeds(1,i),seeds(5,i)))*c15 ...
              +p(seeds(5,i),seeds(6,i))*c56;

pconfchamp(6,i)=(1-p(seeds(2,i),seeds(6,i)))*c26 ...
              +(1-p(seeds(4,i),seeds(6,i)))*c46 ...
              +(1-p(seeds(5,i),seeds(6,i)))*c56;
           
end
        
for i=1:2
   if i==1
      k=2;
   else
      k=1;
   end
   for j=1:6
      sum=0;
      for m=1:6
         sum=sum+q(seeds(j,i),seeds(m,k))*pconfchamp(m,k);
      end
      sbwinprobs(j,i)=pconfchamp(j,i)*sum;
   end
end