diff options
Diffstat (limited to '3513')
68 files changed, 1852 insertions, 0 deletions
diff --git a/3513/CH2/EX2.1/Ex2_1.sce b/3513/CH2/EX2.1/Ex2_1.sce new file mode 100644 index 000000000..cf0b1b4d7 --- /dev/null +++ b/3513/CH2/EX2.1/Ex2_1.sce @@ -0,0 +1,24 @@ +// To determine percentage of specimens
+//page no 21
+clear
+clc;
+m=5600;//mean
+sd=840;//psi
+lwrlmt=5600;
+uprlmt=6200;
+z=-0.7142;
+p=23.89;//probability
+sp=(100-p)-p;
+mprintf("Therefore the percentage of the specimens falling between 5000 and 6200 = %.2f percentage",sp);
+//(b) Percentage of specimens falling below 4000
+z=4000-(lwrlmt/sd);
+p=0.0287; //The probability from tables
+p=p*100;
+pa=100-p
+mprintf("\nTherefore the probability above 4000 = %.2f percentage",pa);
+//(c)
+z=3500-(lwrlmt/sd);
+p=0.0062;//The probability from tables
+pa=p*100;
+mprintf("\nTherefore the percentage of specimen falling below 3500 = %.2f percentage",pa);
+
diff --git a/3513/CH2/EX2.1/ch2_1.sce b/3513/CH2/EX2.1/ch2_1.sce new file mode 100644 index 000000000..cf0b1b4d7 --- /dev/null +++ b/3513/CH2/EX2.1/ch2_1.sce @@ -0,0 +1,24 @@ +// To determine percentage of specimens
+//page no 21
+clear
+clc;
+m=5600;//mean
+sd=840;//psi
+lwrlmt=5600;
+uprlmt=6200;
+z=-0.7142;
+p=23.89;//probability
+sp=(100-p)-p;
+mprintf("Therefore the percentage of the specimens falling between 5000 and 6200 = %.2f percentage",sp);
+//(b) Percentage of specimens falling below 4000
+z=4000-(lwrlmt/sd);
+p=0.0287; //The probability from tables
+p=p*100;
+pa=100-p
+mprintf("\nTherefore the probability above 4000 = %.2f percentage",pa);
+//(c)
+z=3500-(lwrlmt/sd);
+p=0.0062;//The probability from tables
+pa=p*100;
+mprintf("\nTherefore the percentage of specimen falling below 3500 = %.2f percentage",pa);
+
diff --git a/3513/CH2/EX2.2/Ex2_2.sce b/3513/CH2/EX2.2/Ex2_2.sce new file mode 100644 index 000000000..32efc5098 --- /dev/null +++ b/3513/CH2/EX2.2/Ex2_2.sce @@ -0,0 +1,21 @@ +// To determine percentage of specimens
+//page no 22
+clear
+clc;
+m=118.5;// volts mean
+sd=1.2;// volts psi
+p1=0.0188;//The probability from tables
+p1=p1*100;
+x=116;
+z=x-(m/sd);
+p2=0.8944;//The probability from tables
+p2=p2*100;
+p=p2-p1;
+mprintf("Therefore the percentage of specimen falling between 116 and 120 volts = %.2f percentage",p);
+//(b)
+p=1.2; //The probability from tables
+x=115;
+z=1.404;
+m=115+z;
+av=m-z
+mprintf("\nThe adjustment is 116.404 – 115 = %.2f Voltage",av);
diff --git a/3513/CH2/EX2.2/ch2_2.sce b/3513/CH2/EX2.2/ch2_2.sce new file mode 100644 index 000000000..32efc5098 --- /dev/null +++ b/3513/CH2/EX2.2/ch2_2.sce @@ -0,0 +1,21 @@ +// To determine percentage of specimens
+//page no 22
+clear
+clc;
+m=118.5;// volts mean
+sd=1.2;// volts psi
+p1=0.0188;//The probability from tables
+p1=p1*100;
+x=116;
+z=x-(m/sd);
+p2=0.8944;//The probability from tables
+p2=p2*100;
+p=p2-p1;
+mprintf("Therefore the percentage of specimen falling between 116 and 120 volts = %.2f percentage",p);
+//(b)
+p=1.2; //The probability from tables
+x=115;
+z=1.404;
+m=115+z;
+av=m-z
+mprintf("\nThe adjustment is 116.404 – 115 = %.2f Voltage",av);
diff --git a/3513/CH2/EX2.3/Ex2_3.sce b/3513/CH2/EX2.3/Ex2_3.sce new file mode 100644 index 000000000..d683e1edc --- /dev/null +++ b/3513/CH2/EX2.3/Ex2_3.sce @@ -0,0 +1,15 @@ +// To determine percentage of specimens
+//page no 22
+clear
+clc;
+m=10;// cm mean
+sd=0.03;// standard deviation
+x=10.075;
+z=0.9938 //The values of Z from normal tables
+z=z*100;
+t=100-z;
+mprintf("Therefore 100-99.38 = %.2f percent of items are falling greater than 10.075 cms.",t);
+//(b)
+z=1.04;
+x=(z*sd)-m;
+mprintf("\nBelow %.2f value of diameter will 15 percentage piston rings fall",x);
diff --git a/3513/CH2/EX2.3/ch2_3.sce b/3513/CH2/EX2.3/ch2_3.sce new file mode 100644 index 000000000..d683e1edc --- /dev/null +++ b/3513/CH2/EX2.3/ch2_3.sce @@ -0,0 +1,15 @@ +// To determine percentage of specimens
+//page no 22
+clear
+clc;
+m=10;// cm mean
+sd=0.03;// standard deviation
+x=10.075;
+z=0.9938 //The values of Z from normal tables
+z=z*100;
+t=100-z;
+mprintf("Therefore 100-99.38 = %.2f percent of items are falling greater than 10.075 cms.",t);
+//(b)
+z=1.04;
+x=(z*sd)-m;
+mprintf("\nBelow %.2f value of diameter will 15 percentage piston rings fall",x);
diff --git a/3513/CH2/EX2.4/Ex2_4.sce b/3513/CH2/EX2.4/Ex2_4.sce new file mode 100644 index 000000000..22978af9c --- /dev/null +++ b/3513/CH2/EX2.4/Ex2_4.sce @@ -0,0 +1,9 @@ +// To determine percentage of specimens
+//page no 22
+clear
+clc;
+sd=0.01;// standard deviation
+x=0.567;
+z=-2.16//The values of Z from normal tables
+m=(z*sd)-x;
+mprintf("Therefore the mean weight required is %.3f.",abs(m));
diff --git a/3513/CH2/EX2.4/ch2_4.sce b/3513/CH2/EX2.4/ch2_4.sce new file mode 100644 index 000000000..22978af9c --- /dev/null +++ b/3513/CH2/EX2.4/ch2_4.sce @@ -0,0 +1,9 @@ +// To determine percentage of specimens
+//page no 22
+clear
+clc;
+sd=0.01;// standard deviation
+x=0.567;
+z=-2.16//The values of Z from normal tables
+m=(z*sd)-x;
+mprintf("Therefore the mean weight required is %.3f.",abs(m));
diff --git a/3513/CH2/EX2.5/Ex2_5.sce b/3513/CH2/EX2.5/Ex2_5.sce new file mode 100644 index 000000000..bc310e63c --- /dev/null +++ b/3513/CH2/EX2.5/Ex2_5.sce @@ -0,0 +1,9 @@ +// To determine percentage of specimens
+//page no 23
+clear
+clc;
+m=10;//years mean
+sd=2;// years standard deviation
+z=1.64;//The values of Z from normal tables
+x=(z*sd)-m;
+mprintf("Therefore the value of x =%.3f.years",abs(x));
diff --git a/3513/CH2/EX2.5/ch2_5.sce b/3513/CH2/EX2.5/ch2_5.sce new file mode 100644 index 000000000..bc310e63c --- /dev/null +++ b/3513/CH2/EX2.5/ch2_5.sce @@ -0,0 +1,9 @@ +// To determine percentage of specimens
+//page no 23
+clear
+clc;
+m=10;//years mean
+sd=2;// years standard deviation
+z=1.64;//The values of Z from normal tables
+x=(z*sd)-m;
+mprintf("Therefore the value of x =%.3f.years",abs(x));
diff --git a/3513/CH2/EX2.6/Ex2_6.sce b/3513/CH2/EX2.6/Ex2_6.sce new file mode 100644 index 000000000..059ea33a3 --- /dev/null +++ b/3513/CH2/EX2.6/Ex2_6.sce @@ -0,0 +1,27 @@ +// To determine probabilty
+//page no 24
+clear
+clc;
+n=50;
+p=0.25;
+q=0.75;
+sd=3.062; // standard deviation
+x1=14.5;
+x2=15.5;
+x3=11.5
+p1=0.7422*100 // probability from table
+p2=0.8340*100 // probability from table
+//(a) Exactly 15 defective
+m=n*p
+z1=((x1-m)/sd)*100
+z2=((x2-m)/sd)*100
+p=(p2-p1)
+mprintf("\nTherefore the probability Exactly 15 defective = %.2f percentage",p);
+//(b) Less than or equal to 11 defective
+z3=((x3-m)/sd)*100
+p3=37.45
+mprintf("\nTherefore the probabilityLess than or equal to 11 defective = %.2f percentage",p3);
+//(c) More than 15 or 15 defectives
+z4=((x1-m)/sd)*100
+p4=100-p1
+mprintf("\nTherefore the probability More than 15 or 15 defectives = %.2f percentage",p4);
diff --git a/3513/CH2/EX2.6/ch2_6.sce b/3513/CH2/EX2.6/ch2_6.sce new file mode 100644 index 000000000..059ea33a3 --- /dev/null +++ b/3513/CH2/EX2.6/ch2_6.sce @@ -0,0 +1,27 @@ +// To determine probabilty
+//page no 24
+clear
+clc;
+n=50;
+p=0.25;
+q=0.75;
+sd=3.062; // standard deviation
+x1=14.5;
+x2=15.5;
+x3=11.5
+p1=0.7422*100 // probability from table
+p2=0.8340*100 // probability from table
+//(a) Exactly 15 defective
+m=n*p
+z1=((x1-m)/sd)*100
+z2=((x2-m)/sd)*100
+p=(p2-p1)
+mprintf("\nTherefore the probability Exactly 15 defective = %.2f percentage",p);
+//(b) Less than or equal to 11 defective
+z3=((x3-m)/sd)*100
+p3=37.45
+mprintf("\nTherefore the probabilityLess than or equal to 11 defective = %.2f percentage",p3);
+//(c) More than 15 or 15 defectives
+z4=((x1-m)/sd)*100
+p4=100-p1
+mprintf("\nTherefore the probability More than 15 or 15 defectives = %.2f percentage",p4);
diff --git a/3513/CH2/EX2.7/Ex2_7.sce b/3513/CH2/EX2.7/Ex2_7.sce new file mode 100644 index 000000000..e15d32448 --- /dev/null +++ b/3513/CH2/EX2.7/Ex2_7.sce @@ -0,0 +1,27 @@ +// To determine probabilty
+//page no 25
+clear
+clc;
+sd=3.1622; // standard deviation
+x1=7.5;
+x2=8.5;
+x3=13.5;
+x3=6.5;
+m=10;
+p1=0.2148*100; // probability from table
+p2=0.3192*100; // probability from table
+p3=0.8643*100; // probability from table
+p4=0.1357*100; // probability from table
+
+//(a) Exactly 8 defectives
+z1=((x1-m)/sd)*100
+z2=((x2-m)/sd)*100
+pa=(p2-p1)
+mprintf("\nTherefore for exactly 8 defectives = %.2f percentage",pa);
+//(b) For 14 or more defectives
+z3=((x3-m)/sd)*100
+pa2=100-p3
+mprintf("\nFor 14 or more is = %.2f percentage",pa2);
+//(c) For less than 6 defectives
+z4=((x1-m)/sd)*100
+mprintf("\nFor less than 6 defectives = %.2f percentage",p4);
diff --git a/3513/CH2/EX2.7/ch2_7.sce b/3513/CH2/EX2.7/ch2_7.sce new file mode 100644 index 000000000..e15d32448 --- /dev/null +++ b/3513/CH2/EX2.7/ch2_7.sce @@ -0,0 +1,27 @@ +// To determine probabilty
+//page no 25
+clear
+clc;
+sd=3.1622; // standard deviation
+x1=7.5;
+x2=8.5;
+x3=13.5;
+x3=6.5;
+m=10;
+p1=0.2148*100; // probability from table
+p2=0.3192*100; // probability from table
+p3=0.8643*100; // probability from table
+p4=0.1357*100; // probability from table
+
+//(a) Exactly 8 defectives
+z1=((x1-m)/sd)*100
+z2=((x2-m)/sd)*100
+pa=(p2-p1)
+mprintf("\nTherefore for exactly 8 defectives = %.2f percentage",pa);
+//(b) For 14 or more defectives
+z3=((x3-m)/sd)*100
+pa2=100-p3
+mprintf("\nFor 14 or more is = %.2f percentage",pa2);
+//(c) For less than 6 defectives
+z4=((x1-m)/sd)*100
+mprintf("\nFor less than 6 defectives = %.2f percentage",p4);
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);
diff --git a/3513/CH7/EX7.1/Ex7_1.sce b/3513/CH7/EX7.1/Ex7_1.sce new file mode 100644 index 000000000..c5df44549 --- /dev/null +++ b/3513/CH7/EX7.1/Ex7_1.sce @@ -0,0 +1,30 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 139
+clear
+clc;
+d=116;
+n=1000;
+P1=d/n;
+//Control limits for p-chart
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+P1 =(116-15-18)/(18*50);
+//Again control limits
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+//Again one more subgroup is crossing the crossing limits
+P1 =(116-15-18-12)/(17*50);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.1/ch7_1.sce b/3513/CH7/EX7.1/ch7_1.sce new file mode 100644 index 000000000..c5df44549 --- /dev/null +++ b/3513/CH7/EX7.1/ch7_1.sce @@ -0,0 +1,30 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 139
+clear
+clc;
+d=116;
+n=1000;
+P1=d/n;
+//Control limits for p-chart
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+P1 =(116-15-18)/(18*50);
+//Again control limits
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+//Again one more subgroup is crossing the crossing limits
+P1 =(116-15-18-12)/(17*50);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.10/Ex7_10.sce b/3513/CH7/EX7.10/Ex7_10.sce new file mode 100644 index 000000000..cfd8c5444 --- /dev/null +++ b/3513/CH7/EX7.10/Ex7_10.sce @@ -0,0 +1,39 @@ +//Calculate value of standard fraction defective would you recommend for the future
+//page no 153
+clear
+clc;
+P1=0.04;
+n=1600;
+minsg=P1+(25/100)*n;
+maxsg=P1-(25/100)*n;
+
+//(a) Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+//The individual control limits for the 3rd day
+//Control limits for p-chart
+mprintf("The individual control limits for the 3rd day");
+n=900;
+mprintf("\Again for \nP1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.4f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+//The individual control limits for the 6th day
+//Control limits for p-chart
+mprintf("The individual control limits for the 6th day");
+n=2000;
+mprintf("\Again for \nP1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.4f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.10/ch7_10.sce b/3513/CH7/EX7.10/ch7_10.sce new file mode 100644 index 000000000..cfd8c5444 --- /dev/null +++ b/3513/CH7/EX7.10/ch7_10.sce @@ -0,0 +1,39 @@ +//Calculate value of standard fraction defective would you recommend for the future
+//page no 153
+clear
+clc;
+P1=0.04;
+n=1600;
+minsg=P1+(25/100)*n;
+maxsg=P1-(25/100)*n;
+
+//(a) Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+//The individual control limits for the 3rd day
+//Control limits for p-chart
+mprintf("The individual control limits for the 3rd day");
+n=900;
+mprintf("\Again for \nP1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.4f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+//The individual control limits for the 6th day
+//Control limits for p-chart
+mprintf("The individual control limits for the 6th day");
+n=2000;
+mprintf("\Again for \nP1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.4f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.2/Ex7_2.sce b/3513/CH7/EX7.2/Ex7_2.sce new file mode 100644 index 000000000..aca5acedc --- /dev/null +++ b/3513/CH7/EX7.2/Ex7_2.sce @@ -0,0 +1,31 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 143
+clear
+clc;
+d=20;
+n=4000;
+P1=0.116;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+P1 =(116-15-18)/(18*50);
+//Again control limits
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+//Again one more subgroup is crossing the crossing limits
+P1 =(116-15-18-12)/(17*50);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.2/ch7_2.sce b/3513/CH7/EX7.2/ch7_2.sce new file mode 100644 index 000000000..aca5acedc --- /dev/null +++ b/3513/CH7/EX7.2/ch7_2.sce @@ -0,0 +1,31 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 143
+clear
+clc;
+d=20;
+n=4000;
+P1=0.116;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+P1 =(116-15-18)/(18*50);
+//Again control limits
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+//Again one more subgroup is crossing the crossing limits
+P1 =(116-15-18-12)/(17*50);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.3/Ex7_3.sce b/3513/CH7/EX7.3/Ex7_3.sce new file mode 100644 index 000000000..bf29035e1 --- /dev/null +++ b/3513/CH7/EX7.3/Ex7_3.sce @@ -0,0 +1,23 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 144
+clear
+clc;
+d=312;
+n=65;
+P1=d/(n*15);
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+P1 =(312-6)/(65*14);
+//Again control limits
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.3/ch7_3.sce b/3513/CH7/EX7.3/ch7_3.sce new file mode 100644 index 000000000..bf29035e1 --- /dev/null +++ b/3513/CH7/EX7.3/ch7_3.sce @@ -0,0 +1,23 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 144
+clear
+clc;
+d=312;
+n=65;
+P1=d/(n*15);
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+P1 =(312-6)/(65*14);
+//Again control limits
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.4/Ex7_4.sce b/3513/CH7/EX7.4/Ex7_4.sce new file mode 100644 index 000000000..682c71ec6 --- /dev/null +++ b/3513/CH7/EX7.4/Ex7_4.sce @@ -0,0 +1,16 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 147
+clear
+clc;
+d=26;
+ns=447;
+P1=d/ns;
+n=48;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.4/ch7_4.sce b/3513/CH7/EX7.4/ch7_4.sce new file mode 100644 index 000000000..682c71ec6 --- /dev/null +++ b/3513/CH7/EX7.4/ch7_4.sce @@ -0,0 +1,16 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 147
+clear
+clc;
+d=26;
+ns=447;
+P1=d/ns;
+n=48;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.5/Ex7_5.sce b/3513/CH7/EX7.5/Ex7_5.sce new file mode 100644 index 000000000..e693a650e --- /dev/null +++ b/3513/CH7/EX7.5/Ex7_5.sce @@ -0,0 +1,16 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 148
+clear
+clc;
+ds=2349;
+ns=19720;
+P1=ds/ns;
+n=960;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.5/ch7_5.sce b/3513/CH7/EX7.5/ch7_5.sce new file mode 100644 index 000000000..e693a650e --- /dev/null +++ b/3513/CH7/EX7.5/ch7_5.sce @@ -0,0 +1,16 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 148
+clear
+clc;
+ds=2349;
+ns=19720;
+P1=ds/ns;
+n=960;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.6/Ex7_6.sce b/3513/CH7/EX7.6/Ex7_6.sce new file mode 100644 index 000000000..40673fb1d --- /dev/null +++ b/3513/CH7/EX7.6/Ex7_6.sce @@ -0,0 +1,16 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 148
+clear
+clc;
+ds=614;
+ns=33725;
+P1=ds/ns;
+n=6000;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.6/ch7_6.sce b/3513/CH7/EX7.6/ch7_6.sce new file mode 100644 index 000000000..40673fb1d --- /dev/null +++ b/3513/CH7/EX7.6/ch7_6.sce @@ -0,0 +1,16 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 148
+clear
+clc;
+ds=614;
+ns=33725;
+P1=ds/ns;
+n=6000;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.7/Ex7_7.sce b/3513/CH7/EX7.7/Ex7_7.sce new file mode 100644 index 000000000..dca766204 --- /dev/null +++ b/3513/CH7/EX7.7/Ex7_7.sce @@ -0,0 +1,16 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 150
+clear
+clc;
+ds=158;
+ns=2196;
+P1=ds/ns;
+n=54;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.7/ch7_7.sce b/3513/CH7/EX7.7/ch7_7.sce new file mode 100644 index 000000000..dca766204 --- /dev/null +++ b/3513/CH7/EX7.7/ch7_7.sce @@ -0,0 +1,16 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 150
+clear
+clc;
+ds=158;
+ns=2196;
+P1=ds/ns;
+n=54;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.8/Ex7_8.sce b/3513/CH7/EX7.8/Ex7_8.sce new file mode 100644 index 000000000..5a5d104b7 --- /dev/null +++ b/3513/CH7/EX7.8/Ex7_8.sce @@ -0,0 +1,16 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 151
+clear
+clc;
+ds=127;
+ns=1547;
+P1=ds/ns;
+n=188;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.8/ch7_8.sce b/3513/CH7/EX7.8/ch7_8.sce new file mode 100644 index 000000000..5a5d104b7 --- /dev/null +++ b/3513/CH7/EX7.8/ch7_8.sce @@ -0,0 +1,16 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 151
+clear
+clc;
+ds=127;
+ns=1547;
+P1=ds/ns;
+n=188;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.9/Ex7_9.sce b/3513/CH7/EX7.9/Ex7_9.sce new file mode 100644 index 000000000..ea801870c --- /dev/null +++ b/3513/CH7/EX7.9/Ex7_9.sce @@ -0,0 +1,25 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 152
+clear
+clc;
+ds=230;
+ns=4150;
+P1=ds/ns;
+n=200;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+P1=(230-30)/(4150-200);
+//Control limits for p-chart
+mprintf("\Again for \nP1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.4f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH7/EX7.9/ch7_9.sce b/3513/CH7/EX7.9/ch7_9.sce new file mode 100644 index 000000000..ea801870c --- /dev/null +++ b/3513/CH7/EX7.9/ch7_9.sce @@ -0,0 +1,25 @@ +//Calculate the control limits for the p-chart using 3σ limits
+//page no 152
+clear
+clc;
+ds=230;
+ns=4150;
+P1=ds/ns;
+n=200;
+//Control limits for p-chart
+mprintf("\P1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.2f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
+P1=(230-30)/(4150-200);
+//Control limits for p-chart
+mprintf("\Again for \nP1 = %.4f \n",P1);
+UCL=P1+3*sqrt((P1*(1-P1))/n);
+mprintf("\UCL = %.2f \n",UCL);
+LCL=P1-3*sqrt((P1*(1-P1))/n);
+mprintf("\LCL = %.4f \n",LCL);
+CL=P1;
+mprintf("\CL = %.2f \n",CL);
diff --git a/3513/CH9/EX9.1/Ex9_1.sce b/3513/CH9/EX9.1/Ex9_1.sce new file mode 100644 index 000000000..81823dc4e --- /dev/null +++ b/3513/CH9/EX9.1/Ex9_1.sce @@ -0,0 +1,10 @@ +//Calculate the reliability of system RA = 0.8
+//page no 218
+clear
+clc;
+RA = 0.8;
+RB = 0.7;
+RC = 0.9;
+RAB = 1-((1-RA)*(1-RB));
+RS=RAB*RC;
+mprintf("\RS = %.2f \n",RS);
diff --git a/3513/CH9/EX9.1/ch9_1.sce b/3513/CH9/EX9.1/ch9_1.sce new file mode 100644 index 000000000..81823dc4e --- /dev/null +++ b/3513/CH9/EX9.1/ch9_1.sce @@ -0,0 +1,10 @@ +//Calculate the reliability of system RA = 0.8
+//page no 218
+clear
+clc;
+RA = 0.8;
+RB = 0.7;
+RC = 0.9;
+RAB = 1-((1-RA)*(1-RB));
+RS=RAB*RC;
+mprintf("\RS = %.2f \n",RS);
diff --git a/3513/CH9/EX9.2/Ex9_2.sce b/3513/CH9/EX9.2/Ex9_2.sce new file mode 100644 index 000000000..3f9355d79 --- /dev/null +++ b/3513/CH9/EX9.2/Ex9_2.sce @@ -0,0 +1,15 @@ +//Calculate the reliability of the configuration
+//page no 218
+clear
+clc;
+RA = 0.7;
+RB = 0.7;
+RC = 0.9;
+RD = 0.8;
+RE = 0.9;
+
+RAB = 1-((1-RA)*(1-RB));
+RABC=RAB*RC;
+RABCD = 1-((1-RABC)*(1-RD));
+RABCDE = RABCD*RE;
+mprintf("\RABCDE = %.4f \n",RABCDE);
diff --git a/3513/CH9/EX9.2/ch9_2.sce b/3513/CH9/EX9.2/ch9_2.sce new file mode 100644 index 000000000..3f9355d79 --- /dev/null +++ b/3513/CH9/EX9.2/ch9_2.sce @@ -0,0 +1,15 @@ +//Calculate the reliability of the configuration
+//page no 218
+clear
+clc;
+RA = 0.7;
+RB = 0.7;
+RC = 0.9;
+RD = 0.8;
+RE = 0.9;
+
+RAB = 1-((1-RA)*(1-RB));
+RABC=RAB*RC;
+RABCD = 1-((1-RABC)*(1-RD));
+RABCDE = RABCD*RE;
+mprintf("\RABCDE = %.4f \n",RABCDE);
diff --git a/3513/CH9/EX9.3/Ex9_3.sce b/3513/CH9/EX9.3/Ex9_3.sce new file mode 100644 index 000000000..91f8cb150 --- /dev/null +++ b/3513/CH9/EX9.3/Ex9_3.sce @@ -0,0 +1,21 @@ +//Calculate the reliability of the configuration
+//page no 219
+clear
+clc;
+RA = 0.8;
+RB = 0.8;
+RC = 0.8;
+RD = 0.95;
+RE = 0.85;
+
+RABC = 1-((1-RA)*(1-RB)*(1-RC));
+RABCDE1 = RABC*RD*RE;
+mprintf("\RABC = %.4f \n",RABC);
+mprintf("\RABCDE = %.4f \n",RABCDE1);
+RABC=0.992;
+RD=0.95;
+RE = 1-((1-RE)*(1-RE));
+RABCDE2 = RABC*RD*RE;
+mprintf("\(b) WKT RABCDE = %.4f \n",RABCDE2);
+I=RABCDE2-RABCDE1;
+mprintf("\Improvement in R = %.4f \n",I);
diff --git a/3513/CH9/EX9.3/ch9_3.sce b/3513/CH9/EX9.3/ch9_3.sce new file mode 100644 index 000000000..91f8cb150 --- /dev/null +++ b/3513/CH9/EX9.3/ch9_3.sce @@ -0,0 +1,21 @@ +//Calculate the reliability of the configuration
+//page no 219
+clear
+clc;
+RA = 0.8;
+RB = 0.8;
+RC = 0.8;
+RD = 0.95;
+RE = 0.85;
+
+RABC = 1-((1-RA)*(1-RB)*(1-RC));
+RABCDE1 = RABC*RD*RE;
+mprintf("\RABC = %.4f \n",RABC);
+mprintf("\RABCDE = %.4f \n",RABCDE1);
+RABC=0.992;
+RD=0.95;
+RE = 1-((1-RE)*(1-RE));
+RABCDE2 = RABC*RD*RE;
+mprintf("\(b) WKT RABCDE = %.4f \n",RABCDE2);
+I=RABCDE2-RABCDE1;
+mprintf("\Improvement in R = %.4f \n",I);
diff --git a/3513/CH9/EX9.4/Ex9_4.sce b/3513/CH9/EX9.4/Ex9_4.sce new file mode 100644 index 000000000..39a7c683a --- /dev/null +++ b/3513/CH9/EX9.4/Ex9_4.sce @@ -0,0 +1,26 @@ +//Determine the probability of success for the foll.
+//page no 219
+clear
+clc;
+PA=90;
+PB=80;
+PC=75;
+PD=80;
+
+RA = 0.8;
+RB = 0.8;
+RC = 0.8;
+RD = 0.95;
+RE = 0.85;
+
+RABC = 1-((1-RA)*(1-RB)*(1-RC));
+RABCDE1 = RABC*RD*RE;
+mprintf("\RABC = %.4f \n",RABC);
+mprintf("\RABCDE = %.4f \n",RABCDE1);
+RABC=0.992;
+RD=0.95;
+RE = 1-((1-RE)*(1-RE));
+RABCDE2 = RABC*RD*RE;
+mprintf("\(b) WKT RABCDE = %.4f \n",RABCDE2);
+I=RABCDE2-RABCDE1;
+mprintf("\Improvement in R = %.4f \n",I);
diff --git a/3513/CH9/EX9.4/ch9_4.sce b/3513/CH9/EX9.4/ch9_4.sce new file mode 100644 index 000000000..39a7c683a --- /dev/null +++ b/3513/CH9/EX9.4/ch9_4.sce @@ -0,0 +1,26 @@ +//Determine the probability of success for the foll.
+//page no 219
+clear
+clc;
+PA=90;
+PB=80;
+PC=75;
+PD=80;
+
+RA = 0.8;
+RB = 0.8;
+RC = 0.8;
+RD = 0.95;
+RE = 0.85;
+
+RABC = 1-((1-RA)*(1-RB)*(1-RC));
+RABCDE1 = RABC*RD*RE;
+mprintf("\RABC = %.4f \n",RABC);
+mprintf("\RABCDE = %.4f \n",RABCDE1);
+RABC=0.992;
+RD=0.95;
+RE = 1-((1-RE)*(1-RE));
+RABCDE2 = RABC*RD*RE;
+mprintf("\(b) WKT RABCDE = %.4f \n",RABCDE2);
+I=RABCDE2-RABCDE1;
+mprintf("\Improvement in R = %.4f \n",I);
diff --git a/3513/CH9/EX9.5/Ex9_5.sce b/3513/CH9/EX9.5/Ex9_5.sce new file mode 100644 index 000000000..ca13118e5 --- /dev/null +++ b/3513/CH9/EX9.5/Ex9_5.sce @@ -0,0 +1,23 @@ +//Determine the probability of survival for an operating period of 800 hrs
+//page no 219
+clear
+clc;
+SSA = 12500;
+SSC = 11000;
+SSE = 15550;
+SSB = 2830;
+SSD = 9850;
+t=800;
+MTBFA=12500;
+MTBFB=2830;
+MTBFC=11000;
+MTBFD=9850;
+MTBFE=15550;
+dA=1/MTBFA;
+dB=1/MTBFB;
+dC=1/MTBFC;
+dD=1/MTBFD;
+dE=1/MTBFE;
+dS=dA+dB+dC+dD+dE;
+RS =%e^(-dS*t);
+mprintf("\Rs = %.4f \n",RS);
diff --git a/3513/CH9/EX9.5/ch9_5.sce b/3513/CH9/EX9.5/ch9_5.sce new file mode 100644 index 000000000..ca13118e5 --- /dev/null +++ b/3513/CH9/EX9.5/ch9_5.sce @@ -0,0 +1,23 @@ +//Determine the probability of survival for an operating period of 800 hrs
+//page no 219
+clear
+clc;
+SSA = 12500;
+SSC = 11000;
+SSE = 15550;
+SSB = 2830;
+SSD = 9850;
+t=800;
+MTBFA=12500;
+MTBFB=2830;
+MTBFC=11000;
+MTBFD=9850;
+MTBFE=15550;
+dA=1/MTBFA;
+dB=1/MTBFB;
+dC=1/MTBFC;
+dD=1/MTBFD;
+dE=1/MTBFE;
+dS=dA+dB+dC+dD+dE;
+RS =%e^(-dS*t);
+mprintf("\Rs = %.4f \n",RS);
diff --git a/3513/CH9/EX9.6/Ex9_6.sce b/3513/CH9/EX9.6/Ex9_6.sce new file mode 100644 index 000000000..49f6109ba --- /dev/null +++ b/3513/CH9/EX9.6/Ex9_6.sce @@ -0,0 +1,11 @@ +//Determine the probability of survival of the system
+//page no 220
+clear
+clc;
+t=1500;
+dt=1.56*10^-5;
+dr=2*10^-5;
+df=1.7*10^-5;
+ds=dt+dr+df;
+RS =%e^(-ds*t);
+mprintf("\Rs = %.4f \n",RS);
diff --git a/3513/CH9/EX9.6/ch9_6.sce b/3513/CH9/EX9.6/ch9_6.sce new file mode 100644 index 000000000..49f6109ba --- /dev/null +++ b/3513/CH9/EX9.6/ch9_6.sce @@ -0,0 +1,11 @@ +//Determine the probability of survival of the system
+//page no 220
+clear
+clc;
+t=1500;
+dt=1.56*10^-5;
+dr=2*10^-5;
+df=1.7*10^-5;
+ds=dt+dr+df;
+RS =%e^(-ds*t);
+mprintf("\Rs = %.4f \n",RS);
diff --git a/3513/CH9/EX9.7/Ex9_7.sce b/3513/CH9/EX9.7/Ex9_7.sce new file mode 100644 index 000000000..f6f0112e0 --- /dev/null +++ b/3513/CH9/EX9.7/Ex9_7.sce @@ -0,0 +1,21 @@ +//Determine the reliability of the system for 20 hrs. of operating period
+//page no 220
+clear
+clc;
+dA = 0.01;
+dB = 0.015;
+dC = 0.02;
+dD = 0.02;
+dE = 0.025;
+t=20;
+RA =%e^(-dA*t);
+RB =%e^(-dB*t);
+RC =%e^(-dC*t);
+RD =%e^(-dD*t);
+RE =%e^(-dE*t);
+RBC = 1-((1-RB)*(1-RC));
+RABC = RA*RBC;
+RABCD = 1-((1-RABC)*(1-RD));
+RABCDE = RE*RABCD;
+
+mprintf("RABCDE = Rs = %.4f \n",RABCDE);
diff --git a/3513/CH9/EX9.7/ch9_7.sce b/3513/CH9/EX9.7/ch9_7.sce new file mode 100644 index 000000000..f6f0112e0 --- /dev/null +++ b/3513/CH9/EX9.7/ch9_7.sce @@ -0,0 +1,21 @@ +//Determine the reliability of the system for 20 hrs. of operating period
+//page no 220
+clear
+clc;
+dA = 0.01;
+dB = 0.015;
+dC = 0.02;
+dD = 0.02;
+dE = 0.025;
+t=20;
+RA =%e^(-dA*t);
+RB =%e^(-dB*t);
+RC =%e^(-dC*t);
+RD =%e^(-dD*t);
+RE =%e^(-dE*t);
+RBC = 1-((1-RB)*(1-RC));
+RABC = RA*RBC;
+RABCD = 1-((1-RABC)*(1-RD));
+RABCDE = RE*RABCD;
+
+mprintf("RABCDE = Rs = %.4f \n",RABCDE);
diff --git a/3513/CH9/EX9.8/Ex9_8.sce b/3513/CH9/EX9.8/Ex9_8.sce new file mode 100644 index 000000000..0b6efe991 --- /dev/null +++ b/3513/CH9/EX9.8/Ex9_8.sce @@ -0,0 +1,17 @@ +//determine the reliability of the total system for 150 hrs. of operation
+//page no 221
+clear
+clc;
+RS=0.9;
+RM=0.92;
+RC=0.85;
+RB=0.95;
+RP=0.99;
+Rf=0.95;
+Rb1=0.6686;
+Rb2=0.6686;
+Rblower1=RS*RM*RC*RB;
+Rblower2=Rblower1;
+Rblower=1-((1-Rb1)*(1-Rb2));
+Rsystem=Rblower*RP*Rf;
+mprintf("Rsystem = %.4f \n",Rsystem);
diff --git a/3513/CH9/EX9.8/ch9_8.sce b/3513/CH9/EX9.8/ch9_8.sce new file mode 100644 index 000000000..0b6efe991 --- /dev/null +++ b/3513/CH9/EX9.8/ch9_8.sce @@ -0,0 +1,17 @@ +//determine the reliability of the total system for 150 hrs. of operation
+//page no 221
+clear
+clc;
+RS=0.9;
+RM=0.92;
+RC=0.85;
+RB=0.95;
+RP=0.99;
+Rf=0.95;
+Rb1=0.6686;
+Rb2=0.6686;
+Rblower1=RS*RM*RC*RB;
+Rblower2=Rblower1;
+Rblower=1-((1-Rb1)*(1-Rb2));
+Rsystem=Rblower*RP*Rf;
+mprintf("Rsystem = %.4f \n",Rsystem);
|