GRASS GIS 7 Programmer's Manual  7.5.svn(2018)-r73122
 All Data Structures Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
shapiro1.c
Go to the documentation of this file.
1 #include <stdio.h>
2 #include <stdlib.h>
3 #include <math.h>
4 #include "local_proto.h"
5 
6 
7 double *Cdhc_shapiro_wilk(double *x, int n)
8 {
9  static double y[2];
10  double a[25], s2, *xcopy;
11  double sumb = 0.0, sumx = 0.0, sumx2 = 0.0;
12  int i, k;
13 
14  if ((xcopy = (double *)malloc(n * sizeof(double))) == NULL) {
15  fprintf(stderr, "Memory error in Cdhc_shapiro_wilk\n");
16  exit(EXIT_FAILURE);
17  }
18 
19  k = n / 2;
20  for (i = 0; i < n; ++i) {
21  xcopy[i] = x[i];
22  sumx += x[i];
23  sumx2 += x[i] * x[i];
24  }
25  s2 = sumx2 - sumx * sumx / n;
26 
27  qsort(xcopy, n, sizeof(double), Cdhc_dcmp);
28 
29  if (n == 3) {
30  a[0] = (double).7071;
31  }
32  else if (n == 4) {
33  a[0] = (double).6872;
34  a[1] = (double).1677;
35  }
36  else if (n == 5) {
37  a[0] = (double).6646;
38  a[1] = (double).2413;
39  }
40  else if (n == 6) {
41  a[0] = (double).6431;
42  a[1] = (double).2806;
43  a[2] = (double).0875;
44  }
45  else if (n == 7) {
46  a[0] = (double).6233;
47  a[1] = (double).3031;
48  a[2] = (double).1401;
49  }
50  else if (n == 8) {
51  a[0] = (double).6052;
52  a[1] = (double).3164;
53  a[2] = (double).1743;
54  a[3] = (double).0561;
55  }
56  else if (n == 9) {
57  a[0] = (double).5888;
58  a[1] = (double).3244;
59  a[2] = (double).1976;
60  a[3] = (double).0947;
61  }
62  else if (n == 10) {
63  a[0] = (double).5739;
64  a[1] = (double).3291;
65  a[2] = (double).2141;
66  a[3] = (double).1224;
67  a[4] = (double).0399;
68  }
69  else if (n == 11) {
70  a[0] = (double).5601;
71  a[1] = (double).3315;
72  a[2] = (double).226;
73  a[3] = (double).1429;
74  a[4] = (double).0695;
75  }
76  else if (n == 12) {
77  a[0] = (double).5475;
78  a[1] = (double).3325;
79  a[2] = (double).2347;
80  a[3] = (double).1586;
81  a[4] = (double).0922;
82  a[5] = (double).0303;
83  }
84  else if (n == 13) {
85  a[0] = (double).5359;
86  a[1] = (double).3325;
87  a[2] = (double).2412;
88  a[3] = (double).1707;
89  a[4] = (double).1099;
90  a[5] = (double).0539;
91  }
92  else if (n == 14) {
93  a[0] = (double).5251;
94  a[1] = (double).3318;
95  a[2] = (double).246;
96  a[3] = (double).1802;
97  a[4] = (double).124;
98  a[5] = (double).0727;
99  a[6] = (double).024;
100  }
101  else if (n == 15) {
102  a[0] = (double).515;
103  a[1] = (double).3306;
104  a[2] = (double).2495;
105  a[3] = (double).1878;
106  a[4] = (double).1353;
107  a[5] = (double).088;
108  a[6] = (double).0433;
109  }
110  else if (n == 16) {
111  a[0] = (double).5056;
112  a[1] = (double).329;
113  a[2] = (double).2521;
114  a[3] = (double).1939;
115  a[4] = (double).1447;
116  a[5] = (double).1005;
117  a[6] = (double).0593;
118  a[7] = (double).0196;
119  }
120  else if (n == 17) {
121  a[0] = (double).4968;
122  a[1] = (double).3273;
123  a[2] = (double).254;
124  a[3] = (double).1988;
125  a[4] = (double).1524;
126  a[5] = (double).1109;
127  a[6] = (double).0725;
128  a[7] = (double).0359;
129  }
130  else if (n == 18) {
131  a[0] = (double).4886;
132  a[1] = (double).3253;
133  a[2] = (double).2553;
134  a[3] = (double).2027;
135  a[4] = (double).1587;
136  a[5] = (double).1197;
137  a[6] = (double).0837;
138  a[7] = (double).0496;
139  a[8] = (double).0163;
140  }
141  else if (n == 19) {
142  a[0] = (double).4808;
143  a[1] = (double).3232;
144  a[2] = (double).2561;
145  a[3] = (double).2059;
146  a[4] = (double).1641;
147  a[5] = (double).1271;
148  a[6] = (double).0932;
149  a[7] = (double).0612;
150  a[8] = (double).0303;
151  }
152  else if (n == 20) {
153  a[0] = (double).4734;
154  a[1] = (double).3211;
155  a[2] = (double).2565;
156  a[3] = (double).2085;
157  a[4] = (double).1686;
158  a[5] = (double).1334;
159  a[6] = (double).1013;
160  a[7] = (double).0711;
161  a[8] = (double).0422;
162  a[9] = (double).014;
163  }
164  else if (n == 21) {
165  a[0] = (double).4643;
166  a[1] = (double).3185;
167  a[2] = (double).2578;
168  a[3] = (double).2119;
169  a[4] = (double).1736;
170  a[5] = (double).1399;
171  a[6] = (double).1092;
172  a[7] = (double).0804;
173  a[8] = (double).053;
174  a[9] = (double).0263;
175  }
176  else if (n == 22) {
177  a[0] = (double).459;
178  a[1] = (double).3156;
179  a[2] = (double).2571;
180  a[3] = (double).2131;
181  a[4] = (double).1764;
182  a[5] = (double).1443;
183  a[6] = (double).115;
184  a[7] = (double).0878;
185  a[8] = (double).0618;
186  a[9] = (double).0368;
187  a[10] = (double).0122;
188  }
189  else if (n == 23) {
190  a[0] = (double).4542;
191  a[1] = (double).3126;
192  a[2] = (double).2563;
193  a[3] = (double).2139;
194  a[4] = (double).1787;
195  a[5] = (double).148;
196  a[6] = (double).1201;
197  a[7] = (double).0941;
198  a[8] = (double).0696;
199  a[9] = (double).0459;
200  a[10] = (double).0228;
201  }
202  else if (n == 24) {
203  a[0] = (double).4493;
204  a[1] = (double).3098;
205  a[2] = (double).2554;
206  a[3] = (double).2145;
207  a[4] = (double).1807;
208  a[5] = (double).1512;
209  a[6] = (double).1245;
210  a[7] = (double).0997;
211  a[8] = (double).0764;
212  a[9] = (double).0539;
213  a[10] = (double).0321;
214  a[11] = (double).0107;
215  }
216  else if (n == 25) {
217  a[0] = (double).445;
218  a[1] = (double).3069;
219  a[2] = (double).2543;
220  a[3] = (double).2148;
221  a[4] = (double).1822;
222  a[5] = (double).1539;
223  a[6] = (double).1283;
224  a[7] = (double).1046;
225  a[8] = (double).0823;
226  a[9] = (double).061;
227  a[10] = (double).0403;
228  a[11] = (double).02;
229  }
230  else if (n == 26) {
231  a[0] = (double).4407;
232  a[1] = (double).3043;
233  a[2] = (double).2533;
234  a[3] = (double).2151;
235  a[4] = (double).1836;
236  a[5] = (double).1563;
237  a[6] = (double).1316;
238  a[7] = (double).1089;
239  a[8] = (double).0876;
240  a[9] = (double).0672;
241  a[10] = (double).0476;
242  a[11] = (double).0284;
243  a[12] = (double).0094;
244  }
245  else if (n == 27) {
246  a[0] = (double).4366;
247  a[1] = (double).3018;
248  a[2] = (double).2522;
249  a[3] = (double).2152;
250  a[4] = (double).1848;
251  a[5] = (double).1584;
252  a[6] = (double).1346;
253  a[7] = (double).1128;
254  a[8] = (double).0923;
255  a[9] = (double).0728;
256  a[10] = (double).054;
257  a[11] = (double).0358;
258  a[12] = (double).0178;
259  }
260  else if (n == 28) {
261  a[0] = (double).4328;
262  a[1] = (double).2992;
263  a[2] = (double).251;
264  a[3] = (double).2151;
265  a[4] = (double).1857;
266  a[5] = (double).1601;
267  a[6] = (double).1372;
268  a[7] = (double).1162;
269  a[8] = (double).0965;
270  a[9] = (double).0778;
271  a[10] = (double).0598;
272  a[11] = (double).0424;
273  a[12] = (double).0253;
274  a[13] = (double).0084;
275  }
276  else if (n == 29) {
277  a[0] = (double).4291;
278  a[1] = (double).2968;
279  a[2] = (double).2499;
280  a[3] = (double).215;
281  a[4] = (double).1864;
282  a[5] = (double).1616;
283  a[6] = (double).1395;
284  a[7] = (double).1192;
285  a[8] = (double).1002;
286  a[9] = (double).0822;
287  a[10] = (double).065;
288  a[11] = (double).0483;
289  a[12] = (double).032;
290  a[13] = (double).0159;
291  }
292  else if (n == 30) {
293  a[0] = (double).4254;
294  a[1] = (double).2944;
295  a[2] = (double).2487;
296  a[3] = (double).2148;
297  a[4] = (double).187;
298  a[5] = (double).163;
299  a[6] = (double).1415;
300  a[7] = (double).1219;
301  a[8] = (double).1036;
302  a[9] = (double).0862;
303  a[10] = (double).0697;
304  a[11] = (double).0537;
305  a[12] = (double).0381;
306  a[13] = (double).0227;
307  a[14] = (double).0076;
308  }
309  else if (n == 31) {
310  a[0] = (double).422;
311  a[1] = (double).2921;
312  a[2] = (double).2475;
313  a[3] = (double).2145;
314  a[4] = (double).1874;
315  a[5] = (double).1641;
316  a[6] = (double).1433;
317  a[7] = (double).1243;
318  a[8] = (double).1066;
319  a[9] = (double).0899;
320  a[10] = (double).0739;
321  a[11] = (double).0585;
322  a[12] = (double).0435;
323  a[13] = (double).0289;
324  a[14] = (double).0144;
325  }
326  else if (n == 32) {
327  a[0] = (double).4188;
328  a[1] = (double).2898;
329  a[2] = (double).2463;
330  a[3] = (double).2141;
331  a[4] = (double).1878;
332  a[5] = (double).1651;
333  a[6] = (double).1449;
334  a[7] = (double).1265;
335  a[8] = (double).1093;
336  a[9] = (double).0931;
337  a[10] = (double).0777;
338  a[11] = (double).0629;
339  a[12] = (double).0485;
340  a[13] = (double).0344;
341  a[14] = (double).0206;
342  a[15] = (double).0068;
343  }
344  else if (n == 33) {
345  a[0] = (double).4156;
346  a[1] = (double).2876;
347  a[2] = (double).2451;
348  a[3] = (double).2137;
349  a[4] = (double).188;
350  a[5] = (double).166;
351  a[6] = (double).1463;
352  a[7] = (double).1284;
353  a[8] = (double).1118;
354  a[9] = (double).0961;
355  a[10] = (double).0812;
356  a[11] = (double).0669;
357  a[12] = (double).053;
358  a[13] = (double).0395;
359  a[14] = (double).0262;
360  a[15] = (double).0131;
361  }
362  else if (n == 34) {
363  a[0] = (double).4127;
364  a[1] = (double).2854;
365  a[2] = (double).2439;
366  a[3] = (double).2132;
367  a[4] = (double).1882;
368  a[5] = (double).1667;
369  a[6] = (double).1475;
370  a[7] = (double).1301;
371  a[8] = (double).114;
372  a[9] = (double).0988;
373  a[10] = (double).0844;
374  a[11] = (double).0706;
375  a[12] = (double).0572;
376  a[13] = (double).0441;
377  a[14] = (double).0314;
378  a[15] = (double).0187;
379  a[16] = (double).0062;
380  }
381  else if (n == 35) {
382  a[0] = (double).4096;
383  a[1] = (double).2834;
384  a[2] = (double).2427;
385  a[3] = (double).2127;
386  a[4] = (double).1883;
387  a[5] = (double).1673;
388  a[6] = (double).1487;
389  a[7] = (double).1317;
390  a[8] = (double).116;
391  a[9] = (double).1013;
392  a[10] = (double).0873;
393  a[11] = (double).0739;
394  a[12] = (double).061;
395  a[13] = (double).0484;
396  a[14] = (double).0361;
397  a[15] = (double).0239;
398  a[16] = (double).0119;
399  }
400  else if (n == 36) {
401  a[0] = (double).4068;
402  a[1] = (double).2813;
403  a[2] = (double).2415;
404  a[3] = (double).2121;
405  a[4] = (double).1883;
406  a[5] = (double).1678;
407  a[6] = (double).1496;
408  a[7] = (double).1331;
409  a[8] = (double).1179;
410  a[9] = (double).1036;
411  a[10] = (double).09;
412  a[11] = (double).077;
413  a[12] = (double).0645;
414  a[13] = (double).0523;
415  a[14] = (double).0404;
416  a[15] = (double).0287;
417  a[16] = (double).0172;
418  a[17] = (double).0057;
419  }
420  else if (n == 37) {
421  a[0] = (double).404;
422  a[1] = (double).2794;
423  a[2] = (double).2403;
424  a[3] = (double).2116;
425  a[4] = (double).1883;
426  a[5] = (double).1683;
427  a[6] = (double).1505;
428  a[7] = (double).1344;
429  a[8] = (double).1196;
430  a[9] = (double).1056;
431  a[10] = (double).0924;
432  a[11] = (double).0798;
433  a[12] = (double).0677;
434  a[13] = (double).0559;
435  a[14] = (double).0444;
436  a[15] = (double).0331;
437  a[16] = (double).022;
438  a[17] = (double).011;
439  }
440  else if (n == 38) {
441  a[0] = (double).4015;
442  a[1] = (double).2774;
443  a[2] = (double).2391;
444  a[3] = (double).211;
445  a[4] = (double).1881;
446  a[5] = (double).1686;
447  a[6] = (double).1513;
448  a[7] = (double).1356;
449  a[8] = (double).1211;
450  a[9] = (double).1075;
451  a[10] = (double).0947;
452  a[11] = (double).0824;
453  a[12] = (double).0706;
454  a[13] = (double).0592;
455  a[14] = (double).0481;
456  a[15] = (double).0372;
457  a[16] = (double).0264;
458  a[17] = (double).0158;
459  a[18] = (double).0053;
460  }
461  else if (n == 39) {
462  a[0] = (double).3989;
463  a[1] = (double).2755;
464  a[2] = (double).238;
465  a[3] = (double).2104;
466  a[4] = (double).188;
467  a[5] = (double).1689;
468  a[6] = (double).152;
469  a[7] = (double).1366;
470  a[8] = (double).1225;
471  a[9] = (double).1092;
472  a[10] = (double).0967;
473  a[11] = (double).0848;
474  a[12] = (double).0733;
475  a[13] = (double).0622;
476  a[14] = (double).0515;
477  a[15] = (double).0409;
478  a[16] = (double).0305;
479  a[17] = (double).0203;
480  a[18] = (double).0101;
481  }
482  else if (n == 40) {
483  a[0] = (double).3964;
484  a[1] = (double).2737;
485  a[2] = (double).2368;
486  a[3] = (double).2098;
487  a[4] = (double).1878;
488  a[5] = (double).1691;
489  a[6] = (double).1526;
490  a[7] = (double).1376;
491  a[8] = (double).1237;
492  a[9] = (double).1108;
493  a[10] = (double).0986;
494  a[11] = (double).087;
495  a[12] = (double).0759;
496  a[13] = (double).0651;
497  a[14] = (double).0546;
498  a[15] = (double).0444;
499  a[16] = (double).0343;
500  a[17] = (double).0244;
501  a[18] = (double).0146;
502  a[19] = (double).0049;
503  }
504  else if (n == 41) {
505  a[0] = (double).394;
506  a[1] = (double).2719;
507  a[2] = (double).2357;
508  a[3] = (double).2091;
509  a[4] = (double).1876;
510  a[5] = (double).1693;
511  a[6] = (double).1531;
512  a[7] = (double).1384;
513  a[8] = (double).1249;
514  a[9] = (double).1123;
515  a[10] = (double).1004;
516  a[11] = (double).0891;
517  a[12] = (double).0782;
518  a[13] = (double).0677;
519  a[14] = (double).0575;
520  a[15] = (double).0476;
521  a[16] = (double).0379;
522  a[17] = (double).0283;
523  a[18] = (double).0188;
524  a[19] = (double).0094;
525  }
526  else if (n == 42) {
527  a[0] = (double).3917;
528  a[1] = (double).2701;
529  a[2] = (double).2345;
530  a[3] = (double).2085;
531  a[4] = (double).1874;
532  a[5] = (double).1694;
533  a[6] = (double).1535;
534  a[7] = (double).1392;
535  a[8] = (double).1259;
536  a[9] = (double).1136;
537  a[10] = (double).102;
538  a[11] = (double).0909;
539  a[12] = (double).0804;
540  a[13] = (double).0701;
541  a[14] = (double).0602;
542  a[15] = (double).0506;
543  a[16] = (double).0411;
544  a[17] = (double).0318;
545  a[18] = (double).0227;
546  a[19] = (double).0136;
547  a[20] = (double).0045;
548  }
549  else if (n == 43) {
550  a[0] = (double).3894;
551  a[1] = (double).2684;
552  a[2] = (double).2334;
553  a[3] = (double).2078;
554  a[4] = (double).1871;
555  a[5] = (double).1695;
556  a[6] = (double).1539;
557  a[7] = (double).1398;
558  a[8] = (double).1269;
559  a[9] = (double).1149;
560  a[10] = (double).1035;
561  a[11] = (double).0927;
562  a[12] = (double).0824;
563  a[13] = (double).0724;
564  a[14] = (double).0628;
565  a[15] = (double).0534;
566  a[16] = (double).0442;
567  a[17] = (double).0352;
568  a[18] = (double).0263;
569  a[19] = (double).0175;
570  a[20] = (double).0087;
571  }
572  else if (n == 44) {
573  a[0] = (double).3872;
574  a[1] = (double).2667;
575  a[2] = (double).2323;
576  a[3] = (double).2072;
577  a[4] = (double).1868;
578  a[5] = (double).1695;
579  a[6] = (double).1542;
580  a[7] = (double).1405;
581  a[8] = (double).1278;
582  a[9] = (double).116;
583  a[10] = (double).1049;
584  a[11] = (double).0943;
585  a[12] = (double).0842;
586  a[13] = (double).0745;
587  a[14] = (double).0651;
588  a[15] = (double).056;
589  a[16] = (double).0471;
590  a[17] = (double).0383;
591  a[18] = (double).0296;
592  a[19] = (double).0211;
593  a[20] = (double).0126;
594  a[21] = (double).0042;
595  }
596  else if (n == 45) {
597  a[0] = (double).385;
598  a[1] = (double).2651;
599  a[2] = (double).2313;
600  a[3] = (double).2065;
601  a[4] = (double).1865;
602  a[5] = (double).1695;
603  a[6] = (double).1545;
604  a[7] = (double).141;
605  a[8] = (double).1286;
606  a[9] = (double).117;
607  a[10] = (double).1062;
608  a[11] = (double).0959;
609  a[12] = (double).086;
610  a[13] = (double).0765;
611  a[14] = (double).0673;
612  a[15] = (double).0584;
613  a[16] = (double).0497;
614  a[17] = (double).0412;
615  a[18] = (double).0328;
616  a[19] = (double).0245;
617  a[20] = (double).0163;
618  a[21] = (double).0081;
619  }
620  else if (n == 46) {
621  a[0] = (double).383;
622  a[1] = (double).2635;
623  a[2] = (double).2302;
624  a[3] = (double).2058;
625  a[4] = (double).1862;
626  a[5] = (double).1695;
627  a[6] = (double).1548;
628  a[7] = (double).1415;
629  a[8] = (double).1293;
630  a[9] = (double).118;
631  a[10] = (double).1073;
632  a[11] = (double).0972;
633  a[12] = (double).0876;
634  a[13] = (double).0783;
635  a[14] = (double).0694;
636  a[15] = (double).0607;
637  a[16] = (double).0522;
638  a[17] = (double).0439;
639  a[18] = (double).0357;
640  a[19] = (double).0277;
641  a[20] = (double).0197;
642  a[21] = (double).0118;
643  a[22] = (double).0039;
644  }
645  else if (n == 47) {
646  a[0] = (double).3808;
647  a[1] = (double).262;
648  a[2] = (double).2291;
649  a[3] = (double).2052;
650  a[4] = (double).1859;
651  a[5] = (double).1695;
652  a[6] = (double).155;
653  a[7] = (double).142;
654  a[8] = (double).13;
655  a[9] = (double).1189;
656  a[10] = (double).1085;
657  a[11] = (double).0986;
658  a[12] = (double).0892;
659  a[13] = (double).0801;
660  a[14] = (double).0713;
661  a[15] = (double).0628;
662  a[16] = (double).0546;
663  a[17] = (double).0465;
664  a[18] = (double).0385;
665  a[19] = (double).0307;
666  a[20] = (double).0229;
667  a[21] = (double).0153;
668  a[22] = (double).0076;
669  }
670  else if (n == 48) {
671  a[0] = (double).3789;
672  a[1] = (double).2604;
673  a[2] = (double).2281;
674  a[3] = (double).2045;
675  a[4] = (double).1855;
676  a[5] = (double).1693;
677  a[6] = (double).1551;
678  a[7] = (double).1423;
679  a[8] = (double).1306;
680  a[9] = (double).1197;
681  a[10] = (double).1095;
682  a[11] = (double).0998;
683  a[12] = (double).0906;
684  a[13] = (double).0817;
685  a[14] = (double).0731;
686  a[15] = (double).0648;
687  a[16] = (double).0568;
688  a[17] = (double).0489;
689  a[18] = (double).0411;
690  a[19] = (double).0335;
691  a[20] = (double).0259;
692  a[21] = (double).0185;
693  a[22] = (double).0111;
694  a[23] = (double).0037;
695  }
696  else if (n == 49) {
697  a[0] = (double).377;
698  a[1] = (double).2589;
699  a[2] = (double).2271;
700  a[3] = (double).2038;
701  a[4] = (double).1851;
702  a[5] = (double).1692;
703  a[6] = (double).1553;
704  a[7] = (double).1427;
705  a[8] = (double).1312;
706  a[9] = (double).1205;
707  a[10] = (double).1105;
708  a[11] = (double).101;
709  a[12] = (double).0919;
710  a[13] = (double).0832;
711  a[14] = (double).0748;
712  a[15] = (double).0667;
713  a[16] = (double).0588;
714  a[17] = (double).0511;
715  a[18] = (double).0436;
716  a[19] = (double).0361;
717  a[20] = (double).0288;
718  a[21] = (double).0215;
719  a[22] = (double).0143;
720  a[23] = (double).0071;
721  }
722  else if (n == 50) {
723  a[0] = (double).3751;
724  a[1] = (double).2574;
725  a[2] = (double).226;
726  a[3] = (double).2032;
727  a[4] = (double).1847;
728  a[5] = (double).1691;
729  a[6] = (double).1554;
730  a[7] = (double).143;
731  a[8] = (double).1317;
732  a[9] = (double).1212;
733  a[10] = (double).1113;
734  a[11] = (double).102;
735  a[12] = (double).0932;
736  a[13] = (double).0846;
737  a[14] = (double).0764;
738  a[15] = (double).0685;
739  a[16] = (double).0608;
740  a[17] = (double).0532;
741  a[18] = (double).0459;
742  a[19] = (double).0386;
743  a[20] = (double).0314;
744  a[21] = (double).0244;
745  a[22] = (double).0174;
746  a[23] = (double).0104;
747  a[24] = (double).0035;
748  }
749 
750  if (n > 50 || n < 3) {
751 #ifdef NOISY
752  fprintf(stdout,
753  " THIS IS THE SHAPIRO-WILK TEST FOR SMALL SAMPLES\n");
754  fprintf(stdout,
755  " THE SAMPLE SIZE MUST BE LESS THAN OR EQUAL TO 50\n");
756 #endif /* NOISY */
757 
758  y[0] = y[1] = 0.0;
759  }
760  else {
761  for (i = 1; i <= k; ++i)
762  sumb += a[i - 1] * (x[n - i + 1] - x[i]);
763 
764  y[0] = sumb * sumb / s2;
765  y[1] = s2;
766 
767 #ifdef NOISY
768  fprintf(stdout, " TEST13 SW(N) =%10.4f\n", y[0]);
769 #endif /* NOISY */
770  }
771 
772  free(xcopy);
773 
774  return y;
775 }
void free(void *)
#define NULL
Definition: ccmath.h:32
#define x
void * malloc(YYSIZE_T)
double * Cdhc_shapiro_wilk(double *x, int n)
Definition: shapiro1.c:7
int Cdhc_dcmp(const void *i, const void *j)
Definition: dcmp.c:1