18 files changed, 844 insertions, 0 deletions
diff --git a/3513/CH6/EX6.1/Ex6_1.sce b/3513/CH6/EX6.1/Ex6_1.sce
new file mode 100644
index 000000000..858c979a3
--- /dev/null
+++ b/3513/CH6/EX6.1/Ex6_1.sce
@@ -0,0 +1,48 @@
+// percentage of rework
+//page no 103
+clear
+clc;
+x=41.283;
+r=0.280;
+k=20;
+ucl=2.07;        //Upper specification limit USL
+lcl=2.03 ;       //Lower specification limit LSL
+A2 = 0.73;
+d2 = 2.059;
+D3 = 0.0;
+D4 = 2.28;
+X1=x/k;
+R1=r/k;
+
+//a) Control limit for R – chart
+ucl1 = D4*R1;
+lcl1 = D3*R1;
+cl1=R1;
+mprintf("\nControl limit for R – chart\n ucl = %.2f ",ucl1);
+mprintf("\nlcl = %.2f ",lcl1);
+mprintf("\nCL = %.2f ",cl1);
+//Control limits for X – chart
+ucl2=X1+A2*R1;
+lcl2=2.05393;
+cl2=X1;
+mprintf("\nControl limit for X – chart\n ucl = %.2f ",ucl2);
+mprintf("\nlcl = %.2f ",lcl2);
+mprintf("\nCL = %.2f ",cl2);
+
+//(b) Process Capability
+cl3=X1;
+sd=R1/d2;
+pc=6*sd
+mprintf("\n Process Capability = %.2f ",pc);
+
+//(d) UNTL (upper natural tolerance limit)
+UNTL = X1 + 3*sd;
+LNTL = X1 - 3*sd;
+USL = 2.07;
+LSL = 2.03;
+z=(USL-X1)/sd
+p=0.8051*100    //The probability from the tables for z
+rw=100-p;
+mprintf("\n Therefore the rework is = %.2f percent]",rw);
+
+
diff --git a/3513/CH6/EX6.1/ch6_1.sce b/3513/CH6/EX6.1/ch6_1.sce
new file mode 100644
index 000000000..858c979a3
--- /dev/null
+++ b/3513/CH6/EX6.1/ch6_1.sce
@@ -0,0 +1,48 @@
+// percentage of rework
+//page no 103
+clear
+clc;
+x=41.283;
+r=0.280;
+k=20;
+ucl=2.07;        //Upper specification limit USL
+lcl=2.03 ;       //Lower specification limit LSL
+A2 = 0.73;
+d2 = 2.059;
+D3 = 0.0;
+D4 = 2.28;
+X1=x/k;
+R1=r/k;
+
+//a) Control limit for R – chart
+ucl1 = D4*R1;
+lcl1 = D3*R1;
+cl1=R1;
+mprintf("\nControl limit for R – chart\n ucl = %.2f ",ucl1);
+mprintf("\nlcl = %.2f ",lcl1);
+mprintf("\nCL = %.2f ",cl1);
+//Control limits for X – chart
+ucl2=X1+A2*R1;
+lcl2=2.05393;
+cl2=X1;
+mprintf("\nControl limit for X – chart\n ucl = %.2f ",ucl2);
+mprintf("\nlcl = %.2f ",lcl2);
+mprintf("\nCL = %.2f ",cl2);
+
+//(b) Process Capability
+cl3=X1;
+sd=R1/d2;
+pc=6*sd
+mprintf("\n Process Capability = %.2f ",pc);
+
+//(d) UNTL (upper natural tolerance limit)
+UNTL = X1 + 3*sd;
+LNTL = X1 - 3*sd;
+USL = 2.07;
+LSL = 2.03;
+z=(USL-X1)/sd
+p=0.8051*100    //The probability from the tables for z
+rw=100-p;
+mprintf("\n Therefore the rework is = %.2f percent]",rw);
+
+
diff --git a/3513/CH6/EX6.2/Ex6_2.sce b/3513/CH6/EX6.2/Ex6_2.sce
new file mode 100644
index 000000000..3c3e7ebab
--- /dev/null
+++ b/3513/CH6/EX6.2/Ex6_2.sce
@@ -0,0 +1,50 @@
+//Compute the control limits
+//page no 105
+clear
+clc;
+K = 25;
+n = 5;
+X1 = 357.5;
+R1 = 8.8;
+USL=14.8;
+LSL= 14.0;
+X2=X1/K;
+R2=R1/K;
+A2 = 0.58;
+d2 = 2.326;
+d3 = 0.0;
+D4 = 2.11;
+//Control limits for R-chart
+UCLR = D4*R2;
+LCL = d3*R2;
+CL = R2;
+//(a) Control limits for X -chart.
+UCL = X2 + A2*R2;
+mprintf("\nControl limit for X – chart\n ucl = %.2f ",UCL);
+LCL = X2 - A2*R2;
+mprintf("\nControl limit for X – chart\n lcl = %.2f ",LCL);
+CL = X2
+mprintf("\nControl limit for X – chart\n cl = %.2f ",CL);
+
+//(b) Since the process is in a state of statistical control
+X21=14.3;
+sd=R2/d2;
+pc=sd*6;
+mprintf("\nProcess capability = %.2f ",pc);
+
+//(c)
+mprintf("\nSince 6σ1 > (USL – LSL), the process is not capable of meeting the specification limits. i.e., 0.907 > 0.8. Rejections are inevitable");
+UNTL = X21 + 3*sd;
+LNTL = X21 - 3*sd;
+CL = X21;
+X=14;
+Z=(X-X21)/sd;
+p = 0.0239*100    //=Probability from tables
+mprintf("\nPercentage of rejection = %.2f ",p);
+//(e) To minimise the percentage of rejection the possible ways are : change the process
+//centre to the specification mean i.e., 14.3 to 14.4. The calculations are shown
+//below.
+X21=14.4;           // To minimise the percentage of rejection
+Z=(X-X21)/sd;
+p = 0.0041*100    //=Probability from tables
+mprintf("\nPercentage of rejection = %.2f ",p);
diff --git a/3513/CH6/EX6.2/ch6_2.sce b/3513/CH6/EX6.2/ch6_2.sce
new file mode 100644
index 000000000..3c3e7ebab
--- /dev/null
+++ b/3513/CH6/EX6.2/ch6_2.sce
@@ -0,0 +1,50 @@
+//Compute the control limits
+//page no 105
+clear
+clc;
+K = 25;
+n = 5;
+X1 = 357.5;
+R1 = 8.8;
+USL=14.8;
+LSL= 14.0;
+X2=X1/K;
+R2=R1/K;
+A2 = 0.58;
+d2 = 2.326;
+d3 = 0.0;
+D4 = 2.11;
+//Control limits for R-chart
+UCLR = D4*R2;
+LCL = d3*R2;
+CL = R2;
+//(a) Control limits for X -chart.
+UCL = X2 + A2*R2;
+mprintf("\nControl limit for X – chart\n ucl = %.2f ",UCL);
+LCL = X2 - A2*R2;
+mprintf("\nControl limit for X – chart\n lcl = %.2f ",LCL);
+CL = X2
+mprintf("\nControl limit for X – chart\n cl = %.2f ",CL);
+
+//(b) Since the process is in a state of statistical control
+X21=14.3;
+sd=R2/d2;
+pc=sd*6;
+mprintf("\nProcess capability = %.2f ",pc);
+
+//(c)
+mprintf("\nSince 6σ1 > (USL – LSL), the process is not capable of meeting the specification limits. i.e., 0.907 > 0.8. Rejections are inevitable");
+UNTL = X21 + 3*sd;
+LNTL = X21 - 3*sd;
+CL = X21;
+X=14;
+Z=(X-X21)/sd;
+p = 0.0239*100    //=Probability from tables
+mprintf("\nPercentage of rejection = %.2f ",p);
+//(e) To minimise the percentage of rejection the possible ways are : change the process
+//centre to the specification mean i.e., 14.3 to 14.4. The calculations are shown
+//below.
+X21=14.4;           // To minimise the percentage of rejection
+Z=(X-X21)/sd;
+p = 0.0041*100    //=Probability from tables
+mprintf("\nPercentage of rejection = %.2f ",p);
diff --git a/3513/CH6/EX6.3/Ex6_3.sce b/3513/CH6/EX6.3/Ex6_3.sce
new file mode 100644
index 000000000..f051af969
--- /dev/null
+++ b/3513/CH6/EX6.3/Ex6_3.sce
@@ -0,0 +1,29 @@
+//Compute percentage of the product will not meet the specification
+//page no 107
+clear
+clc;
+//a
+UCL=129;
+LCL=121;
+X1=125;
+CL=X1;
+USL=135;
+LSL=121;
+A2 = 0.73;
+d2 = 2.059;
+D3 = 0.0;
+D4 = 2.28;
+R2=(UCL-X1)/A2;
+mprintf("\nR2 = %.2f ",R2);
+sd=R2/d2;
+PC=6*sd;
+UNTL = X1 + 3*sd;
+LNTL = X1 + 3*sd;
+Z=(LSL-CL)/sd   //Percentage of rejection
+mprintf("\nPercentage of rejection = %.2f ",Z);
+mprintf("\n a) Probability from table 0.0122 = 1.22 ");
+Z=(LSL-127)/sd   //Percentage of rejection
+mprintf("\nPercentage of rejection = %.2f ",Z);
+mprintf("\n b) Probability from table 0.00135 = 0.135 ");
+TR=0.135*2;
+mprintf("\n c)Since it is symmetric the total percentage of rejection = %.2f ",TR);
diff --git a/3513/CH6/EX6.3/ch6_3.sce b/3513/CH6/EX6.3/ch6_3.sce
new file mode 100644
index 000000000..f051af969
--- /dev/null
+++ b/3513/CH6/EX6.3/ch6_3.sce
@@ -0,0 +1,29 @@
+//Compute percentage of the product will not meet the specification
+//page no 107
+clear
+clc;
+//a
+UCL=129;
+LCL=121;
+X1=125;
+CL=X1;
+USL=135;
+LSL=121;
+A2 = 0.73;
+d2 = 2.059;
+D3 = 0.0;
+D4 = 2.28;
+R2=(UCL-X1)/A2;
+mprintf("\nR2 = %.2f ",R2);
+sd=R2/d2;
+PC=6*sd;
+UNTL = X1 + 3*sd;
+LNTL = X1 + 3*sd;
+Z=(LSL-CL)/sd   //Percentage of rejection
+mprintf("\nPercentage of rejection = %.2f ",Z);
+mprintf("\n a) Probability from table 0.0122 = 1.22 ");
+Z=(LSL-127)/sd   //Percentage of rejection
+mprintf("\nPercentage of rejection = %.2f ",Z);
+mprintf("\n b) Probability from table 0.00135 = 0.135 ");
+TR=0.135*2;
+mprintf("\n c)Since it is symmetric the total percentage of rejection = %.2f ",TR);
diff --git a/3513/CH6/EX6.4/Ex6_4.sce b/3513/CH6/EX6.4/Ex6_4.sce
new file mode 100644
index 000000000..c59eaf4ad
--- /dev/null
+++ b/3513/CH6/EX6.4/Ex6_4.sce
@@ -0,0 +1,47 @@
+//Compute the control limits
+//page no 109
+clear
+clc;
+//(a)
+X1=15350;
+R1 = 411.4;
+SR = 411.4
+K = 25;
+n = 4;
+X2=X1/K;
+R2=R1/K;
+A2 = 0.73;
+d2 = 2.059;
+D3 = 0.0;
+D4 = 2.28;
+USL = X2 + A2*R1;
+LCL = X2 - A2*R1;
+CL=X2;
+//Control limits for R-chart
+UCL1 = D4*R1;
+LCL1 = D3*R1;
+CL1 = R1;
+mprintf("\UCL = %.2f ",UCL1);
+mprintf("\LCL = %.2f ",LCL1);
+mprintf("\CL = %.2f ",CL1);
+//(b) Specification limits are
+USL = 625;
+LCL = 595;
+sd=R1/d2;
+UNTL = X1 + 3*sd;
+LNTL = X1 - 3*sd;
+PC=6*sd;
+mprintf("\Process capability = %.2f ",PC);
+Z=(595-614)/sd;
+mprintf("\nPercentage of rejection = %.2f ",Z);
+mprintf("\n a) Probability from table 0.947 = 91.47 percent ");
+mprintf("\n Rework = 8.53 percent ");
+Z1=(625-618.97)/sd;
+mprintf("\nPercentage of rejection = %.2f ",Z1);
+mprintf("\n  Probability from table 0.7734 = 77.34 percent ");
+mprintf("\n Rework = 22.66 percent ");
+
+
+
+
+
diff --git a/3513/CH6/EX6.4/ch6_4.sce b/3513/CH6/EX6.4/ch6_4.sce
new file mode 100644
index 000000000..c59eaf4ad
--- /dev/null
+++ b/3513/CH6/EX6.4/ch6_4.sce
@@ -0,0 +1,47 @@
+//Compute the control limits
+//page no 109
+clear
+clc;
+//(a)
+X1=15350;
+R1 = 411.4;
+SR = 411.4
+K = 25;
+n = 4;
+X2=X1/K;
+R2=R1/K;
+A2 = 0.73;
+d2 = 2.059;
+D3 = 0.0;
+D4 = 2.28;
+USL = X2 + A2*R1;
+LCL = X2 - A2*R1;
+CL=X2;
+//Control limits for R-chart
+UCL1 = D4*R1;
+LCL1 = D3*R1;
+CL1 = R1;
+mprintf("\UCL = %.2f ",UCL1);
+mprintf("\LCL = %.2f ",LCL1);
+mprintf("\CL = %.2f ",CL1);
+//(b) Specification limits are
+USL = 625;
+LCL = 595;
+sd=R1/d2;
+UNTL = X1 + 3*sd;
+LNTL = X1 - 3*sd;
+PC=6*sd;
+mprintf("\Process capability = %.2f ",PC);
+Z=(595-614)/sd;
+mprintf("\nPercentage of rejection = %.2f ",Z);
+mprintf("\n a) Probability from table 0.947 = 91.47 percent ");
+mprintf("\n Rework = 8.53 percent ");
+Z1=(625-618.97)/sd;
+mprintf("\nPercentage of rejection = %.2f ",Z1);
+mprintf("\n  Probability from table 0.7734 = 77.34 percent ");
+mprintf("\n Rework = 22.66 percent ");
+
+
+
+
+
diff --git a/3513/CH6/EX6.5/Ex6_5.sce b/3513/CH6/EX6.5/Ex6_5.sce
new file mode 100644
index 000000000..a860ce32f
--- /dev/null
+++ b/3513/CH6/EX6.5/Ex6_5.sce
@@ -0,0 +1,31 @@
+//Compute process average and std. deviation
+//page no 112
+clear
+clc;
+
+USL = 51;
+LSL = 45;
+UNTL=USL;
+LNTL=LSL;
+sd=1;
+X2=48;
+mprintf("\The process average = %.2f \n",X2);
+A2 = 0.73;
+d2 = 2.059;
+D3 = 0.0;
+D4 = 2.28;
+R2= sd*d2;
+//Control limits for R-chart
+UCL = D4*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = D3*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R2;
+mprintf("\CL = %.2f \n",CL);
+//Control limits for X -chart
+UCL = X2 + A2*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X2 - A2*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL=X2;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH6/EX6.5/ch6_5.sce b/3513/CH6/EX6.5/ch6_5.sce
new file mode 100644
index 000000000..a860ce32f
--- /dev/null
+++ b/3513/CH6/EX6.5/ch6_5.sce
@@ -0,0 +1,31 @@
+//Compute process average and std. deviation
+//page no 112
+clear
+clc;
+
+USL = 51;
+LSL = 45;
+UNTL=USL;
+LNTL=LSL;
+sd=1;
+X2=48;
+mprintf("\The process average = %.2f \n",X2);
+A2 = 0.73;
+d2 = 2.059;
+D3 = 0.0;
+D4 = 2.28;
+R2= sd*d2;
+//Control limits for R-chart
+UCL = D4*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = D3*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R2;
+mprintf("\CL = %.2f \n",CL);
+//Control limits for X -chart
+UCL = X2 + A2*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X2 - A2*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL=X2;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH6/EX6.6/Ex6_6.sce b/3513/CH6/EX6.6/Ex6_6.sce
new file mode 100644
index 000000000..673ea2dc7
--- /dev/null
+++ b/3513/CH6/EX6.6/Ex6_6.sce
@@ -0,0 +1,45 @@
+//Compute the 3σ limits
+//page no 113
+clear
+clc;
+USL = 125;
+LSL = 115;
+X1 = 120;
+sd = 1.5;
+X2=X1;
+n = 4;
+A2 = 0.73;
+d2 = 2.059;
+D3 = 0.0;
+D4 = 2.28;
+R1=sd*d2;
+n=4;
+//Control limits for R-chart
+UCL = D4*R1;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = D3*R1;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R1;
+mprintf("\CL = %.2f \n",CL);
+//Control limits for X -chart
+UCL = X1+A2*R1;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X1-A2*R1;
+mprintf("\LCL = %.2f \n",LCL);
+CL = X2;
+mprintf("\CL = %.2f \n",CL);
+Xnew=122.4;
+Xnew2=117.6;
+Zb=(Xnew/Xnew2)/sd;
+pr = 0.0418; //probability
+Za =(125-122.4)/sd;
+pr = 95.82; //probability
+mprintf("\probability = %.2f \n",pr);
+sdn=sd/sqrt(n);
+Za =(117.74-117.6)/sdn;
+pr = 57.14; //probability
+Zb =(122.26-122.4)/sdn;
+pr = 42.86; //probability
+mprintf("\probability = %.2f \n",pr);
+P=100-pr;
+mprintf("\final probability = %.2f \n",P);
diff --git a/3513/CH6/EX6.6/ch6_6.sce b/3513/CH6/EX6.6/ch6_6.sce
new file mode 100644
index 000000000..673ea2dc7
--- /dev/null
+++ b/3513/CH6/EX6.6/ch6_6.sce
@@ -0,0 +1,45 @@
+//Compute the 3σ limits
+//page no 113
+clear
+clc;
+USL = 125;
+LSL = 115;
+X1 = 120;
+sd = 1.5;
+X2=X1;
+n = 4;
+A2 = 0.73;
+d2 = 2.059;
+D3 = 0.0;
+D4 = 2.28;
+R1=sd*d2;
+n=4;
+//Control limits for R-chart
+UCL = D4*R1;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = D3*R1;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R1;
+mprintf("\CL = %.2f \n",CL);
+//Control limits for X -chart
+UCL = X1+A2*R1;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X1-A2*R1;
+mprintf("\LCL = %.2f \n",LCL);
+CL = X2;
+mprintf("\CL = %.2f \n",CL);
+Xnew=122.4;
+Xnew2=117.6;
+Zb=(Xnew/Xnew2)/sd;
+pr = 0.0418; //probability
+Za =(125-122.4)/sd;
+pr = 95.82; //probability
+mprintf("\probability = %.2f \n",pr);
+sdn=sd/sqrt(n);
+Za =(117.74-117.6)/sdn;
+pr = 57.14; //probability
+Zb =(122.26-122.4)/sdn;
+pr = 42.86; //probability
+mprintf("\probability = %.2f \n",pr);
+P=100-pr;
+mprintf("\final probability = %.2f \n",P);
diff --git a/3513/CH6/EX6.7/Ex6_7.sce b/3513/CH6/EX6.7/Ex6_7.sce
new file mode 100644
index 000000000..17932af52
--- /dev/null
+++ b/3513/CH6/EX6.7/Ex6_7.sce
@@ -0,0 +1,48 @@
+//Compute rejections percentage
+//page no 117
+clear
+clc;
+n = 2;
+K = 2;
+R1 = 0.81;
+X = 27.635;
+USL = 1.207;
+LSL = 1.033;
+d2 = 1.128;
+A2 = 1.88;
+D3 = 0.0;
+D4 = 3.27;
+X2=X/K;
+X2=R1/K;
+sd=R1/d2
+mprintf("\(a) Standard Deviation = %.2f \n",sd);
+UNTL = X2+3*sd;
+LNTL = X2-3*sd;
+mprintf("\UNTL = %.2f \n",UNTL);
+mprintf("\LNTL = %.2f \n",LNTL);
+CL = X2;
+mprintf("\CL = %.2f \n",CL);
+//b) The percentage of rejection
+Z =(1.033-1.1054)/sd;
+mprintf("\Z = %.2f \n",Z);
+pr = 0.59; //probability
+mprintf("\n Hence probability = %.2f \n",pr);
+T=USL-LSL;
+T2=6*sd;
+mprintf("\n (c) 0.1722 (T2) < 0.174(T) the process is capable of meeting the specification limits.");
+//(d)
+X2=1.12;
+//Control limits for R-chart
+UCL = D4*R1;
+mprintf("\ (b) UCL = %.2f \n",UCL);
+LCL = D3*R1;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R1;
+mprintf("\CL = %.2f \n",CL);
+//Control limits for X -chart
+UCL = X2+A2*R1;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X2-A2*R1;
+mprintf("\LCL = %.2f \n",LCL);
+CL = X2;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH6/EX6.7/ch6_7.sce b/3513/CH6/EX6.7/ch6_7.sce
new file mode 100644
index 000000000..17932af52
--- /dev/null
+++ b/3513/CH6/EX6.7/ch6_7.sce
@@ -0,0 +1,48 @@
+//Compute rejections percentage
+//page no 117
+clear
+clc;
+n = 2;
+K = 2;
+R1 = 0.81;
+X = 27.635;
+USL = 1.207;
+LSL = 1.033;
+d2 = 1.128;
+A2 = 1.88;
+D3 = 0.0;
+D4 = 3.27;
+X2=X/K;
+X2=R1/K;
+sd=R1/d2
+mprintf("\(a) Standard Deviation = %.2f \n",sd);
+UNTL = X2+3*sd;
+LNTL = X2-3*sd;
+mprintf("\UNTL = %.2f \n",UNTL);
+mprintf("\LNTL = %.2f \n",LNTL);
+CL = X2;
+mprintf("\CL = %.2f \n",CL);
+//b) The percentage of rejection
+Z =(1.033-1.1054)/sd;
+mprintf("\Z = %.2f \n",Z);
+pr = 0.59; //probability
+mprintf("\n Hence probability = %.2f \n",pr);
+T=USL-LSL;
+T2=6*sd;
+mprintf("\n (c) 0.1722 (T2) < 0.174(T) the process is capable of meeting the specification limits.");
+//(d)
+X2=1.12;
+//Control limits for R-chart
+UCL = D4*R1;
+mprintf("\ (b) UCL = %.2f \n",UCL);
+LCL = D3*R1;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R1;
+mprintf("\CL = %.2f \n",CL);
+//Control limits for X -chart
+UCL = X2+A2*R1;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X2-A2*R1;
+mprintf("\LCL = %.2f \n",LCL);
+CL = X2;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH6/EX6.8/Ex6_8.sce b/3513/CH6/EX6.8/Ex6_8.sce
new file mode 100644
index 000000000..35fc08ec3
--- /dev/null
+++ b/3513/CH6/EX6.8/Ex6_8.sce
@@ -0,0 +1,58 @@
+//Determine the control limits for X -R chart.
+//page no 119
+clear
+clc;
+USL = 175;
+LSL = 135;
+n = 5;
+K = 50;
+X1 = 7660;
+R1=880;
+X2=X1/K;
+R2=R1/K;
+A2 = 0.58;
+d2 = 2.326;
+D3 = 0.0;
+D4 = 2.11;
+//Control limits for R-chart
+UCL = D4*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = D3*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R1;
+mprintf("\CL = %.2f \n",CL);
+//Control limits for X -chart
+UCL = X2+A2*R1;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X2-A2*R1;
+mprintf("\LCL = %.2f \n",LCL);
+CL = X2;
+mprintf("\CL = %.2f \n",CL);
+X2 = 153.2;
+sd=R2/d2;
+UNTL = X2+3*sd;
+LNTL = X2-3*sd;
+mprintf("\UNTL = %.2f \n",UNTL);
+mprintf("\LNTL = %.2f \n",LNTL);
+USL = 175;
+LSL = 135;
+Z=(LSL-CL)/sd;
+pr1 = 0.81; //probability
+mprintf("\n Hence probability = %.2f \n",pr1);
+Z=(USL-CL)/sd;
+pr2 = 99.8; //probability
+mprintf("\n Hence probability = %.2f \n",pr2);
+R=100-pr2;
+R2=pr1+R;
+mprintf("\n Total portion of the product which does not confirm with the specifications = %.2f \n",R2);
+//(c)
+Cl=155;
+Z=(LSL-CL)/sd;
+pr1 = 0.41; //probability
+mprintf("\n Hence probability = %.2f \n",pr1);
+Z=(USL-CL)/sd;
+pr2 = 99.59; //probability
+mprintf("\n Hence probability = %.2f \n",pr2);
+R=100-pr2;
+R2=pr1+R;
+mprintf("\n The total rejections will be = %.2f \n",R2);
diff --git a/3513/CH6/EX6.8/ch6_8.sce b/3513/CH6/EX6.8/ch6_8.sce
new file mode 100644
index 000000000..35fc08ec3
--- /dev/null
+++ b/3513/CH6/EX6.8/ch6_8.sce
@@ -0,0 +1,58 @@
+//Determine the control limits for X -R chart.
+//page no 119
+clear
+clc;
+USL = 175;
+LSL = 135;
+n = 5;
+K = 50;
+X1 = 7660;
+R1=880;
+X2=X1/K;
+R2=R1/K;
+A2 = 0.58;
+d2 = 2.326;
+D3 = 0.0;
+D4 = 2.11;
+//Control limits for R-chart
+UCL = D4*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = D3*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R1;
+mprintf("\CL = %.2f \n",CL);
+//Control limits for X -chart
+UCL = X2+A2*R1;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X2-A2*R1;
+mprintf("\LCL = %.2f \n",LCL);
+CL = X2;
+mprintf("\CL = %.2f \n",CL);
+X2 = 153.2;
+sd=R2/d2;
+UNTL = X2+3*sd;
+LNTL = X2-3*sd;
+mprintf("\UNTL = %.2f \n",UNTL);
+mprintf("\LNTL = %.2f \n",LNTL);
+USL = 175;
+LSL = 135;
+Z=(LSL-CL)/sd;
+pr1 = 0.81; //probability
+mprintf("\n Hence probability = %.2f \n",pr1);
+Z=(USL-CL)/sd;
+pr2 = 99.8; //probability
+mprintf("\n Hence probability = %.2f \n",pr2);
+R=100-pr2;
+R2=pr1+R;
+mprintf("\n Total portion of the product which does not confirm with the specifications = %.2f \n",R2);
+//(c)
+Cl=155;
+Z=(LSL-CL)/sd;
+pr1 = 0.41; //probability
+mprintf("\n Hence probability = %.2f \n",pr1);
+Z=(USL-CL)/sd;
+pr2 = 99.59; //probability
+mprintf("\n Hence probability = %.2f \n",pr2);
+R=100-pr2;
+R2=pr1+R;
+mprintf("\n The total rejections will be = %.2f \n",R2);
diff --git a/3513/CH6/EX6.9/Ex6_9.sce b/3513/CH6/EX6.9/Ex6_9.sce
new file mode 100644
index 000000000..88613e0f2
--- /dev/null
+++ b/3513/CH6/EX6.9/Ex6_9.sce
@@ -0,0 +1,66 @@
+//Determine the control limits for X and R chart..
+//page no 121
+clear
+clc;
+
+K = 20;
+X1 = 669.2;
+R1=126.0;
+X2=X1/K;
+R2=R1/K;
+A2 = 0.58;
+d2 = 2.326;
+D3 = 0.0;
+D4 = 2.11;
+//(a) Control limits for R-chart
+mprintf("Control limits for R-chart");
+UCL = D4*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = D3*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R1;
+mprintf("\CL = %.2f \n",CL);
+R2=(R1-14-19)/(K-2);
+//Again Control limits for R-chart
+UCL = D4*R2;
+mprintf("Again Control limits for R-chart");
+mprintf("\UCL = %.2f \n",UCL);
+LCL = D3*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R1;
+mprintf("\CL = %.2f \n",CL);
+R2=(R1-14-19-13)/(K-3);
+//Control limits for X -chart
+mprintf("Control limits for X-chart");
+UCL = X2+A2*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X2-A2*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = X2;
+X2=33.22;
+mprintf("Again Control limits for X-chart");
+UCL = X2+A2*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X2-A2*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = X2;
+//(b)
+sd=R2/d2;
+mprintf("\Standard Deviation = %.2f \n",sd);
+pc=6*sd;
+mprintf("\process capability = %.2f \n",pc);
+UNTL = X2+3*sd;
+LNTL = X2-3*sd;
+mprintf("\UNTL = %.2f \n",UNTL);
+mprintf("\LNTL = %.2f \n",LNTL);
+USL = 38;
+LSL = 28;
+Z=(LSL-CL)/sd;
+pr1 = 0.52; //probability
+mprintf("\n Hence probability = %.2f \n",pr1);
+Z=(USL-CL)/sd;
+pr2 = 99.09; //probability
+mprintf("\n Hence probability = %.2f \n",pr2);
+R=100-pr2;
+R2=pr1+R;
+mprintf("\n The total rejections will be = %.2f \n",R2);
diff --git a/3513/CH6/EX6.9/ch6_9.sce b/3513/CH6/EX6.9/ch6_9.sce
new file mode 100644
index 000000000..88613e0f2
--- /dev/null
+++ b/3513/CH6/EX6.9/ch6_9.sce
@@ -0,0 +1,66 @@
+//Determine the control limits for X and R chart..
+//page no 121
+clear
+clc;
+
+K = 20;
+X1 = 669.2;
+R1=126.0;
+X2=X1/K;
+R2=R1/K;
+A2 = 0.58;
+d2 = 2.326;
+D3 = 0.0;
+D4 = 2.11;
+//(a) Control limits for R-chart
+mprintf("Control limits for R-chart");
+UCL = D4*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = D3*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R1;
+mprintf("\CL = %.2f \n",CL);
+R2=(R1-14-19)/(K-2);
+//Again Control limits for R-chart
+UCL = D4*R2;
+mprintf("Again Control limits for R-chart");
+mprintf("\UCL = %.2f \n",UCL);
+LCL = D3*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = R1;
+mprintf("\CL = %.2f \n",CL);
+R2=(R1-14-19-13)/(K-3);
+//Control limits for X -chart
+mprintf("Control limits for X-chart");
+UCL = X2+A2*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X2-A2*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = X2;
+X2=33.22;
+mprintf("Again Control limits for X-chart");
+UCL = X2+A2*R2;
+mprintf("\UCL = %.2f \n",UCL);
+LCL = X2-A2*R2;
+mprintf("\LCL = %.2f \n",LCL);
+CL = X2;
+//(b)
+sd=R2/d2;
+mprintf("\Standard Deviation = %.2f \n",sd);
+pc=6*sd;
+mprintf("\process capability = %.2f \n",pc);
+UNTL = X2+3*sd;
+LNTL = X2-3*sd;
+mprintf("\UNTL = %.2f \n",UNTL);
+mprintf("\LNTL = %.2f \n",LNTL);
+USL = 38;
+LSL = 28;
+Z=(LSL-CL)/sd;
+pr1 = 0.52; //probability
+mprintf("\n Hence probability = %.2f \n",pr1);
+Z=(USL-CL)/sd;
+pr2 = 99.09; //probability
+mprintf("\n Hence probability = %.2f \n",pr2);
+R=100-pr2;
+R2=pr1+R;
+mprintf("\n The total rejections will be = %.2f \n",R2);