function A=SCHNE(n) p=1/4; matrix=0; safe=0; good=0; sched=0; for j=0:2^n-1 for k=0:2^n-1 matrix=matrix+1; x=j; y=k; for j1=n:-1:1 a(j1)=rem(x,2); x=floor(x/2); end for j2=n:-1:1 b(j2)=rem(y,2); y=floor(y/2); end sum=0; jelg=0; jels=0; for j3=1:n sum=sum+a(j3)+b(j3); if sum>j3 jels=1; end if a(j3)+b(j3)>1 jelg=1; end end if jelg==0 val=p^(sum)*(1-p)^(2*n-sum); good=good+val; sched=sched+val; safe=safe+val; end if jelg>0&jels==0 val=p^(sum)*(1-p)^(2*n-sum); safe=safe+val; for i2=1:n+1 for j2=1:n+1 s(i2,j2)=0; end end s(1,1)=1; for k3=1:2n for i2=1:n for j2=1:n if (s(i2,j2)==1)&(a(i2)==0) s(i2+1,j2)=1; end if (s(i2,j2)==1)&(b(j2)==0) s(i2,j2+1)=1; end if (s(i2,j2)==1)&(a(i2)+b(j2)==1) s(i2+1,j2+1)=1; end end end end board=0; for i3=1:n+1 board=board+s(n+1,i3); end for i4=1:n board=board+s(i4,n+1); end if board>0 sched=sched+val; end end end end disp([matrix,good,sched,safe]); end