Context Navigation

source: branches/3.2/sources/HeuristicLab.ExtLibs/HeuristicLab.ALGLIB/2.3.0/ALGLIB-2.3.0/wsr.cs @ 13783

Visit:

Last change on this file since 13783 was 2806, checked in by gkronber, 15 years ago
Added plugin for new version of ALGLIB. #875 (Update ALGLIB sources)
File size: 147.0 KB

Line
1	/*************************************************************************
2	Copyright (c) 2007, Sergey Bochkanov (ALGLIB project).
3
4	>>> SOURCE LICENSE >>>
5	This program is free software; you can redistribute it and/or modify
6	it under the terms of the GNU General Public License as published by
7	the Free Software Foundation (www.fsf.org); either version 2 of the
8	License, or (at your option) any later version.
9
10	This program is distributed in the hope that it will be useful,
11	but WITHOUT ANY WARRANTY; without even the implied warranty of
12	MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13	GNU General Public License for more details.
14
15	A copy of the GNU General Public License is available at
16	http://www.fsf.org/licensing/licenses
17
18	>>> END OF LICENSE >>>
19	*************************************************************************/
20
21	using System;
22
23	namespace alglib
24	{
25	public class wsr
26	{
27	/*************************************************************************
28	Wilcoxon signed-rank test
29
30	This test checks three hypotheses about the median of the given sample.
31	The following tests are performed:
32	* two-tailed test (null hypothesis - the median is equal to the given
33	value)
34	* left-tailed test (null hypothesis - the median is greater than or
35	equal to the given value)
36	* right-tailed test (null hypothesis - the median is less than or
37	equal to the given value)
38
39	Requirements:
40	* the scale of measurement should be ordinal, interval or ratio (i.e.
41	the test could not be applied to nominal variables).
42	* the distribution should be continuous and symmetric relative to its
43	median.
44	* number of distinct values in the X array should be greater than 4
45
46	The test is non-parametric and doesn't require distribution X to be normal
47
48	Input parameters:
49	X - sample. Array whose index goes from 0 to N-1.
50	N - size of the sample.
51	Median - assumed median value.
52
53	Output parameters:
54	BothTails - p-value for two-tailed test.
55	If BothTails is less than the given significance level
56	the null hypothesis is rejected.
57	LeftTail - p-value for left-tailed test.
58	If LeftTail is less than the given significance level,
59	the null hypothesis is rejected.
60	RightTail - p-value for right-tailed test.
61	If RightTail is less than the given significance level
62	the null hypothesis is rejected.
63
64	To calculate p-values, special approximation is used. This method lets us
65	calculate p-values with two decimal places in interval [0.0001, 1].
66
67	"Two decimal places" does not sound very impressive, but in practice the
68	relative error of less than 1% is enough to make a decision.
69
70	There is no approximation outside the [0.0001, 1] interval. Therefore, if
71	the significance level outlies this interval, the test returns 0.0001.
72
73	-- ALGLIB --
74	Copyright 08.09.2006 by Bochkanov Sergey
75	*************************************************************************/
76	public static void wilcoxonsignedranktest(double[] x,
77	int n,
78	double e,
79	ref double bothtails,
80	ref double lefttail,
81	ref double righttail)
82	{
83	int i = 0;
84	int j = 0;
85	int k = 0;
86	int t = 0;
87	double tmp = 0;
88	int tmpi = 0;
89	int ns = 0;
90	double[] r = new double[0];
91	int[] c = new int[0];
92	double w = 0;
93	double p = 0;
94	double mp = 0;
95	double s = 0;
96	double sigma = 0;
97	double mu = 0;
98
99	x = (double[])x.Clone();
100
101
102	//
103	// Prepare
104	//
105	if( n<5 )
106	{
107	bothtails = 1.0;
108	lefttail = 1.0;
109	righttail = 1.0;
110	return;
111	}
112	ns = 0;
113	for(i=0; i<=n-1; i++)
114	{
115	if( (double)(x[i])==(double)(e) )
116	{
117	continue;
118	}
119	x[ns] = x[i];
120	ns = ns+1;
121	}
122	if( ns<5 )
123	{
124	bothtails = 1.0;
125	lefttail = 1.0;
126	righttail = 1.0;
127	return;
128	}
129	r = new double[ns-1+1];
130	c = new int[ns-1+1];
131	for(i=0; i<=ns-1; i++)
132	{
133	r[i] = Math.Abs(x[i]-e);
134	c[i] = i;
135	}
136
137	//
138	// sort {R, C}
139	//
140	if( ns!=1 )
141	{
142	i = 2;
143	do
144	{
145	t = i;
146	while( t!=1 )
147	{
148	k = t/2;
149	if( (double)(r[k-1])>=(double)(r[t-1]) )
150	{
151	t = 1;
152	}
153	else
154	{
155	tmp = r[k-1];
156	r[k-1] = r[t-1];
157	r[t-1] = tmp;
158	tmpi = c[k-1];
159	c[k-1] = c[t-1];
160	c[t-1] = tmpi;
161	t = k;
162	}
163	}
164	i = i+1;
165	}
166	while( i<=ns );
167	i = ns-1;
168	do
169	{
170	tmp = r[i];
171	r[i] = r[0];
172	r[0] = tmp;
173	tmpi = c[i];
174	c[i] = c[0];
175	c[0] = tmpi;
176	t = 1;
177	while( t!=0 )
178	{
179	k = 2*t;
180	if( k>i )
181	{
182	t = 0;
183	}
184	else
185	{
186	if( k<i )
187	{
188	if( (double)(r[k])>(double)(r[k-1]) )
189	{
190	k = k+1;
191	}
192	}
193	if( (double)(r[t-1])>=(double)(r[k-1]) )
194	{
195	t = 0;
196	}
197	else
198	{
199	tmp = r[k-1];
200	r[k-1] = r[t-1];
201	r[t-1] = tmp;
202	tmpi = c[k-1];
203	c[k-1] = c[t-1];
204	c[t-1] = tmpi;
205	t = k;
206	}
207	}
208	}
209	i = i-1;
210	}
211	while( i>=1 );
212	}
213
214	//
215	// compute tied ranks
216	//
217	i = 0;
218	while( i<=ns-1 )
219	{
220	j = i+1;
221	while( j<=ns-1 )
222	{
223	if( (double)(r[j])!=(double)(r[i]) )
224	{
225	break;
226	}
227	j = j+1;
228	}
229	for(k=i; k<=j-1; k++)
230	{
231	r[k] = 1+((double)(i+j-1))/(double)(2);
232	}
233	i = j;
234	}
235
236	//
237	// Compute W+
238	//
239	w = 0;
240	for(i=0; i<=ns-1; i++)
241	{
242	if( (double)(x[c[i]])>(double)(e) )
243	{
244	w = w+r[i];
245	}
246	}
247
248	//
249	// Result
250	//
251	mu = (double)(ns*(ns+1))/(double)(4);
252	sigma = Math.Sqrt((double)(ns(ns+1)(2*ns+1))/(double)(24));
253	s = (w-mu)/sigma;
254	if( (double)(s)<=(double)(0) )
255	{
256	p = Math.Exp(wsigma(-((w-mu)/sigma), ns));
257	mp = 1-Math.Exp(wsigma(-((w-1-mu)/sigma), ns));
258	}
259	else
260	{
261	mp = Math.Exp(wsigma((w-mu)/sigma, ns));
262	p = 1-Math.Exp(wsigma((w+1-mu)/sigma, ns));
263	}
264	bothtails = Math.Max(2*Math.Min(p, mp), 1.0E-4);
265	lefttail = Math.Max(p, 1.0E-4);
266	righttail = Math.Max(mp, 1.0E-4);
267	}
268
269
270	/*************************************************************************
271	Sequential Chebyshev interpolation.
272	*************************************************************************/
273	private static void wcheb(double x,
274	double c,
275	ref double tj,
276	ref double tj1,
277	ref double r)
278	{
279	double t = 0;
280
281	r = r+c*tj;
282	t = 2xtj1-tj;
283	tj = tj1;
284	tj1 = t;
285	}
286
287
288	/*************************************************************************
289	Tail(S, 5)
290	*************************************************************************/
291	private static double w5(double s)
292	{
293	double result = 0;
294	int w = 0;
295	double r = 0;
296
297	w = (int)Math.Round(-(3.708099e+00*s)+7.500000e+00);
298	if( w>=7 )
299	{
300	r = -6.931e-01;
301	}
302	if( w==6 )
303	{
304	r = -9.008e-01;
305	}
306	if( w==5 )
307	{
308	r = -1.163e+00;
309	}
310	if( w==4 )
311	{
312	r = -1.520e+00;
313	}
314	if( w==3 )
315	{
316	r = -1.856e+00;
317	}
318	if( w==2 )
319	{
320	r = -2.367e+00;
321	}
322	if( w==1 )
323	{
324	r = -2.773e+00;
325	}
326	if( w<=0 )
327	{
328	r = -3.466e+00;
329	}
330	result = r;
331	return result;
332	}
333
334
335	/*************************************************************************
336	Tail(S, 6)
337	*************************************************************************/
338	private static double w6(double s)
339	{
340	double result = 0;
341	int w = 0;
342	double r = 0;
343
344	w = (int)Math.Round(-(4.769696e+00*s)+1.050000e+01);
345	if( w>=10 )
346	{
347	r = -6.931e-01;
348	}
349	if( w==9 )
350	{
351	r = -8.630e-01;
352	}
353	if( w==8 )
354	{
355	r = -1.068e+00;
356	}
357	if( w==7 )
358	{
359	r = -1.269e+00;
360	}
361	if( w==6 )
362	{
363	r = -1.520e+00;
364	}
365	if( w==5 )
366	{
367	r = -1.856e+00;
368	}
369	if( w==4 )
370	{
371	r = -2.213e+00;
372	}
373	if( w==3 )
374	{
375	r = -2.549e+00;
376	}
377	if( w==2 )
378	{
379	r = -3.060e+00;
380	}
381	if( w==1 )
382	{
383	r = -3.466e+00;
384	}
385	if( w<=0 )
386	{
387	r = -4.159e+00;
388	}
389	result = r;
390	return result;
391	}
392
393
394	/*************************************************************************
395	Tail(S, 7)
396	*************************************************************************/
397	private static double w7(double s)
398	{
399	double result = 0;
400	int w = 0;
401	double r = 0;
402
403	w = (int)Math.Round(-(5.916080e+00*s)+1.400000e+01);
404	if( w>=14 )
405	{
406	r = -6.325e-01;
407	}
408	if( w==13 )
409	{
410	r = -7.577e-01;
411	}
412	if( w==12 )
413	{
414	r = -9.008e-01;
415	}
416	if( w==11 )
417	{
418	r = -1.068e+00;
419	}
420	if( w==10 )
421	{
422	r = -1.241e+00;
423	}
424	if( w==9 )
425	{
426	r = -1.451e+00;
427	}
428	if( w==8 )
429	{
430	r = -1.674e+00;
431	}
432	if( w==7 )
433	{
434	r = -1.908e+00;
435	}
436	if( w==6 )
437	{
438	r = -2.213e+00;
439	}
440	if( w==5 )
441	{
442	r = -2.549e+00;
443	}
444	if( w==4 )
445	{
446	r = -2.906e+00;
447	}
448	if( w==3 )
449	{
450	r = -3.243e+00;
451	}
452	if( w==2 )
453	{
454	r = -3.753e+00;
455	}
456	if( w==1 )
457	{
458	r = -4.159e+00;
459	}
460	if( w<=0 )
461	{
462	r = -4.852e+00;
463	}
464	result = r;
465	return result;
466	}
467
468
469	/*************************************************************************
470	Tail(S, 8)
471	*************************************************************************/
472	private static double w8(double s)
473	{
474	double result = 0;
475	int w = 0;
476	double r = 0;
477
478	w = (int)Math.Round(-(7.141428e+00*s)+1.800000e+01);
479	if( w>=18 )
480	{
481	r = -6.399e-01;
482	}
483	if( w==17 )
484	{
485	r = -7.494e-01;
486	}
487	if( w==16 )
488	{
489	r = -8.630e-01;
490	}
491	if( w==15 )
492	{
493	r = -9.913e-01;
494	}
495	if( w==14 )
496	{
497	r = -1.138e+00;
498	}
499	if( w==13 )
500	{
501	r = -1.297e+00;
502	}
503	if( w==12 )
504	{
505	r = -1.468e+00;
506	}
507	if( w==11 )
508	{
509	r = -1.653e+00;
510	}
511	if( w==10 )
512	{
513	r = -1.856e+00;
514	}
515	if( w==9 )
516	{
517	r = -2.079e+00;
518	}
519	if( w==8 )
520	{
521	r = -2.326e+00;
522	}
523	if( w==7 )
524	{
525	r = -2.601e+00;
526	}
527	if( w==6 )
528	{
529	r = -2.906e+00;
530	}
531	if( w==5 )
532	{
533	r = -3.243e+00;
534	}
535	if( w==4 )
536	{
537	r = -3.599e+00;
538	}
539	if( w==3 )
540	{
541	r = -3.936e+00;
542	}
543	if( w==2 )
544	{
545	r = -4.447e+00;
546	}
547	if( w==1 )
548	{
549	r = -4.852e+00;
550	}
551	if( w<=0 )
552	{
553	r = -5.545e+00;
554	}
555	result = r;
556	return result;
557	}
558
559
560	/*************************************************************************
561	Tail(S, 9)
562	*************************************************************************/
563	private static double w9(double s)
564	{
565	double result = 0;
566	int w = 0;
567	double r = 0;
568
569	w = (int)Math.Round(-(8.440972e+00*s)+2.250000e+01);
570	if( w>=22 )
571	{
572	r = -6.931e-01;
573	}
574	if( w==21 )
575	{
576	r = -7.873e-01;
577	}
578	if( w==20 )
579	{
580	r = -8.912e-01;
581	}
582	if( w==19 )
583	{
584	r = -1.002e+00;
585	}
586	if( w==18 )
587	{
588	r = -1.120e+00;
589	}
590	if( w==17 )
591	{
592	r = -1.255e+00;
593	}
594	if( w==16 )
595	{
596	r = -1.394e+00;
597	}
598	if( w==15 )
599	{
600	r = -1.547e+00;
601	}
602	if( w==14 )
603	{
604	r = -1.717e+00;
605	}
606	if( w==13 )
607	{
608	r = -1.895e+00;
609	}
610	if( w==12 )
611	{
612	r = -2.079e+00;
613	}
614	if( w==11 )
615	{
616	r = -2.287e+00;
617	}
618	if( w==10 )
619	{
620	r = -2.501e+00;
621	}
622	if( w==9 )
623	{
624	r = -2.742e+00;
625	}
626	if( w==8 )
627	{
628	r = -3.019e+00;
629	}
630	if( w==7 )
631	{
632	r = -3.294e+00;
633	}
634	if( w==6 )
635	{
636	r = -3.599e+00;
637	}
638	if( w==5 )
639	{
640	r = -3.936e+00;
641	}
642	if( w==4 )
643	{
644	r = -4.292e+00;
645	}
646	if( w==3 )
647	{
648	r = -4.629e+00;
649	}
650	if( w==2 )
651	{
652	r = -5.140e+00;
653	}
654	if( w==1 )
655	{
656	r = -5.545e+00;
657	}
658	if( w<=0 )
659	{
660	r = -6.238e+00;
661	}
662	result = r;
663	return result;
664	}
665
666
667	/*************************************************************************
668	Tail(S, 10)
669	*************************************************************************/
670	private static double w10(double s)
671	{
672	double result = 0;
673	int w = 0;
674	double r = 0;
675
676	w = (int)Math.Round(-(9.810708e+00*s)+2.750000e+01);
677	if( w>=27 )
678	{
679	r = -6.931e-01;
680	}
681	if( w==26 )
682	{
683	r = -7.745e-01;
684	}
685	if( w==25 )
686	{
687	r = -8.607e-01;
688	}
689	if( w==24 )
690	{
691	r = -9.551e-01;
692	}
693	if( w==23 )
694	{
695	r = -1.057e+00;
696	}
697	if( w==22 )
698	{
699	r = -1.163e+00;
700	}
701	if( w==21 )
702	{
703	r = -1.279e+00;
704	}
705	if( w==20 )
706	{
707	r = -1.402e+00;
708	}
709	if( w==19 )
710	{
711	r = -1.533e+00;
712	}
713	if( w==18 )
714	{
715	r = -1.674e+00;
716	}
717	if( w==17 )
718	{
719	r = -1.826e+00;
720	}
721	if( w==16 )
722	{
723	r = -1.983e+00;
724	}
725	if( w==15 )
726	{
727	r = -2.152e+00;
728	}
729	if( w==14 )
730	{
731	r = -2.336e+00;
732	}
733	if( w==13 )
734	{
735	r = -2.525e+00;
736	}
737	if( w==12 )
738	{
739	r = -2.727e+00;
740	}
741	if( w==11 )
742	{
743	r = -2.942e+00;
744	}
745	if( w==10 )
746	{
747	r = -3.170e+00;
748	}
749	if( w==9 )
750	{
751	r = -3.435e+00;
752	}
753	if( w==8 )
754	{
755	r = -3.713e+00;
756	}
757	if( w==7 )
758	{
759	r = -3.987e+00;
760	}
761	if( w==6 )
762	{
763	r = -4.292e+00;
764	}
765	if( w==5 )
766	{
767	r = -4.629e+00;
768	}
769	if( w==4 )
770	{
771	r = -4.986e+00;
772	}
773	if( w==3 )
774	{
775	r = -5.322e+00;
776	}
777	if( w==2 )
778	{
779	r = -5.833e+00;
780	}
781	if( w==1 )
782	{
783	r = -6.238e+00;
784	}
785	if( w<=0 )
786	{
787	r = -6.931e+00;
788	}
789	result = r;
790	return result;
791	}
792
793
794	/*************************************************************************
795	Tail(S, 11)
796	*************************************************************************/
797	private static double w11(double s)
798	{
799	double result = 0;
800	int w = 0;
801	double r = 0;
802
803	w = (int)Math.Round(-(1.124722e+01*s)+3.300000e+01);
804	if( w>=33 )
805	{
806	r = -6.595e-01;
807	}
808	if( w==32 )
809	{
810	r = -7.279e-01;
811	}
812	if( w==31 )
813	{
814	r = -8.002e-01;
815	}
816	if( w==30 )
817	{
818	r = -8.782e-01;
819	}
820	if( w==29 )
821	{
822	r = -9.615e-01;
823	}
824	if( w==28 )
825	{
826	r = -1.050e+00;
827	}
828	if( w==27 )
829	{
830	r = -1.143e+00;
831	}
832	if( w==26 )
833	{
834	r = -1.243e+00;
835	}
836	if( w==25 )
837	{
838	r = -1.348e+00;
839	}
840	if( w==24 )
841	{
842	r = -1.459e+00;
843	}
844	if( w==23 )
845	{
846	r = -1.577e+00;
847	}
848	if( w==22 )
849	{
850	r = -1.700e+00;
851	}
852	if( w==21 )
853	{
854	r = -1.832e+00;
855	}
856	if( w==20 )
857	{
858	r = -1.972e+00;
859	}
860	if( w==19 )
861	{
862	r = -2.119e+00;
863	}
864	if( w==18 )
865	{
866	r = -2.273e+00;
867	}
868	if( w==17 )
869	{
870	r = -2.437e+00;
871	}
872	if( w==16 )
873	{
874	r = -2.607e+00;
875	}
876	if( w==15 )
877	{
878	r = -2.788e+00;
879	}
880	if( w==14 )
881	{
882	r = -2.980e+00;
883	}
884	if( w==13 )
885	{
886	r = -3.182e+00;
887	}
888	if( w==12 )
889	{
890	r = -3.391e+00;
891	}
892	if( w==11 )
893	{
894	r = -3.617e+00;
895	}
896	if( w==10 )
897	{
898	r = -3.863e+00;
899	}
900	if( w==9 )
901	{
902	r = -4.128e+00;
903	}
904	if( w==8 )
905	{
906	r = -4.406e+00;
907	}
908	if( w==7 )
909	{
910	r = -4.680e+00;
911	}
912	if( w==6 )
913	{
914	r = -4.986e+00;
915	}
916	if( w==5 )
917	{
918	r = -5.322e+00;
919	}
920	if( w==4 )
921	{
922	r = -5.679e+00;
923	}
924	if( w==3 )
925	{
926	r = -6.015e+00;
927	}
928	if( w==2 )
929	{
930	r = -6.526e+00;
931	}
932	if( w==1 )
933	{
934	r = -6.931e+00;
935	}
936	if( w<=0 )
937	{
938	r = -7.625e+00;
939	}
940	result = r;
941	return result;
942	}
943
944
945	/*************************************************************************
946	Tail(S, 12)
947	*************************************************************************/
948	private static double w12(double s)
949	{
950	double result = 0;
951	int w = 0;
952	double r = 0;
953
954	w = (int)Math.Round(-(1.274755e+01*s)+3.900000e+01);
955	if( w>=39 )
956	{
957	r = -6.633e-01;
958	}
959	if( w==38 )
960	{
961	r = -7.239e-01;
962	}
963	if( w==37 )
964	{
965	r = -7.878e-01;
966	}
967	if( w==36 )
968	{
969	r = -8.556e-01;
970	}
971	if( w==35 )
972	{
973	r = -9.276e-01;
974	}
975	if( w==34 )
976	{
977	r = -1.003e+00;
978	}
979	if( w==33 )
980	{
981	r = -1.083e+00;
982	}
983	if( w==32 )
984	{
985	r = -1.168e+00;
986	}
987	if( w==31 )
988	{
989	r = -1.256e+00;
990	}
991	if( w==30 )
992	{
993	r = -1.350e+00;
994	}
995	if( w==29 )
996	{
997	r = -1.449e+00;
998	}
999	if( w==28 )
1000	{
1001	r = -1.552e+00;
1002	}
1003	if( w==27 )
1004	{
1005	r = -1.660e+00;
1006	}
1007	if( w==26 )
1008	{
1009	r = -1.774e+00;
1010	}
1011	if( w==25 )
1012	{
1013	r = -1.893e+00;
1014	}
1015	if( w==24 )
1016	{
1017	r = -2.017e+00;
1018	}
1019	if( w==23 )
1020	{
1021	r = -2.148e+00;
1022	}
1023	if( w==22 )
1024	{
1025	r = -2.285e+00;
1026	}
1027	if( w==21 )
1028	{
1029	r = -2.429e+00;
1030	}
1031	if( w==20 )
1032	{
1033	r = -2.581e+00;
1034	}
1035	if( w==19 )
1036	{
1037	r = -2.738e+00;
1038	}
1039	if( w==18 )
1040	{
1041	r = -2.902e+00;
1042	}
1043	if( w==17 )
1044	{
1045	r = -3.076e+00;
1046	}
1047	if( w==16 )
1048	{
1049	r = -3.255e+00;
1050	}
1051	if( w==15 )
1052	{
1053	r = -3.443e+00;
1054	}
1055	if( w==14 )
1056	{
1057	r = -3.645e+00;
1058	}
1059	if( w==13 )
1060	{
1061	r = -3.852e+00;
1062	}
1063	if( w==12 )
1064	{
1065	r = -4.069e+00;
1066	}
1067	if( w==11 )
1068	{
1069	r = -4.310e+00;
1070	}
1071	if( w==10 )
1072	{
1073	r = -4.557e+00;
1074	}
1075	if( w==9 )
1076	{
1077	r = -4.821e+00;
1078	}
1079	if( w==8 )
1080	{
1081	r = -5.099e+00;
1082	}
1083	if( w==7 )
1084	{
1085	r = -5.373e+00;
1086	}
1087	if( w==6 )
1088	{
1089	r = -5.679e+00;
1090	}
1091	if( w==5 )
1092	{
1093	r = -6.015e+00;
1094	}
1095	if( w==4 )
1096	{
1097	r = -6.372e+00;
1098	}
1099	if( w==3 )
1100	{
1101	r = -6.708e+00;
1102	}
1103	if( w==2 )
1104	{
1105	r = -7.219e+00;
1106	}
1107	if( w==1 )
1108	{
1109	r = -7.625e+00;
1110	}
1111	if( w<=0 )
1112	{
1113	r = -8.318e+00;
1114	}
1115	result = r;
1116	return result;
1117	}
1118
1119
1120	/*************************************************************************
1121	Tail(S, 13)
1122	*************************************************************************/
1123	private static double w13(double s)
1124	{
1125	double result = 0;
1126	int w = 0;
1127	double r = 0;
1128
1129	w = (int)Math.Round(-(1.430909e+01*s)+4.550000e+01);
1130	if( w>=45 )
1131	{
1132	r = -6.931e-01;
1133	}
1134	if( w==44 )
1135	{
1136	r = -7.486e-01;
1137	}
1138	if( w==43 )
1139	{
1140	r = -8.068e-01;
1141	}
1142	if( w==42 )
1143	{
1144	r = -8.683e-01;
1145	}
1146	if( w==41 )
1147	{
1148	r = -9.328e-01;
1149	}
1150	if( w==40 )
1151	{
1152	r = -1.001e+00;
1153	}
1154	if( w==39 )
1155	{
1156	r = -1.072e+00;
1157	}
1158	if( w==38 )
1159	{
1160	r = -1.146e+00;
1161	}
1162	if( w==37 )
1163	{
1164	r = -1.224e+00;
1165	}
1166	if( w==36 )
1167	{
1168	r = -1.306e+00;
1169	}
1170	if( w==35 )
1171	{
1172	r = -1.392e+00;
1173	}
1174	if( w==34 )
1175	{
1176	r = -1.481e+00;
1177	}
1178	if( w==33 )
1179	{
1180	r = -1.574e+00;
1181	}
1182	if( w==32 )
1183	{
1184	r = -1.672e+00;
1185	}
1186	if( w==31 )
1187	{
1188	r = -1.773e+00;
1189	}
1190	if( w==30 )
1191	{
1192	r = -1.879e+00;
1193	}
1194	if( w==29 )
1195	{
1196	r = -1.990e+00;
1197	}
1198	if( w==28 )
1199	{
1200	r = -2.104e+00;
1201	}
1202	if( w==27 )
1203	{
1204	r = -2.224e+00;
1205	}
1206	if( w==26 )
1207	{
1208	r = -2.349e+00;
1209	}
1210	if( w==25 )
1211	{
1212	r = -2.479e+00;
1213	}
1214	if( w==24 )
1215	{
1216	r = -2.614e+00;
1217	}
1218	if( w==23 )
1219	{
1220	r = -2.755e+00;
1221	}
1222	if( w==22 )
1223	{
1224	r = -2.902e+00;
1225	}
1226	if( w==21 )
1227	{
1228	r = -3.055e+00;
1229	}
1230	if( w==20 )
1231	{
1232	r = -3.215e+00;
1233	}
1234	if( w==19 )
1235	{
1236	r = -3.380e+00;
1237	}
1238	if( w==18 )
1239	{
1240	r = -3.551e+00;
1241	}
1242	if( w==17 )
1243	{
1244	r = -3.733e+00;
1245	}
1246	if( w==16 )
1247	{
1248	r = -3.917e+00;
1249	}
1250	if( w==15 )
1251	{
1252	r = -4.113e+00;
1253	}
1254	if( w==14 )
1255	{
1256	r = -4.320e+00;
1257	}
1258	if( w==13 )
1259	{
1260	r = -4.534e+00;
1261	}
1262	if( w==12 )
1263	{
1264	r = -4.762e+00;
1265	}
1266	if( w==11 )
1267	{
1268	r = -5.004e+00;
1269	}
1270	if( w==10 )
1271	{
1272	r = -5.250e+00;
1273	}
1274	if( w==9 )
1275	{
1276	r = -5.514e+00;
1277	}
1278	if( w==8 )
1279	{
1280	r = -5.792e+00;
1281	}
1282	if( w==7 )
1283	{
1284	r = -6.066e+00;
1285	}
1286	if( w==6 )
1287	{
1288	r = -6.372e+00;
1289	}
1290	if( w==5 )
1291	{
1292	r = -6.708e+00;
1293	}
1294	if( w==4 )
1295	{
1296	r = -7.065e+00;
1297	}
1298	if( w==3 )
1299	{
1300	r = -7.401e+00;
1301	}
1302	if( w==2 )
1303	{
1304	r = -7.912e+00;
1305	}
1306	if( w==1 )
1307	{
1308	r = -8.318e+00;
1309	}
1310	if( w<=0 )
1311	{
1312	r = -9.011e+00;
1313	}
1314	result = r;
1315	return result;
1316	}
1317
1318
1319	/*************************************************************************
1320	Tail(S, 14)
1321	*************************************************************************/
1322	private static double w14(double s)
1323	{
1324	double result = 0;
1325	int w = 0;
1326	double r = 0;
1327
1328	w = (int)Math.Round(-(1.592953e+01*s)+5.250000e+01);
1329	if( w>=52 )
1330	{
1331	r = -6.931e-01;
1332	}
1333	if( w==51 )
1334	{
1335	r = -7.428e-01;
1336	}
1337	if( w==50 )
1338	{
1339	r = -7.950e-01;
1340	}
1341	if( w==49 )
1342	{
1343	r = -8.495e-01;
1344	}
1345	if( w==48 )
1346	{
1347	r = -9.067e-01;
1348	}
1349	if( w==47 )
1350	{
1351	r = -9.664e-01;
1352	}
1353	if( w==46 )
1354	{
1355	r = -1.029e+00;
1356	}
1357	if( w==45 )
1358	{
1359	r = -1.094e+00;
1360	}
1361	if( w==44 )
1362	{
1363	r = -1.162e+00;
1364	}
1365	if( w==43 )
1366	{
1367	r = -1.233e+00;
1368	}
1369	if( w==42 )
1370	{
1371	r = -1.306e+00;
1372	}
1373	if( w==41 )
1374	{
1375	r = -1.383e+00;
1376	}
1377	if( w==40 )
1378	{
1379	r = -1.463e+00;
1380	}
1381	if( w==39 )
1382	{
1383	r = -1.546e+00;
1384	}
1385	if( w==38 )
1386	{
1387	r = -1.632e+00;
1388	}
1389	if( w==37 )
1390	{
1391	r = -1.722e+00;
1392	}
1393	if( w==36 )
1394	{
1395	r = -1.815e+00;
1396	}
1397	if( w==35 )
1398	{
1399	r = -1.911e+00;
1400	}
1401	if( w==34 )
1402	{
1403	r = -2.011e+00;
1404	}
1405	if( w==33 )
1406	{
1407	r = -2.115e+00;
1408	}
1409	if( w==32 )
1410	{
1411	r = -2.223e+00;
1412	}
1413	if( w==31 )
1414	{
1415	r = -2.334e+00;
1416	}
1417	if( w==30 )
1418	{
1419	r = -2.450e+00;
1420	}
1421	if( w==29 )
1422	{
1423	r = -2.570e+00;
1424	}
1425	if( w==28 )
1426	{
1427	r = -2.694e+00;
1428	}
1429	if( w==27 )
1430	{
1431	r = -2.823e+00;
1432	}
1433	if( w==26 )
1434	{
1435	r = -2.956e+00;
1436	}
1437	if( w==25 )
1438	{
1439	r = -3.095e+00;
1440	}
1441	if( w==24 )
1442	{
1443	r = -3.238e+00;
1444	}
1445	if( w==23 )
1446	{
1447	r = -3.387e+00;
1448	}
1449	if( w==22 )
1450	{
1451	r = -3.541e+00;
1452	}
1453	if( w==21 )
1454	{
1455	r = -3.700e+00;
1456	}
1457	if( w==20 )
1458	{
1459	r = -3.866e+00;
1460	}
1461	if( w==19 )
1462	{
1463	r = -4.038e+00;
1464	}
1465	if( w==18 )
1466	{
1467	r = -4.215e+00;
1468	}
1469	if( w==17 )
1470	{
1471	r = -4.401e+00;
1472	}
1473	if( w==16 )
1474	{
1475	r = -4.592e+00;
1476	}
1477	if( w==15 )
1478	{
1479	r = -4.791e+00;
1480	}
1481	if( w==14 )
1482	{
1483	r = -5.004e+00;
1484	}
1485	if( w==13 )
1486	{
1487	r = -5.227e+00;
1488	}
1489	if( w==12 )
1490	{
1491	r = -5.456e+00;
1492	}
1493	if( w==11 )
1494	{
1495	r = -5.697e+00;
1496	}
1497	if( w==10 )
1498	{
1499	r = -5.943e+00;
1500	}
1501	if( w==9 )
1502	{
1503	r = -6.208e+00;
1504	}
1505	if( w==8 )
1506	{
1507	r = -6.485e+00;
1508	}
1509	if( w==7 )
1510	{
1511	r = -6.760e+00;
1512	}
1513	if( w==6 )
1514	{
1515	r = -7.065e+00;
1516	}
1517	if( w==5 )
1518	{
1519	r = -7.401e+00;
1520	}
1521	if( w==4 )
1522	{
1523	r = -7.758e+00;
1524	}
1525	if( w==3 )
1526	{
1527	r = -8.095e+00;
1528	}
1529	if( w==2 )
1530	{
1531	r = -8.605e+00;
1532	}
1533	if( w==1 )
1534	{
1535	r = -9.011e+00;
1536	}
1537	if( w<=0 )
1538	{
1539	r = -9.704e+00;
1540	}
1541	result = r;
1542	return result;
1543	}
1544
1545
1546	/*************************************************************************
1547	Tail(S, 15)
1548	*************************************************************************/
1549	private static double w15(double s)
1550	{
1551	double result = 0;
1552	int w = 0;
1553	double r = 0;
1554
1555	w = (int)Math.Round(-(1.760682e+01*s)+6.000000e+01);
1556	if( w>=60 )
1557	{
1558	r = -6.714e-01;
1559	}
1560	if( w==59 )
1561	{
1562	r = -7.154e-01;
1563	}
1564	if( w==58 )
1565	{
1566	r = -7.613e-01;
1567	}
1568	if( w==57 )
1569	{
1570	r = -8.093e-01;
1571	}
1572	if( w==56 )
1573	{
1574	r = -8.593e-01;
1575	}
1576	if( w==55 )
1577	{
1578	r = -9.114e-01;
1579	}
1580	if( w==54 )
1581	{
1582	r = -9.656e-01;
1583	}
1584	if( w==53 )
1585	{
1586	r = -1.022e+00;
1587	}
1588	if( w==52 )
1589	{
1590	r = -1.081e+00;
1591	}
1592	if( w==51 )
1593	{
1594	r = -1.142e+00;
1595	}
1596	if( w==50 )
1597	{
1598	r = -1.205e+00;
1599	}
1600	if( w==49 )
1601	{
1602	r = -1.270e+00;
1603	}
1604	if( w==48 )
1605	{
1606	r = -1.339e+00;
1607	}
1608	if( w==47 )
1609	{
1610	r = -1.409e+00;
1611	}
1612	if( w==46 )
1613	{
1614	r = -1.482e+00;
1615	}
1616	if( w==45 )
1617	{
1618	r = -1.558e+00;
1619	}
1620	if( w==44 )
1621	{
1622	r = -1.636e+00;
1623	}
1624	if( w==43 )
1625	{
1626	r = -1.717e+00;
1627	}
1628	if( w==42 )
1629	{
1630	r = -1.801e+00;
1631	}
1632	if( w==41 )
1633	{
1634	r = -1.888e+00;
1635	}
1636	if( w==40 )
1637	{
1638	r = -1.977e+00;
1639	}
1640	if( w==39 )
1641	{
1642	r = -2.070e+00;
1643	}
1644	if( w==38 )
1645	{
1646	r = -2.166e+00;
1647	}
1648	if( w==37 )
1649	{
1650	r = -2.265e+00;
1651	}
1652	if( w==36 )
1653	{
1654	r = -2.366e+00;
1655	}
1656	if( w==35 )
1657	{
1658	r = -2.472e+00;
1659	}
1660	if( w==34 )
1661	{
1662	r = -2.581e+00;
1663	}
1664	if( w==33 )
1665	{
1666	r = -2.693e+00;
1667	}
1668	if( w==32 )
1669	{
1670	r = -2.809e+00;
1671	}
1672	if( w==31 )
1673	{
1674	r = -2.928e+00;
1675	}
1676	if( w==30 )
1677	{
1678	r = -3.051e+00;
1679	}
1680	if( w==29 )
1681	{
1682	r = -3.179e+00;
1683	}
1684	if( w==28 )
1685	{
1686	r = -3.310e+00;
1687	}
1688	if( w==27 )
1689	{
1690	r = -3.446e+00;
1691	}
1692	if( w==26 )
1693	{
1694	r = -3.587e+00;
1695	}
1696	if( w==25 )
1697	{
1698	r = -3.732e+00;
1699	}
1700	if( w==24 )
1701	{
1702	r = -3.881e+00;
1703	}
1704	if( w==23 )
1705	{
1706	r = -4.036e+00;
1707	}
1708	if( w==22 )
1709	{
1710	r = -4.195e+00;
1711	}
1712	if( w==21 )
1713	{
1714	r = -4.359e+00;
1715	}
1716	if( w==20 )
1717	{
1718	r = -4.531e+00;
1719	}
1720	if( w==19 )
1721	{
1722	r = -4.707e+00;
1723	}
1724	if( w==18 )
1725	{
1726	r = -4.888e+00;
1727	}
1728	if( w==17 )
1729	{
1730	r = -5.079e+00;
1731	}
1732	if( w==16 )
1733	{
1734	r = -5.273e+00;
1735	}
1736	if( w==15 )
1737	{
1738	r = -5.477e+00;
1739	}
1740	if( w==14 )
1741	{
1742	r = -5.697e+00;
1743	}
1744	if( w==13 )
1745	{
1746	r = -5.920e+00;
1747	}
1748	if( w==12 )
1749	{
1750	r = -6.149e+00;
1751	}
1752	if( w==11 )
1753	{
1754	r = -6.390e+00;
1755	}
1756	if( w==10 )
1757	{
1758	r = -6.636e+00;
1759	}
1760	if( w==9 )
1761	{
1762	r = -6.901e+00;
1763	}
1764	if( w==8 )
1765	{
1766	r = -7.178e+00;
1767	}
1768	if( w==7 )
1769	{
1770	r = -7.453e+00;
1771	}
1772	if( w==6 )
1773	{
1774	r = -7.758e+00;
1775	}
1776	if( w==5 )
1777	{
1778	r = -8.095e+00;
1779	}
1780	if( w==4 )
1781	{
1782	r = -8.451e+00;
1783	}
1784	if( w==3 )
1785	{
1786	r = -8.788e+00;
1787	}
1788	if( w==2 )
1789	{
1790	r = -9.299e+00;
1791	}
1792	if( w==1 )
1793	{
1794	r = -9.704e+00;
1795	}
1796	if( w<=0 )
1797	{
1798	r = -1.040e+01;
1799	}
1800	result = r;
1801	return result;
1802	}
1803
1804
1805	/*************************************************************************
1806	Tail(S, 16)
1807	*************************************************************************/
1808	private static double w16(double s)
1809	{
1810	double result = 0;
1811	int w = 0;
1812	double r = 0;
1813
1814	w = (int)Math.Round(-(1.933908e+01*s)+6.800000e+01);
1815	if( w>=68 )
1816	{
1817	r = -6.733e-01;
1818	}
1819	if( w==67 )
1820	{
1821	r = -7.134e-01;
1822	}
1823	if( w==66 )
1824	{
1825	r = -7.551e-01;
1826	}
1827	if( w==65 )
1828	{
1829	r = -7.986e-01;
1830	}
1831	if( w==64 )
1832	{
1833	r = -8.437e-01;
1834	}
1835	if( w==63 )
1836	{
1837	r = -8.905e-01;
1838	}
1839	if( w==62 )
1840	{
1841	r = -9.391e-01;
1842	}
1843	if( w==61 )
1844	{
1845	r = -9.895e-01;
1846	}
1847	if( w==60 )
1848	{
1849	r = -1.042e+00;
1850	}
1851	if( w==59 )
1852	{
1853	r = -1.096e+00;
1854	}
1855	if( w==58 )
1856	{
1857	r = -1.152e+00;
1858	}
1859	if( w==57 )
1860	{
1861	r = -1.210e+00;
1862	}
1863	if( w==56 )
1864	{
1865	r = -1.270e+00;
1866	}
1867	if( w==55 )
1868	{
1869	r = -1.331e+00;
1870	}
1871	if( w==54 )
1872	{
1873	r = -1.395e+00;
1874	}
1875	if( w==53 )
1876	{
1877	r = -1.462e+00;
1878	}
1879	if( w==52 )
1880	{
1881	r = -1.530e+00;
1882	}
1883	if( w==51 )
1884	{
1885	r = -1.600e+00;
1886	}
1887	if( w==50 )
1888	{
1889	r = -1.673e+00;
1890	}
1891	if( w==49 )
1892	{
1893	r = -1.748e+00;
1894	}
1895	if( w==48 )
1896	{
1897	r = -1.825e+00;
1898	}
1899	if( w==47 )
1900	{
1901	r = -1.904e+00;
1902	}
1903	if( w==46 )
1904	{
1905	r = -1.986e+00;
1906	}
1907	if( w==45 )
1908	{
1909	r = -2.071e+00;
1910	}
1911	if( w==44 )
1912	{
1913	r = -2.158e+00;
1914	}
1915	if( w==43 )
1916	{
1917	r = -2.247e+00;
1918	}
1919	if( w==42 )
1920	{
1921	r = -2.339e+00;
1922	}
1923	if( w==41 )
1924	{
1925	r = -2.434e+00;
1926	}
1927	if( w==40 )
1928	{
1929	r = -2.532e+00;
1930	}
1931	if( w==39 )
1932	{
1933	r = -2.632e+00;
1934	}
1935	if( w==38 )
1936	{
1937	r = -2.735e+00;
1938	}
1939	if( w==37 )
1940	{
1941	r = -2.842e+00;
1942	}
1943	if( w==36 )
1944	{
1945	r = -2.951e+00;
1946	}
1947	if( w==35 )
1948	{
1949	r = -3.064e+00;
1950	}
1951	if( w==34 )
1952	{
1953	r = -3.179e+00;
1954	}
1955	if( w==33 )
1956	{
1957	r = -3.298e+00;
1958	}
1959	if( w==32 )
1960	{
1961	r = -3.420e+00;
1962	}
1963	if( w==31 )
1964	{
1965	r = -3.546e+00;
1966	}
1967	if( w==30 )
1968	{
1969	r = -3.676e+00;
1970	}
1971	if( w==29 )
1972	{
1973	r = -3.810e+00;
1974	}
1975	if( w==28 )
1976	{
1977	r = -3.947e+00;
1978	}
1979	if( w==27 )
1980	{
1981	r = -4.088e+00;
1982	}
1983	if( w==26 )
1984	{
1985	r = -4.234e+00;
1986	}
1987	if( w==25 )
1988	{
1989	r = -4.383e+00;
1990	}
1991	if( w==24 )
1992	{
1993	r = -4.538e+00;
1994	}
1995	if( w==23 )
1996	{
1997	r = -4.697e+00;
1998	}
1999	if( w==22 )
2000	{
2001	r = -4.860e+00;
2002	}
2003	if( w==21 )
2004	{
2005	r = -5.029e+00;
2006	}
2007	if( w==20 )
2008	{
2009	r = -5.204e+00;
2010	}
2011	if( w==19 )
2012	{
2013	r = -5.383e+00;
2014	}
2015	if( w==18 )
2016	{
2017	r = -5.569e+00;
2018	}
2019	if( w==17 )
2020	{
2021	r = -5.762e+00;
2022	}
2023	if( w==16 )
2024	{
2025	r = -5.960e+00;
2026	}
2027	if( w==15 )
2028	{
2029	r = -6.170e+00;
2030	}
2031	if( w==14 )
2032	{
2033	r = -6.390e+00;
2034	}
2035	if( w==13 )
2036	{
2037	r = -6.613e+00;
2038	}
2039	if( w==12 )
2040	{
2041	r = -6.842e+00;
2042	}
2043	if( w==11 )
2044	{
2045	r = -7.083e+00;
2046	}
2047	if( w==10 )
2048	{
2049	r = -7.329e+00;
2050	}
2051	if( w==9 )
2052	{
2053	r = -7.594e+00;
2054	}
2055	if( w==8 )
2056	{
2057	r = -7.871e+00;
2058	}
2059	if( w==7 )
2060	{
2061	r = -8.146e+00;
2062	}
2063	if( w==6 )
2064	{
2065	r = -8.451e+00;
2066	}
2067	if( w==5 )
2068	{
2069	r = -8.788e+00;
2070	}
2071	if( w==4 )
2072	{
2073	r = -9.144e+00;
2074	}
2075	if( w==3 )
2076	{
2077	r = -9.481e+00;
2078	}
2079	if( w==2 )
2080	{
2081	r = -9.992e+00;
2082	}
2083	if( w==1 )
2084	{
2085	r = -1.040e+01;
2086	}
2087	if( w<=0 )
2088	{
2089	r = -1.109e+01;
2090	}
2091	result = r;
2092	return result;
2093	}
2094
2095
2096	/*************************************************************************
2097	Tail(S, 17)
2098	*************************************************************************/
2099	private static double w17(double s)
2100	{
2101	double result = 0;
2102	int w = 0;
2103	double r = 0;
2104
2105	w = (int)Math.Round(-(2.112463e+01*s)+7.650000e+01);
2106	if( w>=76 )
2107	{
2108	r = -6.931e-01;
2109	}
2110	if( w==75 )
2111	{
2112	r = -7.306e-01;
2113	}
2114	if( w==74 )
2115	{
2116	r = -7.695e-01;
2117	}
2118	if( w==73 )
2119	{
2120	r = -8.097e-01;
2121	}
2122	if( w==72 )
2123	{
2124	r = -8.514e-01;
2125	}
2126	if( w==71 )
2127	{
2128	r = -8.946e-01;
2129	}
2130	if( w==70 )
2131	{
2132	r = -9.392e-01;
2133	}
2134	if( w==69 )
2135	{
2136	r = -9.853e-01;
2137	}
2138	if( w==68 )
2139	{
2140	r = -1.033e+00;
2141	}
2142	if( w==67 )
2143	{
2144	r = -1.082e+00;
2145	}
2146	if( w==66 )
2147	{
2148	r = -1.133e+00;
2149	}
2150	if( w==65 )
2151	{
2152	r = -1.185e+00;
2153	}
2154	if( w==64 )
2155	{
2156	r = -1.240e+00;
2157	}
2158	if( w==63 )
2159	{
2160	r = -1.295e+00;
2161	}
2162	if( w==62 )
2163	{
2164	r = -1.353e+00;
2165	}
2166	if( w==61 )
2167	{
2168	r = -1.412e+00;
2169	}
2170	if( w==60 )
2171	{
2172	r = -1.473e+00;
2173	}
2174	if( w==59 )
2175	{
2176	r = -1.536e+00;
2177	}
2178	if( w==58 )
2179	{
2180	r = -1.600e+00;
2181	}
2182	if( w==57 )
2183	{
2184	r = -1.666e+00;
2185	}
2186	if( w==56 )
2187	{
2188	r = -1.735e+00;
2189	}
2190	if( w==55 )
2191	{
2192	r = -1.805e+00;
2193	}
2194	if( w==54 )
2195	{
2196	r = -1.877e+00;
2197	}
2198	if( w==53 )
2199	{
2200	r = -1.951e+00;
2201	}
2202	if( w==52 )
2203	{
2204	r = -2.028e+00;
2205	}
2206	if( w==51 )
2207	{
2208	r = -2.106e+00;
2209	}
2210	if( w==50 )
2211	{
2212	r = -2.186e+00;
2213	}
2214	if( w==49 )
2215	{
2216	r = -2.269e+00;
2217	}
2218	if( w==48 )
2219	{
2220	r = -2.353e+00;
2221	}
2222	if( w==47 )
2223	{
2224	r = -2.440e+00;
2225	}
2226	if( w==46 )
2227	{
2228	r = -2.530e+00;
2229	}
2230	if( w==45 )
2231	{
2232	r = -2.621e+00;
2233	}
2234	if( w==44 )
2235	{
2236	r = -2.715e+00;
2237	}
2238	if( w==43 )
2239	{
2240	r = -2.812e+00;
2241	}
2242	if( w==42 )
2243	{
2244	r = -2.911e+00;
2245	}
2246	if( w==41 )
2247	{
2248	r = -3.012e+00;
2249	}
2250	if( w==40 )
2251	{
2252	r = -3.116e+00;
2253	}
2254	if( w==39 )
2255	{
2256	r = -3.223e+00;
2257	}
2258	if( w==38 )
2259	{
2260	r = -3.332e+00;
2261	}
2262	if( w==37 )
2263	{
2264	r = -3.445e+00;
2265	}
2266	if( w==36 )
2267	{
2268	r = -3.560e+00;
2269	}
2270	if( w==35 )
2271	{
2272	r = -3.678e+00;
2273	}
2274	if( w==34 )
2275	{
2276	r = -3.799e+00;
2277	}
2278	if( w==33 )
2279	{
2280	r = -3.924e+00;
2281	}
2282	if( w==32 )
2283	{
2284	r = -4.052e+00;
2285	}
2286	if( w==31 )
2287	{
2288	r = -4.183e+00;
2289	}
2290	if( w==30 )
2291	{
2292	r = -4.317e+00;
2293	}
2294	if( w==29 )
2295	{
2296	r = -4.456e+00;
2297	}
2298	if( w==28 )
2299	{
2300	r = -4.597e+00;
2301	}
2302	if( w==27 )
2303	{
2304	r = -4.743e+00;
2305	}
2306	if( w==26 )
2307	{
2308	r = -4.893e+00;
2309	}
2310	if( w==25 )
2311	{
2312	r = -5.047e+00;
2313	}
2314	if( w==24 )
2315	{
2316	r = -5.204e+00;
2317	}
2318	if( w==23 )
2319	{
2320	r = -5.367e+00;
2321	}
2322	if( w==22 )
2323	{
2324	r = -5.534e+00;
2325	}
2326	if( w==21 )
2327	{
2328	r = -5.706e+00;
2329	}
2330	if( w==20 )
2331	{
2332	r = -5.884e+00;
2333	}
2334	if( w==19 )
2335	{
2336	r = -6.066e+00;
2337	}
2338	if( w==18 )
2339	{
2340	r = -6.254e+00;
2341	}
2342	if( w==17 )
2343	{
2344	r = -6.451e+00;
2345	}
2346	if( w==16 )
2347	{
2348	r = -6.654e+00;
2349	}
2350	if( w==15 )
2351	{
2352	r = -6.864e+00;
2353	}
2354	if( w==14 )
2355	{
2356	r = -7.083e+00;
2357	}
2358	if( w==13 )
2359	{
2360	r = -7.306e+00;
2361	}
2362	if( w==12 )
2363	{
2364	r = -7.535e+00;
2365	}
2366	if( w==11 )
2367	{
2368	r = -7.776e+00;
2369	}
2370	if( w==10 )
2371	{
2372	r = -8.022e+00;
2373	}
2374	if( w==9 )
2375	{
2376	r = -8.287e+00;
2377	}
2378	if( w==8 )
2379	{
2380	r = -8.565e+00;
2381	}
2382	if( w==7 )
2383	{
2384	r = -8.839e+00;
2385	}
2386	if( w==6 )
2387	{
2388	r = -9.144e+00;
2389	}
2390	if( w==5 )
2391	{
2392	r = -9.481e+00;
2393	}
2394	if( w==4 )
2395	{
2396	r = -9.838e+00;
2397	}
2398	if( w==3 )
2399	{
2400	r = -1.017e+01;
2401	}
2402	if( w==2 )
2403	{
2404	r = -1.068e+01;
2405	}
2406	if( w==1 )
2407	{
2408	r = -1.109e+01;
2409	}
2410	if( w<=0 )
2411	{
2412	r = -1.178e+01;
2413	}
2414	result = r;
2415	return result;
2416	}
2417
2418
2419	/*************************************************************************
2420	Tail(S, 18)
2421	*************************************************************************/
2422	private static double w18(double s)
2423	{
2424	double result = 0;
2425	int w = 0;
2426	double r = 0;
2427
2428	w = (int)Math.Round(-(2.296193e+01*s)+8.550000e+01);
2429	if( w>=85 )
2430	{
2431	r = -6.931e-01;
2432	}
2433	if( w==84 )
2434	{
2435	r = -7.276e-01;
2436	}
2437	if( w==83 )
2438	{
2439	r = -7.633e-01;
2440	}
2441	if( w==82 )
2442	{
2443	r = -8.001e-01;
2444	}
2445	if( w==81 )
2446	{
2447	r = -8.381e-01;
2448	}
2449	if( w==80 )
2450	{
2451	r = -8.774e-01;
2452	}
2453	if( w==79 )
2454	{
2455	r = -9.179e-01;
2456	}
2457	if( w==78 )
2458	{
2459	r = -9.597e-01;
2460	}
2461	if( w==77 )
2462	{
2463	r = -1.003e+00;
2464	}
2465	if( w==76 )
2466	{
2467	r = -1.047e+00;
2468	}
2469	if( w==75 )
2470	{
2471	r = -1.093e+00;
2472	}
2473	if( w==74 )
2474	{
2475	r = -1.140e+00;
2476	}
2477	if( w==73 )
2478	{
2479	r = -1.188e+00;
2480	}
2481	if( w==72 )
2482	{
2483	r = -1.238e+00;
2484	}
2485	if( w==71 )
2486	{
2487	r = -1.289e+00;
2488	}
2489	if( w==70 )
2490	{
2491	r = -1.342e+00;
2492	}
2493	if( w==69 )
2494	{
2495	r = -1.396e+00;
2496	}
2497	if( w==68 )
2498	{
2499	r = -1.452e+00;
2500	}
2501	if( w==67 )
2502	{
2503	r = -1.509e+00;
2504	}
2505	if( w==66 )
2506	{
2507	r = -1.568e+00;
2508	}
2509	if( w==65 )
2510	{
2511	r = -1.628e+00;
2512	}
2513	if( w==64 )
2514	{
2515	r = -1.690e+00;
2516	}
2517	if( w==63 )
2518	{
2519	r = -1.753e+00;
2520	}
2521	if( w==62 )
2522	{
2523	r = -1.818e+00;
2524	}
2525	if( w==61 )
2526	{
2527	r = -1.885e+00;
2528	}
2529	if( w==60 )
2530	{
2531	r = -1.953e+00;
2532	}
2533	if( w==59 )
2534	{
2535	r = -2.023e+00;
2536	}
2537	if( w==58 )
2538	{
2539	r = -2.095e+00;
2540	}
2541	if( w==57 )
2542	{
2543	r = -2.168e+00;
2544	}
2545	if( w==56 )
2546	{
2547	r = -2.244e+00;
2548	}
2549	if( w==55 )
2550	{
2551	r = -2.321e+00;
2552	}
2553	if( w==54 )
2554	{
2555	r = -2.400e+00;
2556	}
2557	if( w==53 )
2558	{
2559	r = -2.481e+00;
2560	}
2561	if( w==52 )
2562	{
2563	r = -2.564e+00;
2564	}
2565	if( w==51 )
2566	{
2567	r = -2.648e+00;
2568	}
2569	if( w==50 )
2570	{
2571	r = -2.735e+00;
2572	}
2573	if( w==49 )
2574	{
2575	r = -2.824e+00;
2576	}
2577	if( w==48 )
2578	{
2579	r = -2.915e+00;
2580	}
2581	if( w==47 )
2582	{
2583	r = -3.008e+00;
2584	}
2585	if( w==46 )
2586	{
2587	r = -3.104e+00;
2588	}
2589	if( w==45 )
2590	{
2591	r = -3.201e+00;
2592	}
2593	if( w==44 )
2594	{
2595	r = -3.301e+00;
2596	}
2597	if( w==43 )
2598	{
2599	r = -3.403e+00;
2600	}
2601	if( w==42 )
2602	{
2603	r = -3.508e+00;
2604	}
2605	if( w==41 )
2606	{
2607	r = -3.615e+00;
2608	}
2609	if( w==40 )
2610	{
2611	r = -3.724e+00;
2612	}
2613	if( w==39 )
2614	{
2615	r = -3.836e+00;
2616	}
2617	if( w==38 )
2618	{
2619	r = -3.950e+00;
2620	}
2621	if( w==37 )
2622	{
2623	r = -4.068e+00;
2624	}
2625	if( w==36 )
2626	{
2627	r = -4.188e+00;
2628	}
2629	if( w==35 )
2630	{
2631	r = -4.311e+00;
2632	}
2633	if( w==34 )
2634	{
2635	r = -4.437e+00;
2636	}
2637	if( w==33 )
2638	{
2639	r = -4.565e+00;
2640	}
2641	if( w==32 )
2642	{
2643	r = -4.698e+00;
2644	}
2645	if( w==31 )
2646	{
2647	r = -4.833e+00;
2648	}
2649	if( w==30 )
2650	{
2651	r = -4.971e+00;
2652	}
2653	if( w==29 )
2654	{
2655	r = -5.113e+00;
2656	}
2657	if( w==28 )
2658	{
2659	r = -5.258e+00;
2660	}
2661	if( w==27 )
2662	{
2663	r = -5.408e+00;
2664	}
2665	if( w==26 )
2666	{
2667	r = -5.561e+00;
2668	}
2669	if( w==25 )
2670	{
2671	r = -5.717e+00;
2672	}
2673	if( w==24 )
2674	{
2675	r = -5.878e+00;
2676	}
2677	if( w==23 )
2678	{
2679	r = -6.044e+00;
2680	}
2681	if( w==22 )
2682	{
2683	r = -6.213e+00;
2684	}
2685	if( w==21 )
2686	{
2687	r = -6.388e+00;
2688	}
2689	if( w==20 )
2690	{
2691	r = -6.569e+00;
2692	}
2693	if( w==19 )
2694	{
2695	r = -6.753e+00;
2696	}
2697	if( w==18 )
2698	{
2699	r = -6.943e+00;
2700	}
2701	if( w==17 )
2702	{
2703	r = -7.144e+00;
2704	}
2705	if( w==16 )
2706	{
2707	r = -7.347e+00;
2708	}
2709	if( w==15 )
2710	{
2711	r = -7.557e+00;
2712	}
2713	if( w==14 )
2714	{
2715	r = -7.776e+00;
2716	}
2717	if( w==13 )
2718	{
2719	r = -7.999e+00;
2720	}
2721	if( w==12 )
2722	{
2723	r = -8.228e+00;
2724	}
2725	if( w==11 )
2726	{
2727	r = -8.469e+00;
2728	}
2729	if( w==10 )
2730	{
2731	r = -8.715e+00;
2732	}
2733	if( w==9 )
2734	{
2735	r = -8.980e+00;
2736	}
2737	if( w==8 )
2738	{
2739	r = -9.258e+00;
2740	}
2741	if( w==7 )
2742	{
2743	r = -9.532e+00;
2744	}
2745	if( w==6 )
2746	{
2747	r = -9.838e+00;
2748	}
2749	if( w==5 )
2750	{
2751	r = -1.017e+01;
2752	}
2753	if( w==4 )
2754	{
2755	r = -1.053e+01;
2756	}
2757	if( w==3 )
2758	{
2759	r = -1.087e+01;
2760	}
2761	if( w==2 )
2762	{
2763	r = -1.138e+01;
2764	}
2765	if( w==1 )
2766	{
2767	r = -1.178e+01;
2768	}
2769	if( w<=0 )
2770	{
2771	r = -1.248e+01;
2772	}
2773	result = r;
2774	return result;
2775	}
2776
2777
2778	/*************************************************************************
2779	Tail(S, 19)
2780	*************************************************************************/
2781	private static double w19(double s)
2782	{
2783	double result = 0;
2784	int w = 0;
2785	double r = 0;
2786
2787	w = (int)Math.Round(-(2.484955e+01*s)+9.500000e+01);
2788	if( w>=95 )
2789	{
2790	r = -6.776e-01;
2791	}
2792	if( w==94 )
2793	{
2794	r = -7.089e-01;
2795	}
2796	if( w==93 )
2797	{
2798	r = -7.413e-01;
2799	}
2800	if( w==92 )
2801	{
2802	r = -7.747e-01;
2803	}
2804	if( w==91 )
2805	{
2806	r = -8.090e-01;
2807	}
2808	if( w==90 )
2809	{
2810	r = -8.445e-01;
2811	}
2812	if( w==89 )
2813	{
2814	r = -8.809e-01;
2815	}
2816	if( w==88 )
2817	{
2818	r = -9.185e-01;
2819	}
2820	if( w==87 )
2821	{
2822	r = -9.571e-01;
2823	}
2824	if( w==86 )
2825	{
2826	r = -9.968e-01;
2827	}
2828	if( w==85 )
2829	{
2830	r = -1.038e+00;
2831	}
2832	if( w==84 )
2833	{
2834	r = -1.080e+00;
2835	}
2836	if( w==83 )
2837	{
2838	r = -1.123e+00;
2839	}
2840	if( w==82 )
2841	{
2842	r = -1.167e+00;
2843	}
2844	if( w==81 )
2845	{
2846	r = -1.213e+00;
2847	}
2848	if( w==80 )
2849	{
2850	r = -1.259e+00;
2851	}
2852	if( w==79 )
2853	{
2854	r = -1.307e+00;
2855	}
2856	if( w==78 )
2857	{
2858	r = -1.356e+00;
2859	}
2860	if( w==77 )
2861	{
2862	r = -1.407e+00;
2863	}
2864	if( w==76 )
2865	{
2866	r = -1.458e+00;
2867	}
2868	if( w==75 )
2869	{
2870	r = -1.511e+00;
2871	}
2872	if( w==74 )
2873	{
2874	r = -1.565e+00;
2875	}
2876	if( w==73 )
2877	{
2878	r = -1.621e+00;
2879	}
2880	if( w==72 )
2881	{
2882	r = -1.678e+00;
2883	}
2884	if( w==71 )
2885	{
2886	r = -1.736e+00;
2887	}
2888	if( w==70 )
2889	{
2890	r = -1.796e+00;
2891	}
2892	if( w==69 )
2893	{
2894	r = -1.857e+00;
2895	}
2896	if( w==68 )
2897	{
2898	r = -1.919e+00;
2899	}
2900	if( w==67 )
2901	{
2902	r = -1.983e+00;
2903	}
2904	if( w==66 )
2905	{
2906	r = -2.048e+00;
2907	}
2908	if( w==65 )
2909	{
2910	r = -2.115e+00;
2911	}
2912	if( w==64 )
2913	{
2914	r = -2.183e+00;
2915	}
2916	if( w==63 )
2917	{
2918	r = -2.253e+00;
2919	}
2920	if( w==62 )
2921	{
2922	r = -2.325e+00;
2923	}
2924	if( w==61 )
2925	{
2926	r = -2.398e+00;
2927	}
2928	if( w==60 )
2929	{
2930	r = -2.472e+00;
2931	}
2932	if( w==59 )
2933	{
2934	r = -2.548e+00;
2935	}
2936	if( w==58 )
2937	{
2938	r = -2.626e+00;
2939	}
2940	if( w==57 )
2941	{
2942	r = -2.706e+00;
2943	}
2944	if( w==56 )
2945	{
2946	r = -2.787e+00;
2947	}
2948	if( w==55 )
2949	{
2950	r = -2.870e+00;
2951	}
2952	if( w==54 )
2953	{
2954	r = -2.955e+00;
2955	}
2956	if( w==53 )
2957	{
2958	r = -3.042e+00;
2959	}
2960	if( w==52 )
2961	{
2962	r = -3.130e+00;
2963	}
2964	if( w==51 )
2965	{
2966	r = -3.220e+00;
2967	}
2968	if( w==50 )
2969	{
2970	r = -3.313e+00;
2971	}
2972	if( w==49 )
2973	{
2974	r = -3.407e+00;
2975	}
2976	if( w==48 )
2977	{
2978	r = -3.503e+00;
2979	}
2980	if( w==47 )
2981	{
2982	r = -3.601e+00;
2983	}
2984	if( w==46 )
2985	{
2986	r = -3.702e+00;
2987	}
2988	if( w==45 )
2989	{
2990	r = -3.804e+00;
2991	}
2992	if( w==44 )
2993	{
2994	r = -3.909e+00;
2995	}
2996	if( w==43 )
2997	{
2998	r = -4.015e+00;
2999	}
3000	if( w==42 )
3001	{
3002	r = -4.125e+00;
3003	}
3004	if( w==41 )
3005	{
3006	r = -4.236e+00;
3007	}
3008	if( w==40 )
3009	{
3010	r = -4.350e+00;
3011	}
3012	if( w==39 )
3013	{
3014	r = -4.466e+00;
3015	}
3016	if( w==38 )
3017	{
3018	r = -4.585e+00;
3019	}
3020	if( w==37 )
3021	{
3022	r = -4.706e+00;
3023	}
3024	if( w==36 )
3025	{
3026	r = -4.830e+00;
3027	}
3028	if( w==35 )
3029	{
3030	r = -4.957e+00;
3031	}
3032	if( w==34 )
3033	{
3034	r = -5.086e+00;
3035	}
3036	if( w==33 )
3037	{
3038	r = -5.219e+00;
3039	}
3040	if( w==32 )
3041	{
3042	r = -5.355e+00;
3043	}
3044	if( w==31 )
3045	{
3046	r = -5.493e+00;
3047	}
3048	if( w==30 )
3049	{
3050	r = -5.634e+00;
3051	}
3052	if( w==29 )
3053	{
3054	r = -5.780e+00;
3055	}
3056	if( w==28 )
3057	{
3058	r = -5.928e+00;
3059	}
3060	if( w==27 )
3061	{
3062	r = -6.080e+00;
3063	}
3064	if( w==26 )
3065	{
3066	r = -6.235e+00;
3067	}
3068	if( w==25 )
3069	{
3070	r = -6.394e+00;
3071	}
3072	if( w==24 )
3073	{
3074	r = -6.558e+00;
3075	}
3076	if( w==23 )
3077	{
3078	r = -6.726e+00;
3079	}
3080	if( w==22 )
3081	{
3082	r = -6.897e+00;
3083	}
3084	if( w==21 )
3085	{
3086	r = -7.074e+00;
3087	}
3088	if( w==20 )
3089	{
3090	r = -7.256e+00;
3091	}
3092	if( w==19 )
3093	{
3094	r = -7.443e+00;
3095	}
3096	if( w==18 )
3097	{
3098	r = -7.636e+00;
3099	}
3100	if( w==17 )
3101	{
3102	r = -7.837e+00;
3103	}
3104	if( w==16 )
3105	{
3106	r = -8.040e+00;
3107	}
3108	if( w==15 )
3109	{
3110	r = -8.250e+00;
3111	}
3112	if( w==14 )
3113	{
3114	r = -8.469e+00;
3115	}
3116	if( w==13 )
3117	{
3118	r = -8.692e+00;
3119	}
3120	if( w==12 )
3121	{
3122	r = -8.921e+00;
3123	}
3124	if( w==11 )
3125	{
3126	r = -9.162e+00;
3127	}
3128	if( w==10 )
3129	{
3130	r = -9.409e+00;
3131	}
3132	if( w==9 )
3133	{
3134	r = -9.673e+00;
3135	}
3136	if( w==8 )
3137	{
3138	r = -9.951e+00;
3139	}
3140	if( w==7 )
3141	{
3142	r = -1.023e+01;
3143	}
3144	if( w==6 )
3145	{
3146	r = -1.053e+01;
3147	}
3148	if( w==5 )
3149	{
3150	r = -1.087e+01;
3151	}
3152	if( w==4 )
3153	{
3154	r = -1.122e+01;
3155	}
3156	if( w==3 )
3157	{
3158	r = -1.156e+01;
3159	}
3160	if( w==2 )
3161	{
3162	r = -1.207e+01;
3163	}
3164	if( w==1 )
3165	{
3166	r = -1.248e+01;
3167	}
3168	if( w<=0 )
3169	{
3170	r = -1.317e+01;
3171	}
3172	result = r;
3173	return result;
3174	}
3175
3176
3177	/*************************************************************************
3178	Tail(S, 20)
3179	*************************************************************************/
3180	private static double w20(double s)
3181	{
3182	double result = 0;
3183	int w = 0;
3184	double r = 0;
3185
3186	w = (int)Math.Round(-(2.678619e+01*s)+1.050000e+02);
3187	if( w>=105 )
3188	{
3189	r = -6.787e-01;
3190	}
3191	if( w==104 )
3192	{
3193	r = -7.078e-01;
3194	}
3195	if( w==103 )
3196	{
3197	r = -7.378e-01;
3198	}
3199	if( w==102 )
3200	{
3201	r = -7.686e-01;
3202	}
3203	if( w==101 )
3204	{
3205	r = -8.004e-01;
3206	}
3207	if( w==100 )
3208	{
3209	r = -8.330e-01;
3210	}
3211	if( w==99 )
3212	{
3213	r = -8.665e-01;
3214	}
3215	if( w==98 )
3216	{
3217	r = -9.010e-01;
3218	}
3219	if( w==97 )
3220	{
3221	r = -9.363e-01;
3222	}
3223	if( w==96 )
3224	{
3225	r = -9.726e-01;
3226	}
3227	if( w==95 )
3228	{
3229	r = -1.010e+00;
3230	}
3231	if( w==94 )
3232	{
3233	r = -1.048e+00;
3234	}
3235	if( w==93 )
3236	{
3237	r = -1.087e+00;
3238	}
3239	if( w==92 )
3240	{
3241	r = -1.128e+00;
3242	}
3243	if( w==91 )
3244	{
3245	r = -1.169e+00;
3246	}
3247	if( w==90 )
3248	{
3249	r = -1.211e+00;
3250	}
3251	if( w==89 )
3252	{
3253	r = -1.254e+00;
3254	}
3255	if( w==88 )
3256	{
3257	r = -1.299e+00;
3258	}
3259	if( w==87 )
3260	{
3261	r = -1.344e+00;
3262	}
3263	if( w==86 )
3264	{
3265	r = -1.390e+00;
3266	}
3267	if( w==85 )
3268	{
3269	r = -1.438e+00;
3270	}
3271	if( w==84 )
3272	{
3273	r = -1.486e+00;
3274	}
3275	if( w==83 )
3276	{
3277	r = -1.536e+00;
3278	}
3279	if( w==82 )
3280	{
3281	r = -1.587e+00;
3282	}
3283	if( w==81 )
3284	{
3285	r = -1.639e+00;
3286	}
3287	if( w==80 )
3288	{
3289	r = -1.692e+00;
3290	}
3291	if( w==79 )
3292	{
3293	r = -1.746e+00;
3294	}
3295	if( w==78 )
3296	{
3297	r = -1.802e+00;
3298	}
3299	if( w==77 )
3300	{
3301	r = -1.859e+00;
3302	}
3303	if( w==76 )
3304	{
3305	r = -1.916e+00;
3306	}
3307	if( w==75 )
3308	{
3309	r = -1.976e+00;
3310	}
3311	if( w==74 )
3312	{
3313	r = -2.036e+00;
3314	}
3315	if( w==73 )
3316	{
3317	r = -2.098e+00;
3318	}
3319	if( w==72 )
3320	{
3321	r = -2.161e+00;
3322	}
3323	if( w==71 )
3324	{
3325	r = -2.225e+00;
3326	}
3327	if( w==70 )
3328	{
3329	r = -2.290e+00;
3330	}
3331	if( w==69 )
3332	{
3333	r = -2.357e+00;
3334	}
3335	if( w==68 )
3336	{
3337	r = -2.426e+00;
3338	}
3339	if( w==67 )
3340	{
3341	r = -2.495e+00;
3342	}
3343	if( w==66 )
3344	{
3345	r = -2.566e+00;
3346	}
3347	if( w==65 )
3348	{
3349	r = -2.639e+00;
3350	}
3351	if( w==64 )
3352	{
3353	r = -2.713e+00;
3354	}
3355	if( w==63 )
3356	{
3357	r = -2.788e+00;
3358	}
3359	if( w==62 )
3360	{
3361	r = -2.865e+00;
3362	}
3363	if( w==61 )
3364	{
3365	r = -2.943e+00;
3366	}
3367	if( w==60 )
3368	{
3369	r = -3.023e+00;
3370	}
3371	if( w==59 )
3372	{
3373	r = -3.104e+00;
3374	}
3375	if( w==58 )
3376	{
3377	r = -3.187e+00;
3378	}
3379	if( w==57 )
3380	{
3381	r = -3.272e+00;
3382	}
3383	if( w==56 )
3384	{
3385	r = -3.358e+00;
3386	}
3387	if( w==55 )
3388	{
3389	r = -3.446e+00;
3390	}
3391	if( w==54 )
3392	{
3393	r = -3.536e+00;
3394	}
3395	if( w==53 )
3396	{
3397	r = -3.627e+00;
3398	}
3399	if( w==52 )
3400	{
3401	r = -3.721e+00;
3402	}
3403	if( w==51 )
3404	{
3405	r = -3.815e+00;
3406	}
3407	if( w==50 )
3408	{
3409	r = -3.912e+00;
3410	}
3411	if( w==49 )
3412	{
3413	r = -4.011e+00;
3414	}
3415	if( w==48 )
3416	{
3417	r = -4.111e+00;
3418	}
3419	if( w==47 )
3420	{
3421	r = -4.214e+00;
3422	}
3423	if( w==46 )
3424	{
3425	r = -4.318e+00;
3426	}
3427	if( w==45 )
3428	{
3429	r = -4.425e+00;
3430	}
3431	if( w==44 )
3432	{
3433	r = -4.534e+00;
3434	}
3435	if( w==43 )
3436	{
3437	r = -4.644e+00;
3438	}
3439	if( w==42 )
3440	{
3441	r = -4.757e+00;
3442	}
3443	if( w==41 )
3444	{
3445	r = -4.872e+00;
3446	}
3447	if( w==40 )
3448	{
3449	r = -4.990e+00;
3450	}
3451	if( w==39 )
3452	{
3453	r = -5.109e+00;
3454	}
3455	if( w==38 )
3456	{
3457	r = -5.232e+00;
3458	}
3459	if( w==37 )
3460	{
3461	r = -5.356e+00;
3462	}
3463	if( w==36 )
3464	{
3465	r = -5.484e+00;
3466	}
3467	if( w==35 )
3468	{
3469	r = -5.614e+00;
3470	}
3471	if( w==34 )
3472	{
3473	r = -5.746e+00;
3474	}
3475	if( w==33 )
3476	{
3477	r = -5.882e+00;
3478	}
3479	if( w==32 )
3480	{
3481	r = -6.020e+00;
3482	}
3483	if( w==31 )
3484	{
3485	r = -6.161e+00;
3486	}
3487	if( w==30 )
3488	{
3489	r = -6.305e+00;
3490	}
3491	if( w==29 )
3492	{
3493	r = -6.453e+00;
3494	}
3495	if( w==28 )
3496	{
3497	r = -6.603e+00;
3498	}
3499	if( w==27 )
3500	{
3501	r = -6.757e+00;
3502	}
3503	if( w==26 )
3504	{
3505	r = -6.915e+00;
3506	}
3507	if( w==25 )
3508	{
3509	r = -7.076e+00;
3510	}
3511	if( w==24 )
3512	{
3513	r = -7.242e+00;
3514	}
3515	if( w==23 )
3516	{
3517	r = -7.411e+00;
3518	}
3519	if( w==22 )
3520	{
3521	r = -7.584e+00;
3522	}
3523	if( w==21 )
3524	{
3525	r = -7.763e+00;
3526	}
3527	if( w==20 )
3528	{
3529	r = -7.947e+00;
3530	}
3531	if( w==19 )
3532	{
3533	r = -8.136e+00;
3534	}
3535	if( w==18 )
3536	{
3537	r = -8.330e+00;
3538	}
3539	if( w==17 )
3540	{
3541	r = -8.530e+00;
3542	}
3543	if( w==16 )
3544	{
3545	r = -8.733e+00;
3546	}
3547	if( w==15 )
3548	{
3549	r = -8.943e+00;
3550	}
3551	if( w==14 )
3552	{
3553	r = -9.162e+00;
3554	}
3555	if( w==13 )
3556	{
3557	r = -9.386e+00;
3558	}
3559	if( w==12 )
3560	{
3561	r = -9.614e+00;
3562	}
3563	if( w==11 )
3564	{
3565	r = -9.856e+00;
3566	}
3567	if( w==10 )
3568	{
3569	r = -1.010e+01;
3570	}
3571	if( w==9 )
3572	{
3573	r = -1.037e+01;
3574	}
3575	if( w==8 )
3576	{
3577	r = -1.064e+01;
3578	}
3579	if( w==7 )
3580	{
3581	r = -1.092e+01;
3582	}
3583	if( w==6 )
3584	{
3585	r = -1.122e+01;
3586	}
3587	if( w==5 )
3588	{
3589	r = -1.156e+01;
3590	}
3591	if( w==4 )
3592	{
3593	r = -1.192e+01;
3594	}
3595	if( w==3 )
3596	{
3597	r = -1.225e+01;
3598	}
3599	if( w==2 )
3600	{
3601	r = -1.276e+01;
3602	}
3603	if( w==1 )
3604	{
3605	r = -1.317e+01;
3606	}
3607	if( w<=0 )
3608	{
3609	r = -1.386e+01;
3610	}
3611	result = r;
3612	return result;
3613	}
3614
3615
3616	/*************************************************************************
3617	Tail(S, 21)
3618	*************************************************************************/
3619	private static double w21(double s)
3620	{
3621	double result = 0;
3622	int w = 0;
3623	double r = 0;
3624
3625	w = (int)Math.Round(-(2.877064e+01*s)+1.155000e+02);
3626	if( w>=115 )
3627	{
3628	r = -6.931e-01;
3629	}
3630	if( w==114 )
3631	{
3632	r = -7.207e-01;
3633	}
3634	if( w==113 )
3635	{
3636	r = -7.489e-01;
3637	}
3638	if( w==112 )
3639	{
3640	r = -7.779e-01;
3641	}
3642	if( w==111 )
3643	{
3644	r = -8.077e-01;
3645	}
3646	if( w==110 )
3647	{
3648	r = -8.383e-01;
3649	}
3650	if( w==109 )
3651	{
3652	r = -8.697e-01;
3653	}
3654	if( w==108 )
3655	{
3656	r = -9.018e-01;
3657	}
3658	if( w==107 )
3659	{
3660	r = -9.348e-01;
3661	}
3662	if( w==106 )
3663	{
3664	r = -9.685e-01;
3665	}
3666	if( w==105 )
3667	{
3668	r = -1.003e+00;
3669	}
3670	if( w==104 )
3671	{
3672	r = -1.039e+00;
3673	}
3674	if( w==103 )
3675	{
3676	r = -1.075e+00;
3677	}
3678	if( w==102 )
3679	{
3680	r = -1.112e+00;
3681	}
3682	if( w==101 )
3683	{
3684	r = -1.150e+00;
3685	}
3686	if( w==100 )
3687	{
3688	r = -1.189e+00;
3689	}
3690	if( w==99 )
3691	{
3692	r = -1.229e+00;
3693	}
3694	if( w==98 )
3695	{
3696	r = -1.269e+00;
3697	}
3698	if( w==97 )
3699	{
3700	r = -1.311e+00;
3701	}
3702	if( w==96 )
3703	{
3704	r = -1.353e+00;
3705	}
3706	if( w==95 )
3707	{
3708	r = -1.397e+00;
3709	}
3710	if( w==94 )
3711	{
3712	r = -1.441e+00;
3713	}
3714	if( w==93 )
3715	{
3716	r = -1.486e+00;
3717	}
3718	if( w==92 )
3719	{
3720	r = -1.533e+00;
3721	}
3722	if( w==91 )
3723	{
3724	r = -1.580e+00;
3725	}
3726	if( w==90 )
3727	{
3728	r = -1.628e+00;
3729	}
3730	if( w==89 )
3731	{
3732	r = -1.677e+00;
3733	}
3734	if( w==88 )
3735	{
3736	r = -1.728e+00;
3737	}
3738	if( w==87 )
3739	{
3740	r = -1.779e+00;
3741	}
3742	if( w==86 )
3743	{
3744	r = -1.831e+00;
3745	}
3746	if( w==85 )
3747	{
3748	r = -1.884e+00;
3749	}
3750	if( w==84 )
3751	{
3752	r = -1.939e+00;
3753	}
3754	if( w==83 )
3755	{
3756	r = -1.994e+00;
3757	}
3758	if( w==82 )
3759	{
3760	r = -2.051e+00;
3761	}
3762	if( w==81 )
3763	{
3764	r = -2.108e+00;
3765	}
3766	if( w==80 )
3767	{
3768	r = -2.167e+00;
3769	}
3770	if( w==79 )
3771	{
3772	r = -2.227e+00;
3773	}
3774	if( w==78 )
3775	{
3776	r = -2.288e+00;
3777	}
3778	if( w==77 )
3779	{
3780	r = -2.350e+00;
3781	}
3782	if( w==76 )
3783	{
3784	r = -2.414e+00;
3785	}
3786	if( w==75 )
3787	{
3788	r = -2.478e+00;
3789	}
3790	if( w==74 )
3791	{
3792	r = -2.544e+00;
3793	}
3794	if( w==73 )
3795	{
3796	r = -2.611e+00;
3797	}
3798	if( w==72 )
3799	{
3800	r = -2.679e+00;
3801	}
3802	if( w==71 )
3803	{
3804	r = -2.748e+00;
3805	}
3806	if( w==70 )
3807	{
3808	r = -2.819e+00;
3809	}
3810	if( w==69 )
3811	{
3812	r = -2.891e+00;
3813	}
3814	if( w==68 )
3815	{
3816	r = -2.964e+00;
3817	}
3818	if( w==67 )
3819	{
3820	r = -3.039e+00;
3821	}
3822	if( w==66 )
3823	{
3824	r = -3.115e+00;
3825	}
3826	if( w==65 )
3827	{
3828	r = -3.192e+00;
3829	}
3830	if( w==64 )
3831	{
3832	r = -3.270e+00;
3833	}
3834	if( w==63 )
3835	{
3836	r = -3.350e+00;
3837	}
3838	if( w==62 )
3839	{
3840	r = -3.432e+00;
3841	}
3842	if( w==61 )
3843	{
3844	r = -3.515e+00;
3845	}
3846	if( w==60 )
3847	{
3848	r = -3.599e+00;
3849	}
3850	if( w==59 )
3851	{
3852	r = -3.685e+00;
3853	}
3854	if( w==58 )
3855	{
3856	r = -3.772e+00;
3857	}
3858	if( w==57 )
3859	{
3860	r = -3.861e+00;
3861	}
3862	if( w==56 )
3863	{
3864	r = -3.952e+00;
3865	}
3866	if( w==55 )
3867	{
3868	r = -4.044e+00;
3869	}
3870	if( w==54 )
3871	{
3872	r = -4.138e+00;
3873	}
3874	if( w==53 )
3875	{
3876	r = -4.233e+00;
3877	}
3878	if( w==52 )
3879	{
3880	r = -4.330e+00;
3881	}
3882	if( w==51 )
3883	{
3884	r = -4.429e+00;
3885	}
3886	if( w==50 )
3887	{
3888	r = -4.530e+00;
3889	}
3890	if( w==49 )
3891	{
3892	r = -4.632e+00;
3893	}
3894	if( w==48 )
3895	{
3896	r = -4.736e+00;
3897	}
3898	if( w==47 )
3899	{
3900	r = -4.842e+00;
3901	}
3902	if( w==46 )
3903	{
3904	r = -4.950e+00;
3905	}
3906	if( w==45 )
3907	{
3908	r = -5.060e+00;
3909	}
3910	if( w==44 )
3911	{
3912	r = -5.172e+00;
3913	}
3914	if( w==43 )
3915	{
3916	r = -5.286e+00;
3917	}
3918	if( w==42 )
3919	{
3920	r = -5.402e+00;
3921	}
3922	if( w==41 )
3923	{
3924	r = -5.520e+00;
3925	}
3926	if( w==40 )
3927	{
3928	r = -5.641e+00;
3929	}
3930	if( w==39 )
3931	{
3932	r = -5.763e+00;
3933	}
3934	if( w==38 )
3935	{
3936	r = -5.889e+00;
3937	}
3938	if( w==37 )
3939	{
3940	r = -6.016e+00;
3941	}
3942	if( w==36 )
3943	{
3944	r = -6.146e+00;
3945	}
3946	if( w==35 )
3947	{
3948	r = -6.278e+00;
3949	}
3950	if( w==34 )
3951	{
3952	r = -6.413e+00;
3953	}
3954	if( w==33 )
3955	{
3956	r = -6.551e+00;
3957	}
3958	if( w==32 )
3959	{
3960	r = -6.692e+00;
3961	}
3962	if( w==31 )
3963	{
3964	r = -6.835e+00;
3965	}
3966	if( w==30 )
3967	{
3968	r = -6.981e+00;
3969	}
3970	if( w==29 )
3971	{
3972	r = -7.131e+00;
3973	}
3974	if( w==28 )
3975	{
3976	r = -7.283e+00;
3977	}
3978	if( w==27 )
3979	{
3980	r = -7.439e+00;
3981	}
3982	if( w==26 )
3983	{
3984	r = -7.599e+00;
3985	}
3986	if( w==25 )
3987	{
3988	r = -7.762e+00;
3989	}
3990	if( w==24 )
3991	{
3992	r = -7.928e+00;
3993	}
3994	if( w==23 )
3995	{
3996	r = -8.099e+00;
3997	}
3998	if( w==22 )
3999	{
4000	r = -8.274e+00;
4001	}
4002	if( w==21 )
4003	{
4004	r = -8.454e+00;
4005	}
4006	if( w==20 )
4007	{
4008	r = -8.640e+00;
4009	}
4010	if( w==19 )
4011	{
4012	r = -8.829e+00;
4013	}
4014	if( w==18 )
4015	{
4016	r = -9.023e+00;
4017	}
4018	if( w==17 )
4019	{
4020	r = -9.223e+00;
4021	}
4022	if( w==16 )
4023	{
4024	r = -9.426e+00;
4025	}
4026	if( w==15 )
4027	{
4028	r = -9.636e+00;
4029	}
4030	if( w==14 )
4031	{
4032	r = -9.856e+00;
4033	}
4034	if( w==13 )
4035	{
4036	r = -1.008e+01;
4037	}
4038	if( w==12 )
4039	{
4040	r = -1.031e+01;
4041	}
4042	if( w==11 )
4043	{
4044	r = -1.055e+01;
4045	}
4046	if( w==10 )
4047	{
4048	r = -1.079e+01;
4049	}
4050	if( w==9 )
4051	{
4052	r = -1.106e+01;
4053	}
4054	if( w==8 )
4055	{
4056	r = -1.134e+01;
4057	}
4058	if( w==7 )
4059	{
4060	r = -1.161e+01;
4061	}
4062	if( w==6 )
4063	{
4064	r = -1.192e+01;
4065	}
4066	if( w==5 )
4067	{
4068	r = -1.225e+01;
4069	}
4070	if( w==4 )
4071	{
4072	r = -1.261e+01;
4073	}
4074	if( w==3 )
4075	{
4076	r = -1.295e+01;
4077	}
4078	if( w==2 )
4079	{
4080	r = -1.346e+01;
4081	}
4082	if( w==1 )
4083	{
4084	r = -1.386e+01;
4085	}
4086	if( w<=0 )
4087	{
4088	r = -1.456e+01;
4089	}
4090	result = r;
4091	return result;
4092	}
4093
4094
4095	/*************************************************************************
4096	Tail(S, 22)
4097	*************************************************************************/
4098	private static double w22(double s)
4099	{
4100	double result = 0;
4101	int w = 0;
4102	double r = 0;
4103
4104	w = (int)Math.Round(-(3.080179e+01*s)+1.265000e+02);
4105	if( w>=126 )
4106	{
4107	r = -6.931e-01;
4108	}
4109	if( w==125 )
4110	{
4111	r = -7.189e-01;
4112	}
4113	if( w==124 )
4114	{
4115	r = -7.452e-01;
4116	}
4117	if( w==123 )
4118	{
4119	r = -7.722e-01;
4120	}
4121	if( w==122 )
4122	{
4123	r = -7.999e-01;
4124	}
4125	if( w==121 )
4126	{
4127	r = -8.283e-01;
4128	}
4129	if( w==120 )
4130	{
4131	r = -8.573e-01;
4132	}
4133	if( w==119 )
4134	{
4135	r = -8.871e-01;
4136	}
4137	if( w==118 )
4138	{
4139	r = -9.175e-01;
4140	}
4141	if( w==117 )
4142	{
4143	r = -9.486e-01;
4144	}
4145	if( w==116 )
4146	{
4147	r = -9.805e-01;
4148	}
4149	if( w==115 )
4150	{
4151	r = -1.013e+00;
4152	}
4153	if( w==114 )
4154	{
4155	r = -1.046e+00;
4156	}
4157	if( w==113 )
4158	{
4159	r = -1.080e+00;
4160	}
4161	if( w==112 )
4162	{
4163	r = -1.115e+00;
4164	}
4165	if( w==111 )
4166	{
4167	r = -1.151e+00;
4168	}
4169	if( w==110 )
4170	{
4171	r = -1.187e+00;
4172	}
4173	if( w==109 )
4174	{
4175	r = -1.224e+00;
4176	}
4177	if( w==108 )
4178	{
4179	r = -1.262e+00;
4180	}
4181	if( w==107 )
4182	{
4183	r = -1.301e+00;
4184	}
4185	if( w==106 )
4186	{
4187	r = -1.340e+00;
4188	}
4189	if( w==105 )
4190	{
4191	r = -1.381e+00;
4192	}
4193	if( w==104 )
4194	{
4195	r = -1.422e+00;
4196	}
4197	if( w==103 )
4198	{
4199	r = -1.464e+00;
4200	}
4201	if( w==102 )
4202	{
4203	r = -1.506e+00;
4204	}
4205	if( w==101 )
4206	{
4207	r = -1.550e+00;
4208	}
4209	if( w==100 )
4210	{
4211	r = -1.594e+00;
4212	}
4213	if( w==99 )
4214	{
4215	r = -1.640e+00;
4216	}
4217	if( w==98 )
4218	{
4219	r = -1.686e+00;
4220	}
4221	if( w==97 )
4222	{
4223	r = -1.733e+00;
4224	}
4225	if( w==96 )
4226	{
4227	r = -1.781e+00;
4228	}
4229	if( w==95 )
4230	{
4231	r = -1.830e+00;
4232	}
4233	if( w==94 )
4234	{
4235	r = -1.880e+00;
4236	}
4237	if( w==93 )
4238	{
4239	r = -1.930e+00;
4240	}
4241	if( w==92 )
4242	{
4243	r = -1.982e+00;
4244	}
4245	if( w==91 )
4246	{
4247	r = -2.034e+00;
4248	}
4249	if( w==90 )
4250	{
4251	r = -2.088e+00;
4252	}
4253	if( w==89 )
4254	{
4255	r = -2.142e+00;
4256	}
4257	if( w==88 )
4258	{
4259	r = -2.198e+00;
4260	}
4261	if( w==87 )
4262	{
4263	r = -2.254e+00;
4264	}
4265	if( w==86 )
4266	{
4267	r = -2.312e+00;
4268	}
4269	if( w==85 )
4270	{
4271	r = -2.370e+00;
4272	}
4273	if( w==84 )
4274	{
4275	r = -2.429e+00;
4276	}
4277	if( w==83 )
4278	{
4279	r = -2.490e+00;
4280	}
4281	if( w==82 )
4282	{
4283	r = -2.551e+00;
4284	}
4285	if( w==81 )
4286	{
4287	r = -2.614e+00;
4288	}
4289	if( w==80 )
4290	{
4291	r = -2.677e+00;
4292	}
4293	if( w==79 )
4294	{
4295	r = -2.742e+00;
4296	}
4297	if( w==78 )
4298	{
4299	r = -2.808e+00;
4300	}
4301	if( w==77 )
4302	{
4303	r = -2.875e+00;
4304	}
4305	if( w==76 )
4306	{
4307	r = -2.943e+00;
4308	}
4309	if( w==75 )
4310	{
4311	r = -3.012e+00;
4312	}
4313	if( w==74 )
4314	{
4315	r = -3.082e+00;
4316	}
4317	if( w==73 )
4318	{
4319	r = -3.153e+00;
4320	}
4321	if( w==72 )
4322	{
4323	r = -3.226e+00;
4324	}
4325	if( w==71 )
4326	{
4327	r = -3.300e+00;
4328	}
4329	if( w==70 )
4330	{
4331	r = -3.375e+00;
4332	}
4333	if( w==69 )
4334	{
4335	r = -3.451e+00;
4336	}
4337	if( w==68 )
4338	{
4339	r = -3.529e+00;
4340	}
4341	if( w==67 )
4342	{
4343	r = -3.607e+00;
4344	}
4345	if( w==66 )
4346	{
4347	r = -3.687e+00;
4348	}
4349	if( w==65 )
4350	{
4351	r = -3.769e+00;
4352	}
4353	if( w==64 )
4354	{
4355	r = -3.851e+00;
4356	}
4357	if( w==63 )
4358	{
4359	r = -3.935e+00;
4360	}
4361	if( w==62 )
4362	{
4363	r = -4.021e+00;
4364	}
4365	if( w==61 )
4366	{
4367	r = -4.108e+00;
4368	}
4369	if( w==60 )
4370	{
4371	r = -4.196e+00;
4372	}
4373	if( w==59 )
4374	{
4375	r = -4.285e+00;
4376	}
4377	if( w==58 )
4378	{
4379	r = -4.376e+00;
4380	}
4381	if( w==57 )
4382	{
4383	r = -4.469e+00;
4384	}
4385	if( w==56 )
4386	{
4387	r = -4.563e+00;
4388	}
4389	if( w==55 )
4390	{
4391	r = -4.659e+00;
4392	}
4393	if( w==54 )
4394	{
4395	r = -4.756e+00;
4396	}
4397	if( w==53 )
4398	{
4399	r = -4.855e+00;
4400	}
4401	if( w==52 )
4402	{
4403	r = -4.955e+00;
4404	}
4405	if( w==51 )
4406	{
4407	r = -5.057e+00;
4408	}
4409	if( w==50 )
4410	{
4411	r = -5.161e+00;
4412	}
4413	if( w==49 )
4414	{
4415	r = -5.266e+00;
4416	}
4417	if( w==48 )
4418	{
4419	r = -5.374e+00;
4420	}
4421	if( w==47 )
4422	{
4423	r = -5.483e+00;
4424	}
4425	if( w==46 )
4426	{
4427	r = -5.594e+00;
4428	}
4429	if( w==45 )
4430	{
4431	r = -5.706e+00;
4432	}
4433	if( w==44 )
4434	{
4435	r = -5.821e+00;
4436	}
4437	if( w==43 )
4438	{
4439	r = -5.938e+00;
4440	}
4441	if( w==42 )
4442	{
4443	r = -6.057e+00;
4444	}
4445	if( w==41 )
4446	{
4447	r = -6.177e+00;
4448	}
4449	if( w==40 )
4450	{
4451	r = -6.300e+00;
4452	}
4453	if( w==39 )
4454	{
4455	r = -6.426e+00;
4456	}
4457	if( w==38 )
4458	{
4459	r = -6.553e+00;
4460	}
4461	if( w==37 )
4462	{
4463	r = -6.683e+00;
4464	}
4465	if( w==36 )
4466	{
4467	r = -6.815e+00;
4468	}
4469	if( w==35 )
4470	{
4471	r = -6.949e+00;
4472	}
4473	if( w==34 )
4474	{
4475	r = -7.086e+00;
4476	}
4477	if( w==33 )
4478	{
4479	r = -7.226e+00;
4480	}
4481	if( w==32 )
4482	{
4483	r = -7.368e+00;
4484	}
4485	if( w==31 )
4486	{
4487	r = -7.513e+00;
4488	}
4489	if( w==30 )
4490	{
4491	r = -7.661e+00;
4492	}
4493	if( w==29 )
4494	{
4495	r = -7.813e+00;
4496	}
4497	if( w==28 )
4498	{
4499	r = -7.966e+00;
4500	}
4501	if( w==27 )
4502	{
4503	r = -8.124e+00;
4504	}
4505	if( w==26 )
4506	{
4507	r = -8.285e+00;
4508	}
4509	if( w==25 )
4510	{
4511	r = -8.449e+00;
4512	}
4513	if( w==24 )
4514	{
4515	r = -8.617e+00;
4516	}
4517	if( w==23 )
4518	{
4519	r = -8.789e+00;
4520	}
4521	if( w==22 )
4522	{
4523	r = -8.965e+00;
4524	}
4525	if( w==21 )
4526	{
4527	r = -9.147e+00;
4528	}
4529	if( w==20 )
4530	{
4531	r = -9.333e+00;
4532	}
4533	if( w==19 )
4534	{
4535	r = -9.522e+00;
4536	}
4537	if( w==18 )
4538	{
4539	r = -9.716e+00;
4540	}
4541	if( w==17 )
4542	{
4543	r = -9.917e+00;
4544	}
4545	if( w==16 )
4546	{
4547	r = -1.012e+01;
4548	}
4549	if( w==15 )
4550	{
4551	r = -1.033e+01;
4552	}
4553	if( w==14 )
4554	{
4555	r = -1.055e+01;
4556	}
4557	if( w==13 )
4558	{
4559	r = -1.077e+01;
4560	}
4561	if( w==12 )
4562	{
4563	r = -1.100e+01;
4564	}
4565	if( w==11 )
4566	{
4567	r = -1.124e+01;
4568	}
4569	if( w==10 )
4570	{
4571	r = -1.149e+01;
4572	}
4573	if( w==9 )
4574	{
4575	r = -1.175e+01;
4576	}
4577	if( w==8 )
4578	{
4579	r = -1.203e+01;
4580	}
4581	if( w==7 )
4582	{
4583	r = -1.230e+01;
4584	}
4585	if( w==6 )
4586	{
4587	r = -1.261e+01;
4588	}
4589	if( w==5 )
4590	{
4591	r = -1.295e+01;
4592	}
4593	if( w==4 )
4594	{
4595	r = -1.330e+01;
4596	}
4597	if( w==3 )
4598	{
4599	r = -1.364e+01;
4600	}
4601	if( w==2 )
4602	{
4603	r = -1.415e+01;
4604	}
4605	if( w==1 )
4606	{
4607	r = -1.456e+01;
4608	}
4609	if( w<=0 )
4610	{
4611	r = -1.525e+01;
4612	}
4613	result = r;
4614	return result;
4615	}
4616
4617
4618	/*************************************************************************
4619	Tail(S, 23)
4620	*************************************************************************/
4621	private static double w23(double s)
4622	{
4623	double result = 0;
4624	int w = 0;
4625	double r = 0;
4626
4627	w = (int)Math.Round(-(3.287856e+01*s)+1.380000e+02);
4628	if( w>=138 )
4629	{
4630	r = -6.813e-01;
4631	}
4632	if( w==137 )
4633	{
4634	r = -7.051e-01;
4635	}
4636	if( w==136 )
4637	{
4638	r = -7.295e-01;
4639	}
4640	if( w==135 )
4641	{
4642	r = -7.544e-01;
4643	}
4644	if( w==134 )
4645	{
4646	r = -7.800e-01;
4647	}
4648	if( w==133 )
4649	{
4650	r = -8.061e-01;
4651	}
4652	if( w==132 )
4653	{
4654	r = -8.328e-01;
4655	}
4656	if( w==131 )
4657	{
4658	r = -8.601e-01;
4659	}
4660	if( w==130 )
4661	{
4662	r = -8.880e-01;
4663	}
4664	if( w==129 )
4665	{
4666	r = -9.166e-01;
4667	}
4668	if( w==128 )
4669	{
4670	r = -9.457e-01;
4671	}
4672	if( w==127 )
4673	{
4674	r = -9.755e-01;
4675	}
4676	if( w==126 )
4677	{
4678	r = -1.006e+00;
4679	}
4680	if( w==125 )
4681	{
4682	r = -1.037e+00;
4683	}
4684	if( w==124 )
4685	{
4686	r = -1.069e+00;
4687	}
4688	if( w==123 )
4689	{
4690	r = -1.101e+00;
4691	}
4692	if( w==122 )
4693	{
4694	r = -1.134e+00;
4695	}
4696	if( w==121 )
4697	{
4698	r = -1.168e+00;
4699	}
4700	if( w==120 )
4701	{
4702	r = -1.202e+00;
4703	}
4704	if( w==119 )
4705	{
4706	r = -1.237e+00;
4707	}
4708	if( w==118 )
4709	{
4710	r = -1.273e+00;
4711	}
4712	if( w==117 )
4713	{
4714	r = -1.309e+00;
4715	}
4716	if( w==116 )
4717	{
4718	r = -1.347e+00;
4719	}
4720	if( w==115 )
4721	{
4722	r = -1.384e+00;
4723	}
4724	if( w==114 )
4725	{
4726	r = -1.423e+00;
4727	}
4728	if( w==113 )
4729	{
4730	r = -1.462e+00;
4731	}
4732	if( w==112 )
4733	{
4734	r = -1.502e+00;
4735	}
4736	if( w==111 )
4737	{
4738	r = -1.543e+00;
4739	}
4740	if( w==110 )
4741	{
4742	r = -1.585e+00;
4743	}
4744	if( w==109 )
4745	{
4746	r = -1.627e+00;
4747	}
4748	if( w==108 )
4749	{
4750	r = -1.670e+00;
4751	}
4752	if( w==107 )
4753	{
4754	r = -1.714e+00;
4755	}
4756	if( w==106 )
4757	{
4758	r = -1.758e+00;
4759	}
4760	if( w==105 )
4761	{
4762	r = -1.804e+00;
4763	}
4764	if( w==104 )
4765	{
4766	r = -1.850e+00;
4767	}
4768	if( w==103 )
4769	{
4770	r = -1.897e+00;
4771	}
4772	if( w==102 )
4773	{
4774	r = -1.944e+00;
4775	}
4776	if( w==101 )
4777	{
4778	r = -1.993e+00;
4779	}
4780	if( w==100 )
4781	{
4782	r = -2.042e+00;
4783	}
4784	if( w==99 )
4785	{
4786	r = -2.093e+00;
4787	}
4788	if( w==98 )
4789	{
4790	r = -2.144e+00;
4791	}
4792	if( w==97 )
4793	{
4794	r = -2.195e+00;
4795	}
4796	if( w==96 )
4797	{
4798	r = -2.248e+00;
4799	}
4800	if( w==95 )
4801	{
4802	r = -2.302e+00;
4803	}
4804	if( w==94 )
4805	{
4806	r = -2.356e+00;
4807	}
4808	if( w==93 )
4809	{
4810	r = -2.412e+00;
4811	}
4812	if( w==92 )
4813	{
4814	r = -2.468e+00;
4815	}
4816	if( w==91 )
4817	{
4818	r = -2.525e+00;
4819	}
4820	if( w==90 )
4821	{
4822	r = -2.583e+00;
4823	}
4824	if( w==89 )
4825	{
4826	r = -2.642e+00;
4827	}
4828	if( w==88 )
4829	{
4830	r = -2.702e+00;
4831	}
4832	if( w==87 )
4833	{
4834	r = -2.763e+00;
4835	}
4836	if( w==86 )
4837	{
4838	r = -2.825e+00;
4839	}
4840	if( w==85 )
4841	{
4842	r = -2.888e+00;
4843	}
4844	if( w==84 )
4845	{
4846	r = -2.951e+00;
4847	}
4848	if( w==83 )
4849	{
4850	r = -3.016e+00;
4851	}
4852	if( w==82 )
4853	{
4854	r = -3.082e+00;
4855	}
4856	if( w==81 )
4857	{
4858	r = -3.149e+00;
4859	}
4860	if( w==80 )
4861	{
4862	r = -3.216e+00;
4863	}
4864	if( w==79 )
4865	{
4866	r = -3.285e+00;
4867	}
4868	if( w==78 )
4869	{
4870	r = -3.355e+00;
4871	}
4872	if( w==77 )
4873	{
4874	r = -3.426e+00;
4875	}
4876	if( w==76 )
4877	{
4878	r = -3.498e+00;
4879	}
4880	if( w==75 )
4881	{
4882	r = -3.571e+00;
4883	}
4884	if( w==74 )
4885	{
4886	r = -3.645e+00;
4887	}
4888	if( w==73 )
4889	{
4890	r = -3.721e+00;
4891	}
4892	if( w==72 )
4893	{
4894	r = -3.797e+00;
4895	}
4896	if( w==71 )
4897	{
4898	r = -3.875e+00;
4899	}
4900	if( w==70 )
4901	{
4902	r = -3.953e+00;
4903	}
4904	if( w==69 )
4905	{
4906	r = -4.033e+00;
4907	}
4908	if( w==68 )
4909	{
4910	r = -4.114e+00;
4911	}
4912	if( w==67 )
4913	{
4914	r = -4.197e+00;
4915	}
4916	if( w==66 )
4917	{
4918	r = -4.280e+00;
4919	}
4920	if( w==65 )
4921	{
4922	r = -4.365e+00;
4923	}
4924	if( w==64 )
4925	{
4926	r = -4.451e+00;
4927	}
4928	if( w==63 )
4929	{
4930	r = -4.539e+00;
4931	}
4932	if( w==62 )
4933	{
4934	r = -4.628e+00;
4935	}
4936	if( w==61 )
4937	{
4938	r = -4.718e+00;
4939	}
4940	if( w==60 )
4941	{
4942	r = -4.809e+00;
4943	}
4944	if( w==59 )
4945	{
4946	r = -4.902e+00;
4947	}
4948	if( w==58 )
4949	{
4950	r = -4.996e+00;
4951	}
4952	if( w==57 )
4953	{
4954	r = -5.092e+00;
4955	}
4956	if( w==56 )
4957	{
4958	r = -5.189e+00;
4959	}
4960	if( w==55 )
4961	{
4962	r = -5.287e+00;
4963	}
4964	if( w==54 )
4965	{
4966	r = -5.388e+00;
4967	}
4968	if( w==53 )
4969	{
4970	r = -5.489e+00;
4971	}
4972	if( w==52 )
4973	{
4974	r = -5.592e+00;
4975	}
4976	if( w==51 )
4977	{
4978	r = -5.697e+00;
4979	}
4980	if( w==50 )
4981	{
4982	r = -5.804e+00;
4983	}
4984	if( w==49 )
4985	{
4986	r = -5.912e+00;
4987	}
4988	if( w==48 )
4989	{
4990	r = -6.022e+00;
4991	}
4992	if( w==47 )
4993	{
4994	r = -6.133e+00;
4995	}
4996	if( w==46 )
4997	{
4998	r = -6.247e+00;
4999	}
5000	if( w==45 )
5001	{
5002	r = -6.362e+00;
5003	}
5004	if( w==44 )
5005	{
5006	r = -6.479e+00;
5007	}
5008	if( w==43 )
5009	{
5010	r = -6.598e+00;
5011	}
5012	if( w==42 )
5013	{
5014	r = -6.719e+00;
5015	}
5016	if( w==41 )
5017	{
5018	r = -6.842e+00;
5019	}
5020	if( w==40 )
5021	{
5022	r = -6.967e+00;
5023	}
5024	if( w==39 )
5025	{
5026	r = -7.094e+00;
5027	}
5028	if( w==38 )
5029	{
5030	r = -7.224e+00;
5031	}
5032	if( w==37 )
5033	{
5034	r = -7.355e+00;
5035	}
5036	if( w==36 )
5037	{
5038	r = -7.489e+00;
5039	}
5040	if( w==35 )
5041	{
5042	r = -7.625e+00;
5043	}
5044	if( w==34 )
5045	{
5046	r = -7.764e+00;
5047	}
5048	if( w==33 )
5049	{
5050	r = -7.905e+00;
5051	}
5052	if( w==32 )
5053	{
5054	r = -8.049e+00;
5055	}
5056	if( w==31 )
5057	{
5058	r = -8.196e+00;
5059	}
5060	if( w==30 )
5061	{
5062	r = -8.345e+00;
5063	}
5064	if( w==29 )
5065	{
5066	r = -8.498e+00;
5067	}
5068	if( w==28 )
5069	{
5070	r = -8.653e+00;
5071	}
5072	if( w==27 )
5073	{
5074	r = -8.811e+00;
5075	}
5076	if( w==26 )
5077	{
5078	r = -8.974e+00;
5079	}
5080	if( w==25 )
5081	{
5082	r = -9.139e+00;
5083	}
5084	if( w==24 )
5085	{
5086	r = -9.308e+00;
5087	}
5088	if( w==23 )
5089	{
5090	r = -9.481e+00;
5091	}
5092	if( w==22 )
5093	{
5094	r = -9.658e+00;
5095	}
5096	if( w==21 )
5097	{
5098	r = -9.840e+00;
5099	}
5100	if( w==20 )
5101	{
5102	r = -1.003e+01;
5103	}
5104	if( w==19 )
5105	{
5106	r = -1.022e+01;
5107	}
5108	if( w==18 )
5109	{
5110	r = -1.041e+01;
5111	}
5112	if( w==17 )
5113	{
5114	r = -1.061e+01;
5115	}
5116	if( w==16 )
5117	{
5118	r = -1.081e+01;
5119	}
5120	if( w==15 )
5121	{
5122	r = -1.102e+01;
5123	}
5124	if( w==14 )
5125	{
5126	r = -1.124e+01;
5127	}
5128	if( w==13 )
5129	{
5130	r = -1.147e+01;
5131	}
5132	if( w==12 )
5133	{
5134	r = -1.169e+01;
5135	}
5136	if( w==11 )
5137	{
5138	r = -1.194e+01;
5139	}
5140	if( w==10 )
5141	{
5142	r = -1.218e+01;
5143	}
5144	if( w==9 )
5145	{
5146	r = -1.245e+01;
5147	}
5148	if( w==8 )
5149	{
5150	r = -1.272e+01;
5151	}
5152	if( w==7 )
5153	{
5154	r = -1.300e+01;
5155	}
5156	if( w==6 )
5157	{
5158	r = -1.330e+01;
5159	}
5160	if( w==5 )
5161	{
5162	r = -1.364e+01;
5163	}
5164	if( w==4 )
5165	{
5166	r = -1.400e+01;
5167	}
5168	if( w==3 )
5169	{
5170	r = -1.433e+01;
5171	}
5172	if( w==2 )
5173	{
5174	r = -1.484e+01;
5175	}
5176	if( w==1 )
5177	{
5178	r = -1.525e+01;
5179	}
5180	if( w<=0 )
5181	{
5182	r = -1.594e+01;
5183	}
5184	result = r;
5185	return result;
5186	}
5187
5188
5189	/*************************************************************************
5190	Tail(S, 24)
5191	*************************************************************************/
5192	private static double w24(double s)
5193	{
5194	double result = 0;
5195	int w = 0;
5196	double r = 0;
5197
5198	w = (int)Math.Round(-(3.500000e+01*s)+1.500000e+02);
5199	if( w>=150 )
5200	{
5201	r = -6.820e-01;
5202	}
5203	if( w==149 )
5204	{
5205	r = -7.044e-01;
5206	}
5207	if( w==148 )
5208	{
5209	r = -7.273e-01;
5210	}
5211	if( w==147 )
5212	{
5213	r = -7.507e-01;
5214	}
5215	if( w==146 )
5216	{
5217	r = -7.746e-01;
5218	}
5219	if( w==145 )
5220	{
5221	r = -7.990e-01;
5222	}
5223	if( w==144 )
5224	{
5225	r = -8.239e-01;
5226	}
5227	if( w==143 )
5228	{
5229	r = -8.494e-01;
5230	}
5231	if( w==142 )
5232	{
5233	r = -8.754e-01;
5234	}
5235	if( w==141 )
5236	{
5237	r = -9.020e-01;
5238	}
5239	if( w==140 )
5240	{
5241	r = -9.291e-01;
5242	}
5243	if( w==139 )
5244	{
5245	r = -9.567e-01;
5246	}
5247	if( w==138 )
5248	{
5249	r = -9.849e-01;
5250	}
5251	if( w==137 )
5252	{
5253	r = -1.014e+00;
5254	}
5255	if( w==136 )
5256	{
5257	r = -1.043e+00;
5258	}
5259	if( w==135 )
5260	{
5261	r = -1.073e+00;
5262	}
5263	if( w==134 )
5264	{
5265	r = -1.103e+00;
5266	}
5267	if( w==133 )
5268	{
5269	r = -1.135e+00;
5270	}
5271	if( w==132 )
5272	{
5273	r = -1.166e+00;
5274	}
5275	if( w==131 )
5276	{
5277	r = -1.198e+00;
5278	}
5279	if( w==130 )
5280	{
5281	r = -1.231e+00;
5282	}
5283	if( w==129 )
5284	{
5285	r = -1.265e+00;
5286	}
5287	if( w==128 )
5288	{
5289	r = -1.299e+00;
5290	}
5291	if( w==127 )
5292	{
5293	r = -1.334e+00;
5294	}
5295	if( w==126 )
5296	{
5297	r = -1.369e+00;
5298	}
5299	if( w==125 )
5300	{
5301	r = -1.405e+00;
5302	}
5303	if( w==124 )
5304	{
5305	r = -1.441e+00;
5306	}
5307	if( w==123 )
5308	{
5309	r = -1.479e+00;
5310	}
5311	if( w==122 )
5312	{
5313	r = -1.517e+00;
5314	}
5315	if( w==121 )
5316	{
5317	r = -1.555e+00;
5318	}
5319	if( w==120 )
5320	{
5321	r = -1.594e+00;
5322	}
5323	if( w==119 )
5324	{
5325	r = -1.634e+00;
5326	}
5327	if( w==118 )
5328	{
5329	r = -1.675e+00;
5330	}
5331	if( w==117 )
5332	{
5333	r = -1.716e+00;
5334	}
5335	if( w==116 )
5336	{
5337	r = -1.758e+00;
5338	}
5339	if( w==115 )
5340	{
5341	r = -1.800e+00;
5342	}
5343	if( w==114 )
5344	{
5345	r = -1.844e+00;
5346	}
5347	if( w==113 )
5348	{
5349	r = -1.888e+00;
5350	}
5351	if( w==112 )
5352	{
5353	r = -1.932e+00;
5354	}
5355	if( w==111 )
5356	{
5357	r = -1.978e+00;
5358	}
5359	if( w==110 )
5360	{
5361	r = -2.024e+00;
5362	}
5363	if( w==109 )
5364	{
5365	r = -2.070e+00;
5366	}
5367	if( w==108 )
5368	{
5369	r = -2.118e+00;
5370	}
5371	if( w==107 )
5372	{
5373	r = -2.166e+00;
5374	}
5375	if( w==106 )
5376	{
5377	r = -2.215e+00;
5378	}
5379	if( w==105 )
5380	{
5381	r = -2.265e+00;
5382	}
5383	if( w==104 )
5384	{
5385	r = -2.316e+00;
5386	}
5387	if( w==103 )
5388	{
5389	r = -2.367e+00;
5390	}
5391	if( w==102 )
5392	{
5393	r = -2.419e+00;
5394	}
5395	if( w==101 )
5396	{
5397	r = -2.472e+00;
5398	}
5399	if( w==100 )
5400	{
5401	r = -2.526e+00;
5402	}
5403	if( w==99 )
5404	{
5405	r = -2.580e+00;
5406	}
5407	if( w==98 )
5408	{
5409	r = -2.636e+00;
5410	}
5411	if( w==97 )
5412	{
5413	r = -2.692e+00;
5414	}
5415	if( w==96 )
5416	{
5417	r = -2.749e+00;
5418	}
5419	if( w==95 )
5420	{
5421	r = -2.806e+00;
5422	}
5423	if( w==94 )
5424	{
5425	r = -2.865e+00;
5426	}
5427	if( w==93 )
5428	{
5429	r = -2.925e+00;
5430	}
5431	if( w==92 )
5432	{
5433	r = -2.985e+00;
5434	}
5435	if( w==91 )
5436	{
5437	r = -3.046e+00;
5438	}
5439	if( w==90 )
5440	{
5441	r = -3.108e+00;
5442	}
5443	if( w==89 )
5444	{
5445	r = -3.171e+00;
5446	}
5447	if( w==88 )
5448	{
5449	r = -3.235e+00;
5450	}
5451	if( w==87 )
5452	{
5453	r = -3.300e+00;
5454	}
5455	if( w==86 )
5456	{
5457	r = -3.365e+00;
5458	}
5459	if( w==85 )
5460	{
5461	r = -3.432e+00;
5462	}
5463	if( w==84 )
5464	{
5465	r = -3.499e+00;
5466	}
5467	if( w==83 )
5468	{
5469	r = -3.568e+00;
5470	}
5471	if( w==82 )
5472	{
5473	r = -3.637e+00;
5474	}
5475	if( w==81 )
5476	{
5477	r = -3.708e+00;
5478	}
5479	if( w==80 )
5480	{
5481	r = -3.779e+00;
5482	}
5483	if( w==79 )
5484	{
5485	r = -3.852e+00;
5486	}
5487	if( w==78 )
5488	{
5489	r = -3.925e+00;
5490	}
5491	if( w==77 )
5492	{
5493	r = -4.000e+00;
5494	}
5495	if( w==76 )
5496	{
5497	r = -4.075e+00;
5498	}
5499	if( w==75 )
5500	{
5501	r = -4.151e+00;
5502	}
5503	if( w==74 )
5504	{
5505	r = -4.229e+00;
5506	}
5507	if( w==73 )
5508	{
5509	r = -4.308e+00;
5510	}
5511	if( w==72 )
5512	{
5513	r = -4.387e+00;
5514	}
5515	if( w==71 )
5516	{
5517	r = -4.468e+00;
5518	}
5519	if( w==70 )
5520	{
5521	r = -4.550e+00;
5522	}
5523	if( w==69 )
5524	{
5525	r = -4.633e+00;
5526	}
5527	if( w==68 )
5528	{
5529	r = -4.718e+00;
5530	}
5531	if( w==67 )
5532	{
5533	r = -4.803e+00;
5534	}
5535	if( w==66 )
5536	{
5537	r = -4.890e+00;
5538	}
5539	if( w==65 )
5540	{
5541	r = -4.978e+00;
5542	}
5543	if( w==64 )
5544	{
5545	r = -5.067e+00;
5546	}
5547	if( w==63 )
5548	{
5549	r = -5.157e+00;
5550	}
5551	if( w==62 )
5552	{
5553	r = -5.249e+00;
5554	}
5555	if( w==61 )
5556	{
5557	r = -5.342e+00;
5558	}
5559	if( w==60 )
5560	{
5561	r = -5.436e+00;
5562	}
5563	if( w==59 )
5564	{
5565	r = -5.531e+00;
5566	}
5567	if( w==58 )
5568	{
5569	r = -5.628e+00;
5570	}
5571	if( w==57 )
5572	{
5573	r = -5.727e+00;
5574	}
5575	if( w==56 )
5576	{
5577	r = -5.826e+00;
5578	}
5579	if( w==55 )
5580	{
5581	r = -5.927e+00;
5582	}
5583	if( w==54 )
5584	{
5585	r = -6.030e+00;
5586	}
5587	if( w==53 )
5588	{
5589	r = -6.134e+00;
5590	}
5591	if( w==52 )
5592	{
5593	r = -6.240e+00;
5594	}
5595	if( w==51 )
5596	{
5597	r = -6.347e+00;
5598	}
5599	if( w==50 )
5600	{
5601	r = -6.456e+00;
5602	}
5603	if( w==49 )
5604	{
5605	r = -6.566e+00;
5606	}
5607	if( w==48 )
5608	{
5609	r = -6.678e+00;
5610	}
5611	if( w==47 )
5612	{
5613	r = -6.792e+00;
5614	}
5615	if( w==46 )
5616	{
5617	r = -6.907e+00;
5618	}
5619	if( w==45 )
5620	{
5621	r = -7.025e+00;
5622	}
5623	if( w==44 )
5624	{
5625	r = -7.144e+00;
5626	}
5627	if( w==43 )
5628	{
5629	r = -7.265e+00;
5630	}
5631	if( w==42 )
5632	{
5633	r = -7.387e+00;
5634	}
5635	if( w==41 )
5636	{
5637	r = -7.512e+00;
5638	}
5639	if( w==40 )
5640	{
5641	r = -7.639e+00;
5642	}
5643	if( w==39 )
5644	{
5645	r = -7.768e+00;
5646	}
5647	if( w==38 )
5648	{
5649	r = -7.899e+00;
5650	}
5651	if( w==37 )
5652	{
5653	r = -8.032e+00;
5654	}
5655	if( w==36 )
5656	{
5657	r = -8.167e+00;
5658	}
5659	if( w==35 )
5660	{
5661	r = -8.305e+00;
5662	}
5663	if( w==34 )
5664	{
5665	r = -8.445e+00;
5666	}
5667	if( w==33 )
5668	{
5669	r = -8.588e+00;
5670	}
5671	if( w==32 )
5672	{
5673	r = -8.733e+00;
5674	}
5675	if( w==31 )
5676	{
5677	r = -8.881e+00;
5678	}
5679	if( w==30 )
5680	{
5681	r = -9.031e+00;
5682	}
5683	if( w==29 )
5684	{
5685	r = -9.185e+00;
5686	}
5687	if( w==28 )
5688	{
5689	r = -9.341e+00;
5690	}
5691	if( w==27 )
5692	{
5693	r = -9.501e+00;
5694	}
5695	if( w==26 )
5696	{
5697	r = -9.664e+00;
5698	}
5699	if( w==25 )
5700	{
5701	r = -9.830e+00;
5702	}
5703	if( w==24 )
5704	{
5705	r = -1.000e+01;
5706	}
5707	if( w==23 )
5708	{
5709	r = -1.017e+01;
5710	}
5711	if( w==22 )
5712	{
5713	r = -1.035e+01;
5714	}
5715	if( w==21 )
5716	{
5717	r = -1.053e+01;
5718	}
5719	if( w==20 )
5720	{
5721	r = -1.072e+01;
5722	}
5723	if( w==19 )
5724	{
5725	r = -1.091e+01;
5726	}
5727	if( w==18 )
5728	{
5729	r = -1.110e+01;
5730	}
5731	if( w==17 )
5732	{
5733	r = -1.130e+01;
5734	}
5735	if( w==16 )
5736	{
5737	r = -1.151e+01;
5738	}
5739	if( w==15 )
5740	{
5741	r = -1.172e+01;
5742	}
5743	if( w==14 )
5744	{
5745	r = -1.194e+01;
5746	}
5747	if( w==13 )
5748	{
5749	r = -1.216e+01;
5750	}
5751	if( w==12 )
5752	{
5753	r = -1.239e+01;
5754	}
5755	if( w==11 )
5756	{
5757	r = -1.263e+01;
5758	}
5759	if( w==10 )
5760	{
5761	r = -1.287e+01;
5762	}
5763	if( w==9 )
5764	{
5765	r = -1.314e+01;
5766	}
5767	if( w==8 )
5768	{
5769	r = -1.342e+01;
5770	}
5771	if( w==7 )
5772	{
5773	r = -1.369e+01;
5774	}
5775	if( w==6 )
5776	{
5777	r = -1.400e+01;
5778	}
5779	if( w==5 )
5780	{
5781	r = -1.433e+01;
5782	}
5783	if( w==4 )
5784	{
5785	r = -1.469e+01;
5786	}
5787	if( w==3 )
5788	{
5789	r = -1.503e+01;
5790	}
5791	if( w==2 )
5792	{
5793	r = -1.554e+01;
5794	}
5795	if( w==1 )
5796	{
5797	r = -1.594e+01;
5798	}
5799	if( w<=0 )
5800	{
5801	r = -1.664e+01;
5802	}
5803	result = r;
5804	return result;
5805	}
5806
5807
5808	/*************************************************************************
5809	Tail(S, 25)
5810	*************************************************************************/
5811	private static double w25(double s)
5812	{
5813	double result = 0;
5814	double x = 0;
5815	double tj = 0;
5816	double tj1 = 0;
5817
5818	result = 0;
5819	x = Math.Min(2*(s-0.000000e+00)/4.000000e+00-1, 1.0);
5820	tj = 1;
5821	tj1 = x;
5822	wcheb(x, -5.150509e+00, ref tj, ref tj1, ref result);
5823	wcheb(x, -5.695528e+00, ref tj, ref tj1, ref result);
5824	wcheb(x, -1.437637e+00, ref tj, ref tj1, ref result);
5825	wcheb(x, -2.611906e-01, ref tj, ref tj1, ref result);
5826	wcheb(x, -7.625722e-02, ref tj, ref tj1, ref result);
5827	wcheb(x, -2.579892e-02, ref tj, ref tj1, ref result);
5828	wcheb(x, -1.086876e-02, ref tj, ref tj1, ref result);
5829	wcheb(x, -2.906543e-03, ref tj, ref tj1, ref result);
5830	wcheb(x, -2.354881e-03, ref tj, ref tj1, ref result);
5831	wcheb(x, 1.007195e-04, ref tj, ref tj1, ref result);
5832	wcheb(x, -8.437327e-04, ref tj, ref tj1, ref result);
5833	return result;
5834	}
5835
5836
5837	/*************************************************************************
5838	Tail(S, 26)
5839	*************************************************************************/
5840	private static double w26(double s)
5841	{
5842	double result = 0;
5843	double x = 0;
5844	double tj = 0;
5845	double tj1 = 0;
5846
5847	result = 0;
5848	x = Math.Min(2*(s-0.000000e+00)/4.000000e+00-1, 1.0);
5849	tj = 1;
5850	tj1 = x;
5851	wcheb(x, -5.117622e+00, ref tj, ref tj1, ref result);
5852	wcheb(x, -5.635159e+00, ref tj, ref tj1, ref result);
5853	wcheb(x, -1.395167e+00, ref tj, ref tj1, ref result);
5854	wcheb(x, -2.382823e-01, ref tj, ref tj1, ref result);
5855	wcheb(x, -6.531987e-02, ref tj, ref tj1, ref result);
5856	wcheb(x, -2.060112e-02, ref tj, ref tj1, ref result);
5857	wcheb(x, -8.203697e-03, ref tj, ref tj1, ref result);
5858	wcheb(x, -1.516523e-03, ref tj, ref tj1, ref result);
5859	wcheb(x, -1.431364e-03, ref tj, ref tj1, ref result);
5860	wcheb(x, 6.384553e-04, ref tj, ref tj1, ref result);
5861	wcheb(x, -3.238369e-04, ref tj, ref tj1, ref result);
5862	return result;
5863	}
5864
5865
5866	/*************************************************************************
5867	Tail(S, 27)
5868	*************************************************************************/
5869	private static double w27(double s)
5870	{
5871	double result = 0;
5872	double x = 0;
5873	double tj = 0;
5874	double tj1 = 0;
5875
5876	result = 0;
5877	x = Math.Min(2*(s-0.000000e+00)/4.000000e+00-1, 1.0);
5878	tj = 1;
5879	tj1 = x;
5880	wcheb(x, -5.089731e+00, ref tj, ref tj1, ref result);
5881	wcheb(x, -5.584248e+00, ref tj, ref tj1, ref result);
5882	wcheb(x, -1.359966e+00, ref tj, ref tj1, ref result);
5883	wcheb(x, -2.203696e-01, ref tj, ref tj1, ref result);
5884	wcheb(x, -5.753344e-02, ref tj, ref tj1, ref result);
5885	wcheb(x, -1.761891e-02, ref tj, ref tj1, ref result);
5886	wcheb(x, -7.096897e-03, ref tj, ref tj1, ref result);
5887	wcheb(x, -1.419108e-03, ref tj, ref tj1, ref result);
5888	wcheb(x, -1.581214e-03, ref tj, ref tj1, ref result);
5889	wcheb(x, 3.033766e-04, ref tj, ref tj1, ref result);
5890	wcheb(x, -5.901441e-04, ref tj, ref tj1, ref result);
5891	return result;
5892	}
5893
5894
5895	/*************************************************************************
5896	Tail(S, 28)
5897	*************************************************************************/
5898	private static double w28(double s)
5899	{
5900	double result = 0;
5901	double x = 0;
5902	double tj = 0;
5903	double tj1 = 0;
5904
5905	result = 0;
5906	x = Math.Min(2*(s-0.000000e+00)/4.000000e+00-1, 1.0);
5907	tj = 1;
5908	tj1 = x;
5909	wcheb(x, -5.065046e+00, ref tj, ref tj1, ref result);
5910	wcheb(x, -5.539163e+00, ref tj, ref tj1, ref result);
5911	wcheb(x, -1.328939e+00, ref tj, ref tj1, ref result);
5912	wcheb(x, -2.046376e-01, ref tj, ref tj1, ref result);
5913	wcheb(x, -5.061515e-02, ref tj, ref tj1, ref result);
5914	wcheb(x, -1.469271e-02, ref tj, ref tj1, ref result);
5915	wcheb(x, -5.711578e-03, ref tj, ref tj1, ref result);
5916	wcheb(x, -8.389153e-04, ref tj, ref tj1, ref result);
5917	wcheb(x, -1.250575e-03, ref tj, ref tj1, ref result);
5918	wcheb(x, 4.047245e-04, ref tj, ref tj1, ref result);
5919	wcheb(x, -5.128555e-04, ref tj, ref tj1, ref result);
5920	return result;
5921	}
5922
5923
5924	/*************************************************************************
5925	Tail(S, 29)
5926	*************************************************************************/
5927	private static double w29(double s)
5928	{
5929	double result = 0;
5930	double x = 0;
5931	double tj = 0;
5932	double tj1 = 0;
5933
5934	result = 0;
5935	x = Math.Min(2*(s-0.000000e+00)/4.000000e+00-1, 1.0);
5936	tj = 1;
5937	tj1 = x;
5938	wcheb(x, -5.043413e+00, ref tj, ref tj1, ref result);
5939	wcheb(x, -5.499756e+00, ref tj, ref tj1, ref result);
5940	wcheb(x, -1.302137e+00, ref tj, ref tj1, ref result);
5941	wcheb(x, -1.915129e-01, ref tj, ref tj1, ref result);
5942	wcheb(x, -4.516329e-02, ref tj, ref tj1, ref result);
5943	wcheb(x, -1.260064e-02, ref tj, ref tj1, ref result);
5944	wcheb(x, -4.817269e-03, ref tj, ref tj1, ref result);
5945	wcheb(x, -5.478130e-04, ref tj, ref tj1, ref result);
5946	wcheb(x, -1.111668e-03, ref tj, ref tj1, ref result);
5947	wcheb(x, 4.093451e-04, ref tj, ref tj1, ref result);
5948	wcheb(x, -5.135860e-04, ref tj, ref tj1, ref result);
5949	return result;
5950	}
5951
5952
5953	/*************************************************************************
5954	Tail(S, 30)
5955	*************************************************************************/
5956	private static double w30(double s)
5957	{
5958	double result = 0;
5959	double x = 0;
5960	double tj = 0;
5961	double tj1 = 0;
5962
5963	result = 0;
5964	x = Math.Min(2*(s-0.000000e+00)/4.000000e+00-1, 1.0);
5965	tj = 1;
5966	tj1 = x;
5967	wcheb(x, -5.024071e+00, ref tj, ref tj1, ref result);
5968	wcheb(x, -5.464515e+00, ref tj, ref tj1, ref result);
5969	wcheb(x, -1.278342e+00, ref tj, ref tj1, ref result);
5970	wcheb(x, -1.800030e-01, ref tj, ref tj1, ref result);
5971	wcheb(x, -4.046294e-02, ref tj, ref tj1, ref result);
5972	wcheb(x, -1.076162e-02, ref tj, ref tj1, ref result);
5973	wcheb(x, -3.968677e-03, ref tj, ref tj1, ref result);
5974	wcheb(x, -1.911679e-04, ref tj, ref tj1, ref result);
5975	wcheb(x, -8.619185e-04, ref tj, ref tj1, ref result);
5976	wcheb(x, 5.125362e-04, ref tj, ref tj1, ref result);
5977	wcheb(x, -3.984370e-04, ref tj, ref tj1, ref result);
5978	return result;
5979	}
5980
5981
5982	/*************************************************************************
5983	Tail(S, 40)
5984	*************************************************************************/
5985	private static double w40(double s)
5986	{
5987	double result = 0;
5988	double x = 0;
5989	double tj = 0;
5990	double tj1 = 0;
5991
5992	result = 0;
5993	x = Math.Min(2*(s-0.000000e+00)/4.000000e+00-1, 1.0);
5994	tj = 1;
5995	tj1 = x;
5996	wcheb(x, -4.904809e+00, ref tj, ref tj1, ref result);
5997	wcheb(x, -5.248327e+00, ref tj, ref tj1, ref result);
5998	wcheb(x, -1.136698e+00, ref tj, ref tj1, ref result);
5999	wcheb(x, -1.170982e-01, ref tj, ref tj1, ref result);
6000	wcheb(x, -1.824427e-02, ref tj, ref tj1, ref result);
6001	wcheb(x, -3.888648e-03, ref tj, ref tj1, ref result);
6002	wcheb(x, -1.344929e-03, ref tj, ref tj1, ref result);
6003	wcheb(x, 2.790407e-04, ref tj, ref tj1, ref result);
6004	wcheb(x, -4.619858e-04, ref tj, ref tj1, ref result);
6005	wcheb(x, 3.359121e-04, ref tj, ref tj1, ref result);
6006	wcheb(x, -2.883026e-04, ref tj, ref tj1, ref result);
6007	return result;
6008	}
6009
6010
6011	/*************************************************************************
6012	Tail(S, 60)
6013	*************************************************************************/
6014	private static double w60(double s)
6015	{
6016	double result = 0;
6017	double x = 0;
6018	double tj = 0;
6019	double tj1 = 0;
6020
6021	result = 0;
6022	x = Math.Min(2*(s-0.000000e+00)/4.000000e+00-1, 1.0);
6023	tj = 1;
6024	tj1 = x;
6025	wcheb(x, -4.809656e+00, ref tj, ref tj1, ref result);
6026	wcheb(x, -5.077191e+00, ref tj, ref tj1, ref result);
6027	wcheb(x, -1.029402e+00, ref tj, ref tj1, ref result);
6028	wcheb(x, -7.507931e-02, ref tj, ref tj1, ref result);
6029	wcheb(x, -6.506226e-03, ref tj, ref tj1, ref result);
6030	wcheb(x, -1.391278e-03, ref tj, ref tj1, ref result);
6031	wcheb(x, -4.263635e-04, ref tj, ref tj1, ref result);
6032	wcheb(x, 2.302271e-04, ref tj, ref tj1, ref result);
6033	wcheb(x, -2.384348e-04, ref tj, ref tj1, ref result);
6034	wcheb(x, 1.865587e-04, ref tj, ref tj1, ref result);
6035	wcheb(x, -1.622355e-04, ref tj, ref tj1, ref result);
6036	return result;
6037	}
6038
6039
6040	/*************************************************************************
6041	Tail(S, 120)
6042	*************************************************************************/
6043	private static double w120(double s)
6044	{
6045	double result = 0;
6046	double x = 0;
6047	double tj = 0;
6048	double tj1 = 0;
6049
6050	result = 0;
6051	x = Math.Min(2*(s-0.000000e+00)/4.000000e+00-1, 1.0);
6052	tj = 1;
6053	tj1 = x;
6054	wcheb(x, -4.729426e+00, ref tj, ref tj1, ref result);
6055	wcheb(x, -4.934426e+00, ref tj, ref tj1, ref result);
6056	wcheb(x, -9.433231e-01, ref tj, ref tj1, ref result);
6057	wcheb(x, -4.492504e-02, ref tj, ref tj1, ref result);
6058	wcheb(x, 1.673948e-05, ref tj, ref tj1, ref result);
6059	wcheb(x, -6.077014e-04, ref tj, ref tj1, ref result);
6060	wcheb(x, -7.215768e-05, ref tj, ref tj1, ref result);
6061	wcheb(x, 9.086734e-05, ref tj, ref tj1, ref result);
6062	wcheb(x, -8.447980e-05, ref tj, ref tj1, ref result);
6063	wcheb(x, 6.705028e-05, ref tj, ref tj1, ref result);
6064	wcheb(x, -5.828507e-05, ref tj, ref tj1, ref result);
6065	return result;
6066	}
6067
6068
6069	/*************************************************************************
6070	Tail(S, 200)
6071	*************************************************************************/
6072	private static double w200(double s)
6073	{
6074	double result = 0;
6075	double x = 0;
6076	double tj = 0;
6077	double tj1 = 0;
6078
6079	result = 0;
6080	x = Math.Min(2*(s-0.000000e+00)/4.000000e+00-1, 1.0);
6081	tj = 1;
6082	tj1 = x;
6083	wcheb(x, -4.700240e+00, ref tj, ref tj1, ref result);
6084	wcheb(x, -4.883080e+00, ref tj, ref tj1, ref result);
6085	wcheb(x, -9.132168e-01, ref tj, ref tj1, ref result);
6086	wcheb(x, -3.512684e-02, ref tj, ref tj1, ref result);
6087	wcheb(x, 1.726342e-03, ref tj, ref tj1, ref result);
6088	wcheb(x, -5.189796e-04, ref tj, ref tj1, ref result);
6089	wcheb(x, -1.628659e-06, ref tj, ref tj1, ref result);
6090	wcheb(x, 4.261786e-05, ref tj, ref tj1, ref result);
6091	wcheb(x, -4.002498e-05, ref tj, ref tj1, ref result);
6092	wcheb(x, 3.146287e-05, ref tj, ref tj1, ref result);
6093	wcheb(x, -2.727576e-05, ref tj, ref tj1, ref result);
6094	return result;
6095	}
6096
6097
6098	/*************************************************************************
6099	Tail(S,N), S>=0
6100	*************************************************************************/
6101	private static double wsigma(double s,
6102	int n)
6103	{
6104	double result = 0;
6105	double f0 = 0;
6106	double f1 = 0;
6107	double f2 = 0;
6108	double f3 = 0;
6109	double f4 = 0;
6110	double x0 = 0;
6111	double x1 = 0;
6112	double x2 = 0;
6113	double x3 = 0;
6114	double x4 = 0;
6115	double x = 0;
6116
6117	if( n==5 )
6118	{
6119	result = w5(s);
6120	}
6121	if( n==6 )
6122	{
6123	result = w6(s);
6124	}
6125	if( n==7 )
6126	{
6127	result = w7(s);
6128	}
6129	if( n==8 )
6130	{
6131	result = w8(s);
6132	}
6133	if( n==9 )
6134	{
6135	result = w9(s);
6136	}
6137	if( n==10 )
6138	{
6139	result = w10(s);
6140	}
6141	if( n==11 )
6142	{
6143	result = w11(s);
6144	}
6145	if( n==12 )
6146	{
6147	result = w12(s);
6148	}
6149	if( n==13 )
6150	{
6151	result = w13(s);
6152	}
6153	if( n==14 )
6154	{
6155	result = w14(s);
6156	}
6157	if( n==15 )
6158	{
6159	result = w15(s);
6160	}
6161	if( n==16 )
6162	{
6163	result = w16(s);
6164	}
6165	if( n==17 )
6166	{
6167	result = w17(s);
6168	}
6169	if( n==18 )
6170	{
6171	result = w18(s);
6172	}
6173	if( n==19 )
6174	{
6175	result = w19(s);
6176	}
6177	if( n==20 )
6178	{
6179	result = w20(s);
6180	}
6181	if( n==21 )
6182	{
6183	result = w21(s);
6184	}
6185	if( n==22 )
6186	{
6187	result = w22(s);
6188	}
6189	if( n==23 )
6190	{
6191	result = w23(s);
6192	}
6193	if( n==24 )
6194	{
6195	result = w24(s);
6196	}
6197	if( n==25 )
6198	{
6199	result = w25(s);
6200	}
6201	if( n==26 )
6202	{
6203	result = w26(s);
6204	}
6205	if( n==27 )
6206	{
6207	result = w27(s);
6208	}
6209	if( n==28 )
6210	{
6211	result = w28(s);
6212	}
6213	if( n==29 )
6214	{
6215	result = w29(s);
6216	}
6217	if( n==30 )
6218	{
6219	result = w30(s);
6220	}
6221	if( n>30 )
6222	{
6223	x = 1.0/n;
6224	x0 = 1.0/30;
6225	f0 = w30(s);
6226	x1 = 1.0/40;
6227	f1 = w40(s);
6228	x2 = 1.0/60;
6229	f2 = w60(s);
6230	x3 = 1.0/120;
6231	f3 = w120(s);
6232	x4 = 1.0/200;
6233	f4 = w200(s);
6234	f1 = ((x-x0)f1-(x-x1)f0)/(x1-x0);
6235	f2 = ((x-x0)f2-(x-x2)f0)/(x2-x0);
6236	f3 = ((x-x0)f3-(x-x3)f0)/(x3-x0);
6237	f4 = ((x-x0)f4-(x-x4)f0)/(x4-x0);
6238	f2 = ((x-x1)f2-(x-x2)f1)/(x2-x1);
6239	f3 = ((x-x1)f3-(x-x3)f1)/(x3-x1);
6240	f4 = ((x-x1)f4-(x-x4)f1)/(x4-x1);
6241	f3 = ((x-x2)f3-(x-x3)f2)/(x3-x2);
6242	f4 = ((x-x2)f4-(x-x4)f2)/(x4-x2);
6243	f4 = ((x-x3)f4-(x-x4)f3)/(x4-x3);
6244	result = f4;
6245	}
6246	return result;
6247	}
6248	}
6249	}

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Update cookies preferences