#if !defined __PFACONST_H
#define __PFACONST_H
#define P120 sin(1.0*M_PI/26.0) // 0.120536680 // P120
#define P142 sin(1.0*M_PI/22.0) // 0.142314838 // P142
#define P173 sin(1.0*M_PI/18.0) // 0.173648178 // P173
#define P222 sin(1.0*M_PI/14.0) // 0.222520934 // P222
#define P239 sin(1.0*M_PI/13.0) // 0.239315664 // P239
#define P281 sin(1.0*M_PI/11.0) // 0.281732557 // P281
#define P342 sin(1.0*M_PI/9.0) // 0.342020143 // P342
#define P354 sin(3.0*M_PI/26.0) // 0.354604887 // P354
#define P382 sin(1.0*M_PI/8.0) // 0.382683432 // P382
#define P415 sin(3.0*M_PI/22.0) // 0.415415013 // P415
#define P433 sin(1.0*M_PI/7.0) // 0.433883739 // P433
#define P464 sin(2.0*M_PI/13.0) // 0.464723172 // P464
#define P540 sin(2.0*M_PI/11.0) // 0.540640817 // P540
#define P559 (1.0/4.0+sin(1.0*M_PI/10.0)) // 0.559016994 // P559
#define P568 sin(5.0*M_PI/26.0) // 0.568064747 // P568
#define P587 sin(1.0*M_PI/5.0) // 0.587785252 // P587
#define P623 sin(3.0*M_PI/14.0) // 0.623489802 // P623
#define P642 sin(2.0*M_PI/9.0) // 0.642787610 // P642
#define P654 sin(5.0*M_PI/22.0) // 0.654860734 // P654
#define P663 sin(3.0*M_PI/13.0) // 0.663122658 // P663
#define P707 sin(1.0*M_PI/4.0) // 0.707106781 // P707
#define P748 sin(7.0*M_PI/26.0) // 0.748510748 // P748
#define P755 sin(3.0*M_PI/11.0) // 0.755749574 // P755
#define P766 sin(5.0*M_PI/18.0) // 0.766044443 // P766
#define P781 sin(2.0*M_PI/7.0) // 0.781831482 // P781
#define P822 sin(4.0*M_PI/13.0) // 0.822983866 // P822
#define P841 sin(7.0*M_PI/22.0) // 0.841253533 // P841
#define P866 sin(1.0*M_PI/3.0) // 0.866025404 // P866
#define P885 sin(9.0*M_PI/26.0) // 0.885456026 // P885
#define P900 sin(5.0*M_PI/14.0) // 0.900968868 // P900
#define P909 sin(4.0*M_PI/11.0) // 0.909631995 // P909
#define P923 sin(3.0*M_PI/8.0) // 0.923879533 // P923
#define P935 sin(5.0*M_PI/13.0) // 0.935016243 // P935
#define P939 sin(7.0*M_PI/18.0) // 0.939692621 // P939
#define P951 sin(2.0*M_PI/5.0) // 0.951056516 // P951
#define P959 sin(9.0*M_PI/22.0) // 0.959492974 // P959
#define P970 sin(11.0*M_PI/26.0) // 0.970941817 // P970
#define P974 sin(3.0*M_PI/7.0) // 0.974927912 // P974
#define P984 sin(4.0*M_PI/9.0) // 0.984807753 // P984
#define P989 sin(5.0*M_PI/11.0) // 0.989821442 // P989
#define P992 sin(6.0*M_PI/13.0) // 0.992708874 // P992
#define NFAX 10 /* # of possible factors */
#define NTAB 240
static struct
{
int n;
double c;
} nctab[NTAB]
=
{
{ 1, 0.000052 },
{ 2, 0.000061 },
{ 3, 0.000030 },
{ 4, 0.000053 },
{ 5, 0.000066 },
{ 6, 0.000067 },
{ 7, 0.000071 },
{ 8, 0.000062 },
{ 9, 0.000079 },
{ 10, 0.000080 },
{ 11, 0.000052 },
{ 12, 0.000069 },
{ 13, 0.000103 },
{ 14, 0.000123 },
{ 15, 0.000050 },
{ 16, 0.000086 },
{ 18, 0.000108 },
{ 20, 0.000101 },
{ 21, 0.000098 },
{ 22, 0.000135 },
{ 24, 0.000090 },
{ 26, 0.000165 },
{ 28, 0.000084 },
{ 30, 0.000132 },
{ 33, 0.000158 },
{ 35, 0.000138 },
{ 36, 0.000147 },
{ 39, 0.000207 },
{ 40, 0.000156 },
{ 42, 0.000158 },
{ 44, 0.000176 },
{ 45, 0.000171 },
{ 48, 0.000185 },
{ 52, 0.000227 },
{ 55, 0.000242 },
{ 56, 0.000194 },
{ 60, 0.000215 },
{ 63, 0.000233 },
{ 65, 0.000288 },
{ 66, 0.000271 },
{ 70, 0.000248 },
{ 72, 0.000247 },
{ 77, 0.000285 },
{ 78, 0.000395 },
{ 80, 0.000285 },
{ 84, 0.000209 },
{ 88, 0.000332 },
{ 90, 0.000321 },
{ 91, 0.000372 },
{ 99, 0.000400 },
{ 104, 0.000391 },
{ 105, 0.000358 },
{ 110, 0.000440 },
{ 112, 0.000367 },
{ 117, 0.000494 },
{ 120, 0.000413 },
{ 126, 0.000424 },
{ 130, 0.000549 },
{ 132, 0.000480 },
{ 140, 0.000450 },
{ 143, 0.000637 },
{ 144, 0.000497 },
{ 154, 0.000590 },
{ 156, 0.000626 },
{ 165, 0.000654 },
{ 168, 0.000536 },
{ 176, 0.000656 },
{ 180, 0.000611 },
{ 182, 0.000730 },
{ 195, 0.000839 },
{ 198, 0.000786 },
{ 208, 0.000835 },
{ 210, 0.000751 },
{ 220, 0.000826 },
{ 231, 0.000926 },
{ 234, 0.000991 },
{ 240, 0.000852 },
{ 252, 0.000820 },
{ 260, 0.001053 },
{ 264, 0.000987 },
{ 273, 0.001152 },
{ 280, 0.000952 },
{ 286, 0.001299 },
{ 308, 0.001155 },
{ 312, 0.001270 },
{ 315, 0.001156 },
{ 330, 0.001397 },
{ 336, 0.001173 },
{ 360, 0.001259 },
{ 364, 0.001471 },
{ 385, 0.001569 },
{ 390, 0.001767 },
{ 396, 0.001552 },
{ 420, 0.001516 },
{ 429, 0.002015 },
{ 440, 0.001748 },
{ 455, 0.001988 },
{ 462, 0.001921 },
{ 468, 0.001956 },
{ 495, 0.002106 },
{ 504, 0.001769 },
{ 520, 0.002196 },
{ 528, 0.002127 },
{ 546, 0.002454 },
{ 560, 0.002099 },
{ 572, 0.002632 },
{ 585, 0.002665 },
{ 616, 0.002397 },
{ 624, 0.002711 },
{ 630, 0.002496 },
{ 660, 0.002812 },
{ 693, 0.002949 },
{ 715, 0.003571 },
{ 720, 0.002783 },
{ 728, 0.003060 },
{ 770, 0.003392 },
{ 780, 0.003553 },
{ 792, 0.003198 },
{ 819, 0.003726 },
{ 840, 0.003234 },
{ 858, 0.004354 },
{ 880, 0.003800 },
{ 910, 0.004304 },
{ 924, 0.003975 },
{ 936, 0.004123 },
{ 990, 0.004517 },
{ 1001, 0.005066 },
{ 1008, 0.003902 },
{ 1040, 0.004785 },
{ 1092, 0.005017 },
{ 1144, 0.005599 },
{ 1155, 0.005380 },
{ 1170, 0.005730 },
{ 1232, 0.005323 },
{ 1260, 0.005112 },
{ 1287, 0.006658 },
{ 1320, 0.005974 },
{ 1365, 0.006781 },
{ 1386, 0.006413 },
{ 1430, 0.007622 },
{ 1456, 0.006679 },
{ 1540, 0.007032 },
{ 1560, 0.007538 },
{ 1584, 0.007126 },
{ 1638, 0.007979 },
{ 1680, 0.007225 },
{ 1716, 0.008961 },
{ 1820, 0.008818 },
{ 1848, 0.008427 },
{ 1872, 0.009004 },
{ 1980, 0.009398 },
{ 2002, 0.010830 },
{ 2145, 0.012010 },
{ 2184, 0.010586 },
{ 2288, 0.012058 },
{ 2310, 0.011673 },
{ 2340, 0.011700 },
{ 2520, 0.011062 },
{ 2574, 0.014313 },
{ 2640, 0.013021 },
{ 2730, 0.014606 },
{ 2772, 0.013216 },
{ 2860, 0.015789 },
{ 3003, 0.016988 },
{ 3080, 0.014911 },
{ 3120, 0.016393 },
{ 3276, 0.016741 },
{ 3432, 0.018821 },
{ 3465, 0.018138 },
{ 3640, 0.018892 },
{ 3696, 0.018634 },
{ 3960, 0.020216 },
{ 4004, 0.022455 },
{ 4095, 0.022523 },
{ 4290, 0.026087 },
{ 4368, 0.023474 },
{ 4620, 0.024590 },
{ 4680, 0.025641 },
{ 5005, 0.030303 },
{ 5040, 0.025253 },
{ 5148, 0.030364 },
{ 5460, 0.031250 },
{ 5544, 0.029412 },
{ 5720, 0.034404 },
{ 6006, 0.037500 },
{ 6160, 0.034091 },
{ 6435, 0.040214 },
{ 6552, 0.037221 },
{ 6864, 0.042735 },
{ 6930, 0.040214 },
{ 7280, 0.042980 },
{ 7920, 0.045872 },
{ 8008, 0.049505 },
{ 8190, 0.049834 },
{ 8580, 0.055762 },
{ 9009, 0.057034 },
{ 9240, 0.054945 },
{ 9360, 0.056818 },
{ 10010, 0.066667 },
{ 10296, 0.065502 },
{ 10920, 0.068182 },
{ 11088, 0.065217 },
{ 11440, 0.075000 },
{ 12012, 0.078534 },
{ 12870, 0.087719 },
{ 13104, 0.081081 },
{ 13860, 0.084270 },
{ 15015, 0.102740 },
{ 16016, 0.106383 },
{ 16380, 0.105634 },
{ 17160, 0.119048 },
{ 18018, 0.123967 },
{ 18480, 0.119048 },
{ 20020, 0.137615 },
{ 20592, 0.140187 },
{ 21840, 0.154639 },
{ 24024, 0.168539 },
{ 25740, 0.180723 },
{ 27720, 0.180723 },
{ 30030, 0.220588 },
{ 32760, 0.241935 },
{ 34320, 0.254237 },
{ 36036, 0.254237 },
{ 40040, 0.288462 },
{ 45045, 0.357143 },
{ 48048, 0.357143 },
{ 51480, 0.384615 },
{ 55440, 0.384615 },
{ 60060, 0.454545 },
{ 65520, 0.517241 },
{ 72072, 0.576923 },
{ 80080, 0.625000 },
{ 90090, 0.833333 },
{ 102960, 0.789474 },
{ 120120, 1.153846 },
{ 144144, 1.153846 },
{ 180180, 1.875000 },
{ 240240, 2.500000 },
{ 360360, 3.750000 },
{ 720720, 7.500000 },
};
#endif // !defined __PFACONST_H