summaryrefslogtreecommitdiff
path: root/3513/CH7/EX7.10/Ex7_10.sce
blob: cfd8c54446afa7f583987d26d91ee2da23e5ea13 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
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);