initial commit / add all books

38 files changed, 209 insertions, 0 deletions

diff --git a/1475/CH2/EX2.10/Example_2_10.sce b/1475/CH2/EX2.10/Example_2_10.sce
new file mode 100755
index 000000000..104de24ed
--- /dev/null
+++ b/1475/CH2/EX2.10/Example_2_10.sce
@@ -0,0 +1,14 @@
+//Example 2.10 the overall percentage of failure    

+clc;

+clear;

+function value = binomial(n, k, p)

+value = factorial(n)*(p^k)*((1-p)^(n-k))/(factorial(k)*factorial(n-k));

+endfunction

+q=40/100;

+n=6;

+p=1-q;

+disp(n,"No. of candidates =",p,"Proabab. of(success) in a single trial",q,"Probab. of failure of a candidate");

+P_4=binomial(n,4,p);

+P_5=binomial(n,5,p);

+P_6=binomial(n,6,p);

+disp(P_4+P_5+P_6,"Required Probability =");

diff --git a/1475/CH2/EX2.10/Output_2_10.PNG b/1475/CH2/EX2.10/Output_2_10.PNG
new file mode 100755
index 000000000..fa01d0f46
--- /dev/null
+++ b/1475/CH2/EX2.10/Output_2_10.PNG
diff --git a/1475/CH2/EX2.13/Example_2_13.sce b/1475/CH2/EX2.13/Example_2_13.sce
new file mode 100755
index 000000000..755a224df
--- /dev/null
+++ b/1475/CH2/EX2.13/Example_2_13.sce
@@ -0,0 +1,8 @@
+// Example 2.13 Arithmetic mean and standard deviation of a binomial distribution

+clc;

+clear;

+np=4;

+npq=8/3;

+q=npq/np;

+p=1-q;

+disp(q,"Value of q",p,"Value of p");

diff --git a/1475/CH2/EX2.13/Output_2_13.PNG b/1475/CH2/EX2.13/Output_2_13.PNG
new file mode 100755
index 000000000..475f762ae
--- /dev/null
+++ b/1475/CH2/EX2.13/Output_2_13.PNG
diff --git a/1475/CH2/EX2.15/Example_2_15.sce b/1475/CH2/EX2.15/Example_2_15.sce
new file mode 100755
index 000000000..d4f84ca14
--- /dev/null
+++ b/1475/CH2/EX2.15/Example_2_15.sce
@@ -0,0 +1,19 @@
+// Example 2.15 A random variable x follow Poisson distribution

+clc;

+clear;

+m=3;

+disp(m,"Parameter =");

+function value = poisson(m, x)

+value = (exp(-m)*(m^x))/(factorial(x));

+endfunction

+X_0=poisson(m,0);

+X_1=poisson(m,1);

+X_2=poisson(m,2);

+X_3=poisson(m,3);

+

+disp(X_0,"Probability when x=0 is ");

+disp(X_1,"Probability when x=1 is ");

+disp(X_2,"Probability when x=2 is ");

+disp(X_3,"Probability when x=3 is ");

+disp(X_0+X_1+X_2,"Probability when the variable assumes the value less than 3 is ");

+disp(1-(X_0+X_1),"Probability when the variable assumes the value at least 2 is ");

diff --git a/1475/CH2/EX2.15/Output_2_15.PNG b/1475/CH2/EX2.15/Output_2_15.PNG
new file mode 100755
index 000000000..88dd40f5d
--- /dev/null
+++ b/1475/CH2/EX2.15/Output_2_15.PNG
diff --git a/1475/CH2/EX2.16/Example_2_16.sce b/1475/CH2/EX2.16/Example_2_16.sce
new file mode 100755
index 000000000..a1256075c
--- /dev/null
+++ b/1475/CH2/EX2.16/Example_2_16.sce
@@ -0,0 +1,12 @@
+// Example 2.16 A business firm receives on an average 

+clc;

+clear;

+m=2.5;

+disp(m,"Average call =");

+function value = poisson(m, x)

+value = (exp(-m)*(m^x))/(factorial(x));

+endfunction

+X_0=poisson(m,0);

+X_4=poisson(m,4);

+disp(X_0,"Probability of no call is ");

+disp(X_4,"Probability of exactly 4 call is ");

diff --git a/1475/CH2/EX2.16/Output_2_16.PNG b/1475/CH2/EX2.16/Output_2_16.PNG
new file mode 100755
index 000000000..c1e497003
--- /dev/null
+++ b/1475/CH2/EX2.16/Output_2_16.PNG
diff --git a/1475/CH2/EX2.2/Example_2_2.sce b/1475/CH2/EX2.2/Example_2_2.sce
new file mode 100755
index 000000000..74c0f2102
--- /dev/null
+++ b/1475/CH2/EX2.2/Example_2_2.sce
@@ -0,0 +1,10 @@
+// Example 2.2 A random variable has the following proabability distribution

+clc;

+clear;

+x=[4 5 6 8];

+Probab=[0.1 0.3 0.4 0.2];

+disp(sum(x .* Probab),"Mean =");

+disp(sum((x .^2) .* Probab), "Expectation =");

+disp(sum((x .^2) .* Probab)-(sum(x .* Probab)^2), "Variance =");

+disp(sqrt(sum((x .^2) .* Probab)-(sum(x .* Probab)^2)),"Standard Deviation = ");

+

diff --git a/1475/CH2/EX2.2/Output_2_2.PNG b/1475/CH2/EX2.2/Output_2_2.PNG
new file mode 100755
index 000000000..5852f2c7e
--- /dev/null
+++ b/1475/CH2/EX2.2/Output_2_2.PNG
diff --git a/1475/CH2/EX2.20/Example_2_20.sce b/1475/CH2/EX2.20/Example_2_20.sce
new file mode 100755
index 000000000..a827bc628
--- /dev/null
+++ b/1475/CH2/EX2.20/Example_2_20.sce
@@ -0,0 +1,11 @@
+// Example 2.20 In turning out certain toys

+clc;

+clear;

+n=10;

+p=0.1;

+m=n*p;

+function value = poisson(m, x)

+value = (exp(-m)*(m^x))/(factorial(x));

+endfunction

+X_3=poisson(m,3);

+disp(X_3,"Probability of 3 defective in the sample is ",);

diff --git a/1475/CH2/EX2.20/Output_2_20.PNG b/1475/CH2/EX2.20/Output_2_20.PNG
new file mode 100755
index 000000000..05a1ed353
--- /dev/null
+++ b/1475/CH2/EX2.20/Output_2_20.PNG
diff --git a/1475/CH2/EX2.21/Example_2_21.sce b/1475/CH2/EX2.21/Example_2_21.sce
new file mode 100755
index 000000000..ffaa2ac73
--- /dev/null
+++ b/1475/CH2/EX2.21/Example_2_21.sce
@@ -0,0 +1,14 @@
+// Example 2.21 2 percent of items made by a machine are defective

+clc;

+clear;

+n=100;

+p=2/100;

+m=n*p;

+disp(m,"Mean is n*p =",p,"probab. of defective items",n,"Total sample of item");

+function value = poisson(m, x)

+value = (exp(-m)*(m^x))/(factorial(x));

+endfunction

+X_0=poisson(m,0);

+X_1=poisson(m,1);

+X_2=poisson(m,2);

+disp(1-(X_0 + X_1 + X_2),"Probability of 3 or more defective in the sample is ",);

diff --git a/1475/CH2/EX2.21/Output_2_21.PNG b/1475/CH2/EX2.21/Output_2_21.PNG
new file mode 100755
index 000000000..49bdfac13
--- /dev/null
+++ b/1475/CH2/EX2.21/Output_2_21.PNG
diff --git a/1475/CH2/EX2.22/Example_2_22.sce b/1475/CH2/EX2.22/Example_2_22.sce
new file mode 100755
index 000000000..db8410d0a
--- /dev/null
+++ b/1475/CH2/EX2.22/Example_2_22.sce
@@ -0,0 +1,20 @@
+// Example 2.22 Two defective tube lights are mixed

+clc;

+clear;

+n=3;

+N=10;

+p=2/10;

+m=n*p;

+q=1-p;

+var=n*p*q;

+S_d=sqrt(var);

+Np=N*p;

+Nq=N*q;

+disp(S_d,"Standard Deviation",var,"Variance",m,"Mean is n*p =",p,"probab. of defective items",N,"Total sample of item");

+function value = hyper(Np,Nq,N,n, x)

+value = (factorial(Np)*factorial(Nq)*factorial(N-n)*factorial(n))/(factorial(Np-x)*factorial(x)*factorial(n-x)*factorial(Nq-(n-x))*factorial(N));

+endfunction

+f_0=hyper(Np,Nq,N,n,0);

+f_1=hyper(Np,Nq,N,n,1);

+f_2=hyper(Np,Nq,N,n,2);

+disp(f_0,"Probability of 0 defective tube light in the sample of 3 is ",f_1,"Probability of 1 defective tube light in the sample of 3 is ",f_2,"Probability of 2 defective tube light in the sample of 3 is ",);

diff --git a/1475/CH2/EX2.22/Output_2_22.PNG b/1475/CH2/EX2.22/Output_2_22.PNG
new file mode 100755
index 000000000..87b8a240c
--- /dev/null
+++ b/1475/CH2/EX2.22/Output_2_22.PNG
diff --git a/1475/CH2/EX2.23/Example_2_23.sce b/1475/CH2/EX2.23/Example_2_23.sce
new file mode 100755
index 000000000..b7a69841a
--- /dev/null
+++ b/1475/CH2/EX2.23/Example_2_23.sce
@@ -0,0 +1,8 @@
+//Example 2.23 A dice is rolled 6 times

+clc;

+clear;

+p=1/6;

+n=6;

+x=1;

+probab=(factorial(n)*p*p*p*p*p*p)/(factorial(x)*factorial(x)*factorial(x)*factorial(x)*factorial(x)*factorial(x));

+disp(probab,"the required probability",n,"No. of independent trials");

diff --git a/1475/CH2/EX2.23/Output_2_23.PNG b/1475/CH2/EX2.23/Output_2_23.PNG
new file mode 100755
index 000000000..48292cad5
--- /dev/null
+++ b/1475/CH2/EX2.23/Output_2_23.PNG
diff --git a/1475/CH2/EX2.24/Example_2_24.sce b/1475/CH2/EX2.24/Example_2_24.sce
new file mode 100755
index 000000000..f9fe4626d
--- /dev/null
+++ b/1475/CH2/EX2.24/Example_2_24.sce
@@ -0,0 +1,11 @@
+//Example 2.24 A box contains 4 white, 3 black and 5 red ball

+clc;

+clear;

+p1=1/3;

+p2=1/6;

+p3=5/12;

+n=6;

+probab1=(factorial(n)*(p1^3)*(p2^1)*(p3^2))/(factorial(3)*factorial(1)*factorial(2));

+probab2=(factorial(n)*(p1^2)*(p2^0)*(p3^4))/(factorial(2)*factorial(0)*factorial(4));

+disp(probab2,"the required probability of 2 are white and 4 are red",probab1,"the required probability of 3 are white , 1 black and 2 red",n,"No. of independent trials");

+

diff --git a/1475/CH2/EX2.24/Output_2_24.PNG b/1475/CH2/EX2.24/Output_2_24.PNG
new file mode 100755
index 000000000..bcb99a55d
--- /dev/null
+++ b/1475/CH2/EX2.24/Output_2_24.PNG
diff --git a/1475/CH2/EX2.35/Example_2_35.sce b/1475/CH2/EX2.35/Example_2_35.sce
new file mode 100755
index 000000000..efc3959d9
--- /dev/null
+++ b/1475/CH2/EX2.35/Example_2_35.sce
@@ -0,0 +1,7 @@
+// Example 2.35 If a continuous random variable x follows rectangular distribution

+clc;

+clear;

+disp((4-2.5)/(7-2),"P(2.5<=x<=4) = (4-2.5)/(7-2)=");

+disp((3.4-2)/(7-2),"P(x<=3.4)= (3.4-2)/(7-2)=");

+disp(1-((6-2)/(7-2)),"P(x>6)=1-P(x<=6)= 1-((6-2)/(7-2))=");

+disp(0,"P(x=4.5)=0");

diff --git a/1475/CH2/EX2.35/Output_2_35.PNG b/1475/CH2/EX2.35/Output_2_35.PNG
new file mode 100755
index 000000000..e84557514
--- /dev/null
+++ b/1475/CH2/EX2.35/Output_2_35.PNG
diff --git a/1475/CH2/EX2.36/Example_2_36.sce b/1475/CH2/EX2.36/Example_2_36.sce
new file mode 100755
index 000000000..42e11b112
--- /dev/null
+++ b/1475/CH2/EX2.36/Example_2_36.sce
@@ -0,0 +1,7 @@
+// Example 2.36 The height distribution of a group of 10,000 men

+clc;

+clear;

+disp(4.5,"s.d =",64.5,"Mean =");

+disp("Hence the total no. of men whose height between 55.5 and 69 inch is 82% of 10,000 i.e 82000",(55.5-64.50)/4.5,"and",(69-64.5)/4.5,"Percentage of men whose height lies between 55.5 and 69 inch = Area under standard normal curve between the vertical lines at the corresponding values of ");

+disp("Hence the total no. of men whose height is less than 55.5 is 20% of 10,000 i.e 2000",(55.5-64.50)/4.5,"Percentage of men whose height is less than 55.5 inch = Area under standard normal curve between the vertical lines at the corresponding values of ");

+disp("Hence the total no. of men whose height is more than 73.5 is 20% of 10,000 i.e 2000",(73.5-64.50)/4.5,"Percentage of men whose height is more than 73.5 = Area under standard normal curve between the vertical lines at the corresponding values of ");

diff --git a/1475/CH2/EX2.36/Output_2_36.PNG b/1475/CH2/EX2.36/Output_2_36.PNG
new file mode 100755
index 000000000..0f93e3f24
--- /dev/null
+++ b/1475/CH2/EX2.36/Output_2_36.PNG
diff --git a/1475/CH2/EX2.37/Example_2_37.sce b/1475/CH2/EX2.37/Example_2_37.sce
new file mode 100755
index 000000000..31e0399a4
--- /dev/null
+++ b/1475/CH2/EX2.37/Example_2_37.sce
@@ -0,0 +1,7 @@
+// Example 2.36 The height distribution of a group of 10,000 men

+clc;

+clear;

+disp(15,"s.d =",151,"Mean =");

+disp("Hence the total no. of student whose weight lies between 120 and 155 lbs is 0.5872*500 i.e 294",(120-151)/15,"and",(155-151)/15,"Proportion of student  whose weight lies between 120 and 155 lbs  = Area under standard normal curve between the vertical lines at the corresponding values of ");

+disp("Hence the total no. of student whose weight is more than  155 lbs is 0.3964*500 i.e 197",(155-151)/15,"Proportion of student  whose weight is more than 155 lbs  = Area under standard normal curve between the vertical lines at the corresponding values ofis (155-151)/15 ");

+

diff --git a/1475/CH2/EX2.37/Output_2_37.PNG b/1475/CH2/EX2.37/Output_2_37.PNG
new file mode 100755
index 000000000..e32723451
--- /dev/null
+++ b/1475/CH2/EX2.37/Output_2_37.PNG
diff --git a/1475/CH2/EX2.38/Example_2_38.sce b/1475/CH2/EX2.38/Example_2_38.sce
new file mode 100755
index 000000000..ef9c79121
--- /dev/null
+++ b/1475/CH2/EX2.38/Example_2_38.sce
@@ -0,0 +1,7 @@
+// Example 2.38 A sample of 100 dry battery 

+clc;

+clear;

+//disp(3,"s.d =",12,"Mean =");

+disp((0.5-0.3413)*100,"Hence percentage of battery cells expected to have life is (0.5-0.3413)*100",(15-12)/3,"Proportion of battery cell with life more than 15 hours= Area under standard normal curve between the vertical lines at the corresponding values of (15-12)/3 = ");

+disp((0.5-0.4772)*100,"Hence percentage of battery cells expected to have life is (0.5-0.4772)*100",(6-12)/3,"Proportion of battery cell with life less than 6 hours= Area under standard normal curve between the vertical lines at the corresponding values of (6-12)/3 = ");

+disp((2*0.2487)*100,"Hence percentage of battery cells expected to have life is (2*0.2487)*100",(10-12)/3,"and",(14-12)/3,"Proportion of battery cell with life between 10 and 14 hours= Area under standard normal curve between the vertical lines at the corresponding values of (10-12)/3 and (14-12)/3= ");

diff --git a/1475/CH2/EX2.38/Output_2_38.PNG b/1475/CH2/EX2.38/Output_2_38.PNG
new file mode 100755
index 000000000..02d2bf422
--- /dev/null
+++ b/1475/CH2/EX2.38/Output_2_38.PNG
diff --git a/1475/CH2/EX2.42/Example_2_42.sce b/1475/CH2/EX2.42/Example_2_42.sce
new file mode 100755
index 000000000..d167b1a36
--- /dev/null
+++ b/1475/CH2/EX2.42/Example_2_42.sce
@@ -0,0 +1,11 @@
+// Example 2.42 A fair coin is tossed 400 times

+clc;

+clear;

+n=400;

+p=1/2;

+q=1-p;

+m=n*p;

+s_d=sqrt(n*p*q);

+disp(2*0.199," = 2*(Area between z= 0 and z= 0.05)",(199.5-m)/s_d,"and",(200.5-m)/s_d,"Probability of obtaining exactly 200 heads  = (199.5-m)/s_d and (200.5-m)/s_d =");

+disp(0.5+0.3289," = (Area between left of z= 0 + Area between z=0 and z= 0.95)",(209.5-m)/s_d,"Probability of obtaining less than 210 heads =  (209.5-m)/s_d and ");

+disp(2*0.3531," = 2*(Area between z= -1.05 and z= 1.05)",(189.5-m)/s_d,"and",(210.5-m)/s_d,"Probability of obtaining betweeen 190 and 210 heads = (189.5-m)/s_d and (210.5-m)/s_d =");

diff --git a/1475/CH2/EX2.42/Output_2_42.PNG b/1475/CH2/EX2.42/Output_2_42.PNG
new file mode 100755
index 000000000..e7a4bc400
--- /dev/null
+++ b/1475/CH2/EX2.42/Output_2_42.PNG
diff --git a/1475/CH2/EX2.43/Example_2_43.sce b/1475/CH2/EX2.43/Example_2_43.sce
new file mode 100755
index 000000000..421cd397d
--- /dev/null
+++ b/1475/CH2/EX2.43/Example_2_43.sce
@@ -0,0 +1,6 @@
+// Example 2.43 A variable x follows Poisson distribution

+clc;

+clear;

+m=16;

+s_d=sqrt(m);

+disp(1-0.8092," = 1-(Area of under standard normal curve )",(19.5-m)/s_d,"P(x>=20 of Poisson Distribution)  = (19.5-m)/s_d =");

diff --git a/1475/CH2/EX2.43/Output_2_43.PNG b/1475/CH2/EX2.43/Output_2_43.PNG
new file mode 100755
index 000000000..e7c8f07f1
--- /dev/null
+++ b/1475/CH2/EX2.43/Output_2_43.PNG
diff --git a/1475/CH2/EX2.6/Example_2_6.sce b/1475/CH2/EX2.6/Example_2_6.sce
new file mode 100755
index 000000000..fe1ac5d75
--- /dev/null
+++ b/1475/CH2/EX2.6/Example_2_6.sce
@@ -0,0 +1,5 @@
+// Example 2.6 If a discrete random varaiable x follows

+clc;

+clear;

+disp(0,"Probab. that it takes a value other than the given",1/6,"Probab. that it takes any particular value of the given","Given that it assumes value 8 , 9 , 11 , 15 , 18 , 20 ");

+disp((1/6)*3,"P(|x-14|<5)=",(1/6)*2,"P(x>15)=",(1/6)*4,"P(x<=15) =",(1/6)*3,"P(x<15)=",0,"P(x=12)",1/6,"P(x=9)");

diff --git a/1475/CH2/EX2.6/Output_2_6.PNG b/1475/CH2/EX2.6/Output_2_6.PNG
new file mode 100755
index 000000000..39b6f69a2
--- /dev/null
+++ b/1475/CH2/EX2.6/Output_2_6.PNG
diff --git a/1475/CH2/EX2.8/Example_2_8.sce b/1475/CH2/EX2.8/Example_2_8.sce
new file mode 100755
index 000000000..c5356260b
--- /dev/null
+++ b/1475/CH2/EX2.8/Example_2_8.sce
@@ -0,0 +1,14 @@
+// Example 2.8 Three coins are tossed

+clc;

+clear;

+function value = binomial(n, k, p)

+value = factorial(n)*(p^k)*((1-p)^(n-k))/(factorial(k)*factorial(n-k));

+endfunction

+p=1/2;

+n=3;

+disp(n,"No. of independent trials =",p,"Proabab. of success in a single trial");

+P_0=binomial(n,0,p);

+P_1=binomial(n,1,p);

+P_2=binomial(n,2,p);

+P_3=binomial(n,3,p);

+disp(1-(P_0),"Prob. of at least one success",P_2 + P_3,"Prob of more than one success",P_3,"Prob. of 3 success",P_2, "Prob. of 2 success",P_1,"Prob. of 1 success",P_0,"Prob. of 0 success");

diff --git a/1475/CH2/EX2.8/Output_2_8.PNG b/1475/CH2/EX2.8/Output_2_8.PNG
new file mode 100755
index 000000000..8a1f89471
--- /dev/null
+++ b/1475/CH2/EX2.8/Output_2_8.PNG
diff --git a/1475/CH2/EX2.9/Example_2_9.sce b/1475/CH2/EX2.9/Example_2_9.sce
new file mode 100755
index 000000000..4b6bc5baa
--- /dev/null
+++ b/1475/CH2/EX2.9/Example_2_9.sce
@@ -0,0 +1,18 @@
+//Example 2.9 Five coins are tossed 3200 times      

+clc;

+clear;

+function value = binomial(n, k, p)

+value = factorial(n)*(p^k)*((1-p)^(n-k))/(factorial(k)*factorial(n-k));

+endfunction

+p=1/2;

+n=5;

+q=1/2;

+disp(n,"No. of coins tossed =",p,"Proabab. of getting head (success) in a single trial");

+P_0=binomial(n,0,p);

+P_1=binomial(n,1,p);

+P_2=binomial(n,2,p);

+P_3=binomial(n,3,p);

+P_4=binomial(n,4,p);

+P_5=binomial(n,5,p);

+EF=sum([P_0 P_1 P_2 P_3 P_4 P_5]*3200);

+disp(sqrt(n*p*q),"Standard Deviation =",n*p,"Mean =",EF,"Efective frequency =");

diff --git a/1475/CH2/EX2.9/Output_2_9.PNG b/1475/CH2/EX2.9/Output_2_9.PNG
new file mode 100755
index 000000000..cdbdbf2d1
--- /dev/null
+++ b/1475/CH2/EX2.9/Output_2_9.PNG