2606/CH6/EX6.17/ex6_17.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25

//Page Number: 6.21
//Example 6.17
clc;
//Given
P1=0.6;
P0=0.4;
n=5; //Five digit sequence
j=2; //two outcomes 0 and 1

//(a)1,1,0,0,0
xf=(factorial(n))/((factorial(j))*(factorial(n-j)));
s=xf*(P1^j)*(P0^(n-j));
disp(s,'Probability for 1,1,0,0,0:');

//(b)atleast 3 1's
//P(X>=3)=1-P(X<=2)
//Here y=1-x
x=0;
for k=0:2
    f=(factorial(n))/((factorial(k))*(factorial(n-k)));
    x=x+{f*((P1^k)*(P0^(n-k)))};
end
y=1-x;
disp(y,'Probability for atleast three 1 s:');