58 static double SN[] = {
59 -8.39167827910303881427E-11,
60 4.62591714427012837309E-8,
61 -9.75759303843632795789E-6,
62 9.76945438170435310816E-4,
63 -4.13470316229406538752E-2,
64 1.00000000000000000302E0,
67 static double SD[] = {
68 2.03269266195951942049E-12,
69 1.27997891179943299903E-9,
70 4.41827842801218905784E-7,
71 9.96412122043875552487E-5,
72 1.42085239326149893930E-2,
73 9.99999999999999996984E-1,
76 static double CN[] = {
77 2.02524002389102268789E-11,
78 -1.35249504915790756375E-8,
79 3.59325051419993077021E-6,
80 -4.74007206873407909465E-4,
81 2.89159652607555242092E-2,
82 -1.00000000000000000080E0,
85 static double CD[] = {
86 4.07746040061880559506E-12,
87 3.06780997581887812692E-9,
88 1.23210355685883423679E-6,
89 3.17442024775032769882E-4,
90 5.10028056236446052392E-2,
91 4.00000000000000000080E0,
94 static double FN4[] = {
95 4.23612862892216586994E0,
96 5.45937717161812843388E0,
97 1.62083287701538329132E0,
98 1.67006611831323023771E-1,
99 6.81020132472518137426E-3,
100 1.08936580650328664411E-4,
101 5.48900223421373614008E-7,
106 8.16496634205391016773E0,
107 7.30828822505564552187E0,
108 1.86792257950184183883E0,
109 1.78792052963149907262E-1,
110 7.01710668322789753610E-3,
111 1.10034357153915731354E-4,
112 5.48900252756255700982E-7,
116 4.55880873470465315206E-1,
117 7.13715274100146711374E-1,
118 1.60300158222319456320E-1,
119 1.16064229408124407915E-2,
120 3.49556442447859055605E-4,
121 4.86215430826454749482E-6,
122 3.20092790091004902806E-8,
123 9.41779576128512936592E-11,
124 9.70507110881952024631E-14,
129 9.17463611873684053703E-1,
130 1.78685545332074536321E-1,
131 1.22253594771971293032E-2,
132 3.58696481881851580297E-4,
133 4.92435064317881464393E-6,
134 3.21956939101046018377E-8,
135 9.43720590350276732376E-11,
136 9.70507110881952025725E-14,
140 8.71001698973114191777E-2,
141 6.11379109952219284151E-1,
142 3.97180296392337498885E-1,
143 7.48527737628469092119E-2,
144 5.38868681462177273157E-3,
145 1.61999794598934024525E-4,
146 1.97963874140963632189E-6,
147 7.82579040744090311069E-9,
152 1.64402202413355338886E0,
153 6.66296701268987968381E-1,
154 9.88771761277688796203E-2,
155 6.22396345441768420760E-3,
156 1.73221081474177119497E-4,
157 2.02659182086343991969E-6,
158 7.82579218933534490868E-9,
162 6.97359953443276214934E-1,
163 3.30410979305632063225E-1,
164 3.84878767649974295920E-2,
165 1.71718239052347903558E-3,
166 3.48941165502279436777E-5,
167 3.47131167084116673800E-7,
168 1.70404452782044526189E-9,
169 3.85945925430276600453E-12,
170 3.14040098946363334640E-15,
175 1.68548898811011640017E0,
176 4.87852258695304967486E-1,
177 4.67913194259625806320E-2,
178 1.90284426674399523638E-3,
179 3.68475504442561108162E-5,
180 3.57043223443740838771E-7,
181 1.72693748966316146736E-9,
182 3.87830166023954706752E-12,
183 3.14040098946363335242E-15,
189 int sici(
double x,
double *si,
double *ci)
191 double z,
c,
s,
f,
g;