1#---------------------------------------------------------------------------#
2# Math::NumberCruncher
3#       Date Written:   30-Aug-2000 02:41:52 PM
4#       Last Modified:  05-Mar-2002 12:29:30 PM
5#       Author:    Kurt Kincaid
6#       Copyright (c) 2002, Kurt Kincaid
7#           All Rights Reserved
8#
9# NOTICE:  Several of the algorithms contained herein are adapted from
10#          _Master Algorithms with Perl_, by John Orway, Jarkko Hietaniemi,
11#          and John Macdonald. Copyright (c) 1999 O'Reilly & Associates, Inc.
12#---------------------------------------------------------------------------#
13
14package Math::NumberCruncher;
15
16use Exporter;
17use constant epsilon => 1E-10;
18use Math::BigFloat;
19use strict;
20no strict 'refs';
21use vars qw( $PI $_e_ $_g_ $_ln2_ $max_ln2p $VERSION @ISA @EXPORT_OK @array $DECIMALS );
22
23@ISA       = qw( Exporter );
24@EXPORT_OK = qw( $PI $_e_ $_g_ $_ln2_ $VERSION log exp sqrt sin cos tan asin acos atan
25                 sec asec csc acsc exsec cot acot vers covers hav );
26
27$VERSION   = '5.00';
28
29$PI       = new Math::BigFloat "3.14159265358979323846264338327950288419716939937510582097494459230781640628620899862803482534211706798214808651328230664709384460955058223172535940812848111745028410270193852110555964462294895493038196442881097566593344612847564823378678316527120190914564856692346034861045432664821339360726024914127372458700660631558817488152092096282925409171536436789259036001133053054882046652138414695194151160943305727036575959195309218611738193261179310511854807446237996274956735188575272489122793818301194912983367336244065664308602139494639522473719070217986094370277053921717629317675238467481846766940513200056812714526356082778577134275778960917363717872146844090122495343014654958537105079227968925892354201995611212902196086403441815981362977477130996051870721134999999837297804995105973173281609631859502445945534690830264252230825334468503526193118817101000313783875288658753320838142061717766914730359825349042875546873115956286388235378759375195778185778053217122680661300192787661119590921642019893809525720106548586327886593615338182796823030195203530185296899577362259941389124972177528347913151557485724245415069595082953311686172785588907509838175463746493931925506040092770167113900984882401285836160356370766010471018194295559619894676783744944825537977472684710404753464620804668425906949129331367702898915210475216205696602405803815019351125338243003558764024749647326391419927260426992279678235478163600934172164121992458631503028618297455570674983850549458858692699569092721079750930295532116534498720275596023648066549911988183479775356636980742654252786255181841757467289097777279380008164706001614524919217321721477235014144197356854816136115735255213347574184946843852332390739414333454776241686251898356948556209921922218427255025425688767179049460165346680498862723279178608578438382796797668145410095388378636095068006422512520511739298489608412848862694560424196528502221066118630674427862203919494504712371378696095636437191728746776465757396241389086583264599581339047802759010";
30$_e_      = new Math::BigFloat "2.71828182845904523536028747135266249775724709369995957496696762772407663035354759457138217852516642742746639193200305992181741359662904357290033429526059563073813232862794349076323382988075319525101901157383418793070215408914993488416750924476146066808226480016847741185374234544243710753907774499206955170276183860626133138458300075204493382656029760673711320070932870912744374704723069697720931014169283681902551510865746377211125238978442505695369677078544996996794686445490598793163688923009879312773617821542499922957635148220826989519366803318252886939849646510582093923982948879332036250944311730123819706841614039701983767932068328237646480429531180232878250981945581530175671736133206981125099618188159304169035159888851934580727386673858942287922849989208680582574927961048419844436346324496848756023362482704197862320900216099023530436994184914631409343173814364054625315209618369088870701676839642437814059271456354906130310720851038375051011574770417189861068739696552126715468895703503540212340784981933432106817012100562788023519303322474501585390473041995777709350366041699732972508868769664035557071622684471625607988265178713419512466520103059212366771943252786753985589448969709640975459185695638023637016211204774272283648961342251644507818244235294863637214174023889344124796357437026375529444833799801612549227850925778256209262264832627793338656648162772516401910590049164499828931505660472580277863186415519565324425869829469593080191529872117255634754639644791014590409058629849679128740687050489585867174798546677575732056812884592054133405392200011378630094556068816674001698420558040336379537645203040243225661352783695117788386387443966253224985065499588623428189970773327617178392803494650143455889707194258639877275471096295374152111513683506275260232648472870392076431005958411661205452970302364725492966693811513732275364509888903136020572481765851180630364428123149655070475102544650117272115551948668508003685322818315219600373562527944951582841882947876108526398140";
31$_g_      = new Math::BigFloat "0.00000000006669531020394004460639036467721593281909711076035470516023410031617030523217887622766564072454302758797214777667825510044012806167171589236837418491695566945822976915357714542230787643797974413526391669997203504054665363724804035965562393645452314150103566079451722764077047780001159557565199157185300966536171598445697039986219614596562560614935883348444842355101781772391448082340919856836998916369518450095951586361599964956411488706712604230707700620231121148199618806694838144697445182374920538586757611082831191857644364217425269590709091396932217720806802248240057817186618417067158155776005255385640775647609250065179135759758507992396189592951980930481320790197399583643553004903638633848904759811394264488026638941042690422162661308367314767966094096732717350331446880659438788763653406174670608730191298841608778254576287486351270250356086586595348083207809542368578263528058113197006120223057460158115004058058346496610022508819904746818401803900084431805495574335423790092057583635142325381072166000392370432529697677943460796262298563755625126232411945538023315664461708646092717341756169185017815867224033437415993177148818991850540373256917580190547317016589079397744669794652571463427871642423999556656701735345089266437031981668333956202291758540353169755113915130162328378846078887754831021679812746264691773036658772927439184228333138447072436567005623282181435417359056639649685854700412554430269668778036840378088488720772714596464415222031008159188450589565739528427684396648817792962196229592100478823034546848482943080104006068776749469639024860584681362418381956826662905081375412550624813087051653933950489372422355734643741215840152422521918411076724517555944789360271997804891709581434658817443482396685154295502152134723282022085702770531064885514782676317483270433573003347808486077588831679942771722544158074564773981752974124512527118869659515788858519390280993535928473379654840896217869772206580519338744773792031249984043415900549966019072266610615037719";
32$max_ln2p = new Math::BigFloat "0.69314718055994530941723212145817656807550013436025525412068000949339362196969471560586332699641868754200148102057068573368552023575813055703267075163507596193072757082837143519030703862389167347112335011536449795523912047517268157493206515552473413952588295045300709532636664265410423915781495204374043038550080194417064167151864471283996817178454695702627163106454615025720740248163777338963855069526066834113727387372292895649354702576265209885969320196505855476470330679365443254763274495125040606943814710468994650622016772042452452961268794654619316517468139267250410380254625965686914419287160829380317271436778265487756648508567407764845146443994046142260319309673540257444607030809608504748663852313818167675143866747664789088143714198549423151997354880375165861275352916610007105355824987941472950929311389715599820565439287170007218085761025236889213244971389320378439353088774825970171559107088236836275898425891853530243634214367061189236789192372314672321720534016492568727477823445353476481149418642386776774406069562657379600867076257199184734022651462837904883062033061144630073719489002743643965002580936519443041191150608094879306786515887090060520346842973619384128965255653968602219412292420757432175748909770675268711581705113700915894266547859596489065305846025866838294002283300538207400567705304678700184162404418833232798386349001563121889560650553151272199398332030751408426091479001265168243443893572472788205486271552741877243002489794540196187233980860831664811490930667519339312890431641370681397776498176974868903887789991296503619270710889264105230924783917373501229842420499568935992206602204654941510613918788574424557751020683703086661948089641218680779020818158858000168811597305618667619918739520076671921459223672060253959543654165531129517598994005600036651356756905124592682574394648316833262490180382424082423145230614096380570070255138770268178516306902551370323405380214501901537402950994226299577964742713815736380172987394070424217997226696297993931270694";
33$_ln2_    = $max_ln2p->copy();
34
35my $max_piconst = $PI->copy();
36my $max_econst  = $_e_->copy();
37
38$DECIMALS = 20;
39
40sub new {
41    my $class = shift;
42    my $self = bless {}, $class;
43    return $self;
44}
45
46sub Range {
47    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
48    my $arrayref = shift;
49    return ( undef, undef ) unless defined $arrayref && @$arrayref > 0;
50    my ( $zzz, $hi, $lo );
51    $hi = $lo = $$arrayref[ 0 ];
52    foreach $zzz ( @$arrayref ) {
53        if ( $zzz > $hi ) {
54            $hi = $zzz;
55        }
56        if ( $zzz < $lo ) {
57            $lo = $zzz;
58        }
59    }
60    if ( $lo eq "" ) { $lo = "0" }
61    return ( $hi, $lo );
62}
63
64sub Mean {
65    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
66    my $arrayref = shift;
67    return undef unless defined $arrayref && @$arrayref > 0;
68    my $result;
69    foreach ( @$arrayref ) { $result += $_ }
70    return $result / @$arrayref;
71}
72
73sub Median {    # median may or may not be an element of the array
74    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
75    my $arrayref = shift;
76    my $P        = shift || $DECIMALS;
77    return undef unless defined $arrayref && @$arrayref > 0;
78    my $median = Math::BigFloat->new()->bfround( -$P );
79    my @array = sort { $a <=> $b } @$arrayref;
80    if ( @array % 2 ) {
81        $median = $array[ @array / 2 ];
82    } else {
83        $median = ( $array[ @array / 2 - 1 ] + $array[ @array / 2 ] ) / 2;
84    }
85    return $median;
86}
87
88sub OddMedian {    # median *is* an element of the array
89    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
90    my $arrayref = shift;
91    return undef unless defined $arrayref && @$arrayref > 0;
92    my @array = sort { $a <=> $b } @$arrayref;
93    return $array[ ( @array - ( 0, 0, 1, 0 )[ @array & 3 ] ) / 2 ];
94}
95
96sub Mode {
97    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
98    my $arrayref = shift;
99    return undef unless defined $arrayref && @$arrayref > 0;
100    my ( %count, @result );
101    foreach ( @$arrayref ) { $count{ $_ }++ }
102    foreach ( sort { $count{ $b } <=> $count{ $a } } keys %count ) {
103        last if @result && $count{ $_ } != $count{ $result[ 0 ] };
104        push ( @result, $_ );
105    }
106    return OddMedian \@result;
107}
108
109sub Covariance {
110    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
111    my $array1ref = shift;
112    my $array2ref = shift;
113    unless ( defined $array1ref && defined $array2ref && @$array1ref > 0 && $array2ref > 0 ) {
114        return undef;
115    }
116    my ( $i, $result );
117    for ( $i = 0 ; $i < @$array1ref ; $i++ ) {
118        $result += $array1ref->[ $i ] * $array2ref->[ $i ];
119    }
120    $result /= @$array1ref;
121    $result -= Mean( $array1ref ) * Mean( $array2ref );
122    return $result;
123}
124
125sub Correlation {
126    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
127    my $array1ref = shift;
128    my $array2ref = shift;
129    unless ( defined $array1ref && defined $array2ref && @$array1ref > 0 && $array2ref > 0 ) {
130        return undef;
131    }
132    my ( $sum1, $sum2, $sum1_sqrd, $sum2_sqrd );
133    foreach ( @$array1ref ) {
134        $sum1      += $_;
135        $sum1_sqrd += $_**2;
136    }
137    foreach ( @$array2ref ) {
138        $sum2      += $_;
139        $sum2_sqrd += $_**2;
140    }
141    return ( @$array1ref ** 2 ) * Covariance( $array1ref, $array2ref ) / SqrRoot(
142        abs( ( ( ( @$array1ref * $sum1_sqrd ) - ( $sum1 ** 2 ) ) * ( ( @$array1ref * $sum2_sqrd ) - ( $sum2 ** 2 ) ) ) ) );
143}
144
145sub BestFit {
146    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
147    my $a_ref = shift;
148    my $b_ref = shift;
149    my $P     = shift || $DECIMALS;
150    unless ( defined $a_ref && defined $b_ref && @$a_ref > 0 && @$b_ref > 0 ) {
151        return ( undef, undef );
152    }
153    my ( $i, $product, $sum1, $sum2, $sum1_sqrs, $a, $b );
154    $a         = Math::BigFloat->new();
155    $b         = $a->copy();
156    $sum1      = $a->copy();
157    $sum2      = $a->copy();
158    $sum1_sqrs = $a->copy();
159    for ( $i = 0 ; $i <= @$a_ref ; $i++ ) {
160        $product   += $a_ref->[ $i ] * $b_ref->[ $i ];
161        $sum1      += $a_ref->[ $i ];
162        $sum1_sqrs += $a_ref->[ $i ] ** 2;
163        $sum2      += $b_ref->[ $i ];
164    }
165    $b = ( ( @$a_ref * $product ) - ( $sum1 * $sum2 ) ) / ( ( @$a_ref * $sum1_sqrs ) - ( $sum1 ** 2 ) );
166    $a = ( $sum2 - $b * $sum1 ) / @$a_ref;
167    $a->bfround( -$P );
168    $b->bfround( -$P );
169    return ( $b, $a );
170}
171
172sub Distance {    # Distance( $x1, $y1, $x2, $y2 );
173    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
174    my @p = @_;
175    my $P;
176    if ( @p % 2 ) {
177        $P = pop @p;
178    } else {
179        $P = $DECIMALS;
180    }
181    return undef unless @p >= 3;
182    my $d = @p / 2;
183    return SqrRoot( abs( ( $_[ 0 ] - $_[ 2 ] ) ** 2 + ( $_[ 1 ] - $_[ 3 ] ) ** 2 ), $P ) if $d == 2;
184    my $S  = 0;
185    my @p0 = splice @p, 0, $d;
186    for ( my $i = 0 ; $i < $d ; $i++ ) {
187        my $di = $p0[ $i ] - $p[ $i ];
188        $S += $di * $di;
189    }
190    return SqrRoot( abs( $S ), $P );
191}
192
193sub ManhattanDistance {
194    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
195    my @p = @_;
196    return undef unless @p >= 3;
197    my $d  = @p / 2;
198    my $S  = 0;
199    my @p0 = splice @p, 0, $d;
200    for ( my $i = 0 ; $i < $d ; $i++ ) {
201        my $di = $p0[ $i ] - $p[ $i ];
202        $S += abs $di;
203    }
204    return $S;
205}
206
207sub AllOf {
208    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
209    my $result = 1;
210    my @array  = @_;
211    return undef unless @array >= 2;
212    while ( @array ) {
213        $result *= shift @array;
214    }
215    return $result;
216}
217
218sub NoneOf {
219    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
220    my $result = 1;
221    @array = @_;
222    foreach my $item ( @array ) {
223        $result *= ( 1 - $item );
224    }
225    return $result;
226}
227
228sub SomeOf {
229    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
230    @array = @_;
231    return undef unless @array >= 2;
232    return 1 - NoneOf( @array );
233}
234
235sub Factorial {
236    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
237    my $n = shift;
238    return undef unless defined $n;
239    my $result = Math::BigFloat->new( 1 );
240    unless ( $n >= 0 && $n == int( $n ) ) {
241        return undef;
242    }
243    while ( $n > 1 ) {
244        $result *= $n--;
245    }
246    return $result;
247}
248
249
250sub Permutation {
251    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
252    my ( $n, $k ) = @_;
253    return undef unless defined $n;
254    my $result = Math::BigFloat->new( 1 );
255    defined $k or $k = $n;
256    while ( $k-- ) { $result *= $n-- }
257    return $result;
258}
259
260sub Dice {
261    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
262    my $number = shift || 1;
263    my $sides  = shift || 6;
264    my $plus   = shift;
265    while ( $number-- ) {
266        $plus += int( rand( $sides ) + 1 );
267    }
268    return $plus;
269}
270
271sub RandInt {
272    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
273    my $low  = shift || 0;
274    my $high = shift || 1;
275    if ( $low > $high ) {
276        ( $low, $high ) = ( $high, $low );
277    }
278    return $low + int( rand( $high - $low + 1 ) );
279}
280
281sub RandomElement {
282    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
283    my $arrayref = shift;
284    $arrayref->[ rand @{ $arrayref } ];
285}
286
287sub ShuffleArray {
288    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
289    my $arrayref = shift;
290    return undef unless defined $arrayref && @$arrayref > 0;
291    for ( my $i = @$arrayref ; --$i ; ) {
292        my $j = int rand( $i + 1 );
293        next if $i == $j;
294        @$arrayref[ $i, $j ] = @$arrayref[ $j, $i ];
295    }
296}
297
298sub Unique {
299    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
300    my $arrayref = shift;
301    my %seen;
302    my $zzz;
303    my @unique;
304    return undef unless defined $arrayref && @$arrayref > 0;
305    foreach $zzz ( @$arrayref ) {
306        push ( @unique, $zzz ) unless $seen{ $zzz }++;
307    }
308    return @unique;
309}
310
311sub Compare {
312    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
313    my ( $arrayref1, $arrayref2 ) = @_;
314    unless ( defined $arrayref1 && defined $arrayref2 && @$arrayref1 > 0 && @$arrayref2 > 0 ) {
315        return undef;
316    }
317    my %seen;
318    my @aonly;
319    my $item;
320    foreach $item ( @$arrayref2 ) { $seen{ $item } = 1 }
321    foreach $item ( @$arrayref1 ) {
322        unless ( $seen{ $item } ) {
323            push ( @aonly, $item );
324        }
325    }
326    return @aonly;
327}
328
329sub Union {
330    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
331    my ( $arrayref1, $arrayref2 ) = @_;
332    unless ( defined $arrayref1 && defined $arrayref2 && @$arrayref1 > 0 && @$arrayref2 > 0 ) {
333        return undef;
334    }
335    my ( @union, @temp );
336    my %union;
337    my $zzz;
338    foreach $zzz ( @$arrayref1 ) { $union{ $zzz } = 1 }
339    foreach $zzz ( @$arrayref2 ) { $union{ $zzz } = 1 }
340    return keys %union;
341}
342
343sub Intersection {
344    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
345    my ( $arrayref1, $arrayref2 ) = @_;
346    unless ( defined $arrayref1 && defined $arrayref2 && @$arrayref1 > 0 && @$arrayref2 > 0 ) {
347        return undef;
348    }
349    my @isect = undef;
350    my ( %isect, %union, %count );
351    my $zzz;
352    foreach $zzz ( @$arrayref1 ) {
353        $union{ $zzz } = 1;
354    }
355    foreach $zzz ( @$arrayref2 ) {
356        if ( $union{ $zzz } ) {
357            $isect{ $zzz } = 1;
358        }
359    }
360    @isect = keys %isect;
361    return @isect;
362}
363
364sub Difference {
365    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
366    my ( $arrayref1, $arrayref2 ) = @_;
367    unless ( defined $arrayref1 && defined $arrayref2 && @$arrayref1 > 0 && @$arrayref2 > 0 ) {
368        return undef;
369    }
370    my ( @isect, @diff, @union ) = undef;
371    my $zzz;
372    my %count;
373    foreach $zzz ( @$arrayref1, @$arrayref2 ) { $count{ $zzz }++ }
374    foreach $zzz ( keys %count ) {
375        push @union, $zzz;
376        push @{ $count{ $zzz } > 1 ? \@isect : \@diff }, $zzz;
377    }
378    return @diff;
379}
380
381sub GaussianRand {
382    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
383    my ( $u1, $u2, $w, $g1, $g2 );
384    do {
385        $u1 = 2 * rand() - 1;
386        $u2 = 2 * rand() - 1;
387        $w  = $u1 * $u1 + $u2 * $u2;
388    } while ( $w >= 1 );
389    $w  = sqrt( abs( ( -2 * log( $w ) ) / $w ) );
390    $g2 = $u1 * $w;
391    $g1 = $u2 * $w;
392    return wantarray ? ( $g1, $g2 ) : $g1;
393}
394
395sub Choose {    # Probability of getting $k heads is $n tosses
396    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
397    my ( $n, $k ) = @_;
398    return undef unless defined $n && defined $k;
399    my ( $result, $j ) = ( 1, 1 );
400    if ( $k > $n || $k < 0 ) {
401        return 0;
402    }
403    while ( $j <= $k ) {
404        $result *= $n--;
405        $result /= $j++;
406    }
407    return $result;
408}
409
410sub Binomial {    # probability of $k successes in $n attempts, given probability of $p
411    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
412    my ( $n, $k, $p ) = @_;
413    return $k == 0 if $p == 0;
414    return $k != $n if $p == 1;
415    return Choose( $n, $k ) * $p ** $k * ( 1 - $p ) ** ( $n - $k );
416}
417
418sub GaussianDist {
419    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
420    use constant two_pi_sqrt_inverse => 1 / sqrt( 8 * atan2( 1, 1 ) );
421    my ( $x, $mean, $variance ) = @_;
422    return two_pi_sqrt_inverse * exp( -( $x - $mean ) ** 2 / ( 2 * $variance ) ) / SqrRoot( abs( $variance ) );
423}
424
425sub StandardDeviation {
426    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
427    my $arrayref = shift;
428    my $P        = shift || $DECIMALS;
429    return undef unless defined $arrayref && @$arrayref > 0;
430    my $mean = Mean( $arrayref );
431    return SqrRoot( abs( Mean( [ map $_**2, @$arrayref ] ) - ( $mean**2 ) ), $P );
432}
433
434sub Variance {
435    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
436    my $arrayref = shift;
437    my $P        = shift || $DECIMALS;
438    my $modP     = modP( $P );
439    return undef unless defined $arrayref && @$arrayref > 0;
440    my $result = StandardDeviation( $arrayref, $modP ) ** 2;
441    return $result->bfround( -$P );
442}
443
444sub StandardScores {    # number of StdDevs above the mean for each element
445    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
446    my $arrayref = shift;
447    my $P        = shift || $DECIMALS;
448    return undef unless defined $arrayref && @$arrayref > 0;
449    my $mean = Mean( $arrayref );
450    my ( $i, @scores );
451    my $deviation = StandardDeviation( $arrayref, $P );
452    return unless $deviation;
453
454    for ( $i = 0 ; $i < @$arrayref ; $i++ ) {
455        push @scores, ( $arrayref->[ $i ] - $mean ) / $deviation;
456    }
457    return @scores;
458}
459
460sub SignSignificance {
461    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
462    my ( $trials, $hits, $probability ) = @_;
463    return undef unless defined $trials && defined $hits && defined $probability;
464    my $confidence;
465    foreach ( $hits .. $trials ) {
466        $confidence += Binomial( $trials, $hits, $probability );
467    }
468    return $confidence;
469}
470
471sub EMC2 {
472    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
473    my $var  = shift;
474    my $unit = shift;
475    my $P    = shift || $DECIMALS;
476    return undef unless defined $var && defined $unit;
477    my $modP = modP( $P );
478    my $C;
479    if ( $unit =~ /^k/i ) {
480        $C = 299792.458;    # km per second
481    } elsif ( $unit =~ /^m/i ) {
482        $C = 186282.056;    # miles per second
483    } else {
484        return undef;
485    }
486    my $result = Math::BigFloat->new();
487    my $sqrd   = Math::BigFloat->new( $C );
488    $sqrd->bpow( 2 );
489    if ( $var =~ /^m(.*)$/i ) {
490        my $val = $1;
491        $result = $sqrd->copy()->bmul( $val, $modP );
492    } elsif ( $var =~ /^e(.*)$/i ) {
493        my $val = Math::BigFloat->new( $1 );
494        $result = $val->fdiv( $sqrd, $modP );
495    } else {
496        return undef;
497    }
498    return $result->bfround( -$P );
499}
500
501sub FMA {
502    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
503    my @vars = @_;
504    my $P;
505    if ( $vars[ -1 ] =~ /^\d*$/ ) {
506        $P = pop @vars;
507    } else {
508        $P = $DECIMALS;
509    }
510    my $modP = modP( $P );
511    @vars = sort @vars;
512    my ( $acc, $force, $mass );
513    my $result;
514    if ( $vars[ 0 ] =~ /^[Aa](.*)$/ ) {
515        $acc = $1;
516    } elsif ( $vars[ 0 ] =~ /^[Ff](.*)$/ ) {
517        $force = $1;
518    }
519    if ( $vars[ 1 ] =~ /^[Ff](.*)$/ ) {
520        $force = $1;
521    } elsif ( $vars[ 1 ] =~ /^[Mm](.*)$/ ) {
522        $mass = $1;
523    }
524    if ( $acc && $force ) {
525        $result = Math::BigFloat->new( $force );
526        $result->bdiv( $acc, $modP );
527    } elsif ( $acc && $mass ) {
528        $result = Math::BigFloat->new( $acc );
529        $result->bmul( $mass, $modP );
530    } elsif ( $force && $mass ) {
531        $result = Math::BigFloat->new( $force );
532        $result->bmul( $mass, $modP );
533    } else {
534        return undef;
535    }
536    return $result->bfround( -$P );
537}
538
539sub Predict {
540    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
541    my $slope       = shift;
542    my $y_intercept = shift;
543    my $proposed    = shift;
544    my $P           = shift || $DECIMALS;
545    my $modP        = modP( $P );
546    my $result = Math::BigFloat->new( $slope );
547    $result->bmul( $proposed, $modP )->badd( $y_intercept );
548    $result->bfround( -$P );
549    return $result;
550}
551
552sub TriangleHeron {
553    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
554    my ( $a, $b, $c, $P );
555    if ( @_ == 4 || @_ == 7 ) {
556        $P = pop @_;
557    } else {
558        $P = $DECIMALS;
559    }
560    my $modP = modP( $P );
561    if ( @_ == 3 ) {
562        ( $a, $b, $c ) = @_;
563    } elsif ( @_ == 6 ) {
564        ( $a, $b, $c ) = (
565          Distance( $_[ 0 ], $_[ 1 ], $_[ 2 ], $_[ 3 ], $modP ), Distance( $_[ 2 ], $_[ 3 ], $_[ 4 ], $_[ 5 ], $modP ),
566          Distance( $_[ 4 ], $_[ 5 ], $_[ 0 ], $_[ 1 ], $modP )
567        );
568    } else {
569        return undef;
570    }
571    my $s    = Math::BigFloat->new();
572    $s       = ( $a + $b + $c ) / 2;
573    my $root = $s * ( $s - $a ) * ( $s - $b ) * ( $s - $c );
574    return SqrRoot( abs( $root ), $P );
575}
576
577sub PolygonPerimeter {
578    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
579    my @xy = @_;
580    my $P;
581    if ( $xy[ -1 ] =~ /p/i ) {
582        ( $P = pop @xy ) =~ s/\D//g;
583    } else {
584        $P = $DECIMALS;
585    }
586    my $modP = modP( $P );
587    my $PP   = Math::BigFloat->new();
588    return undef unless @xy % 2 == 0 && @xy > 0;
589    for ( my ( $xa, $ya ) = @xy[ -2, -1 ] ; my ( $xb, $yb ) = splice @xy, 0, 2 ; ( $xa, $ya ) = ( $xb, $yb ) ) {
590        $PP += Distance( $xa, $ya, $xb, $yb, $modP );
591    }
592    $PP->bfround( -$P );
593    return $PP;
594}
595
596sub Clockwise {
597    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
598    my ( $x0, $y0, $x1, $y1, $x2, $y2 ) = @_;
599    return undef unless defined $x0 && defined $y0 && defined $x1 && defined $y1 && defined $x2 && defined $y2;
600    return ( $x2 - $x0 ) * ( $y1 - $y0 ) - ( $x1 - $x0 ) * ( $y2 - $y0 );
601}
602
603sub InPolygon {
604    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
605    my ( $x, $y, @xy ) = @_;
606    return undef unless defined $x && defined $y && @xy > 0;
607    my $n = @xy / 2;
608    my @i = map { 2 * $_ } 0 .. ( @xy / 2 );
609    my @x = map { $xy[ $_ ] } @i;
610    my @y = map { $xy[ $_ + 1 ] } @i;
611    my ( $i, $j );
612    my $side = 0;
613
614    for ( $i = 0, $j = $n - 1 ; $i < $n ; $j = $i++ ) {
615        if ( ( ( ( $y[ $i ] <= $y ) && ( $y < $y[ $j ] ) ) || ( ( $y[ $j ] <= $y ) && ( $y < $y[ $i ] ) ) )
616            and ( $x < ( $x[ $j ] - $x[ $i ] ) * ( $y - $y[ $i ] ) / ( $y[ $j ] - $y[ $i ] ) + $x[ $i ] ) )
617        {
618            $side = not $side;
619        }
620    }
621    return $side ? 1 : 0;
622}
623
624sub BoundingBox_Points {
625    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
626    my ( $d, @points ) = @_;
627    return undef unless defined $d && @points > 0;
628    my @bb;
629    while ( my @p = splice @points, 0, $d ) {
630        @bb = BoundingBox( $d, @p, @bb );
631    }
632    return @bb;
633}
634
635sub BoundingBox {
636    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
637    my ( $d, @bb ) = @_;
638    return undef unless defined $d && @bb > 0;
639    my @p = splice( @bb, 0, @bb - 2 * $d );
640    @bb = ( @p, @p ) unless @bb;
641    for ( my $i = 0 ; $i < $d ; $i++ ) {
642        for ( my $j = 0 ; $j < @p ; $j += $d ) {
643            my $ij = $i + $j;
644            $bb[ $i ] = $p[ $ij ] if $p[ $ij ] < $bb[ $i ];
645            $bb[ $i + $d ] = $p[ $ij ] if $p[ $ij ] > $bb[ $i + $d ];
646        }
647    }
648    return @bb;
649}
650
651sub InTriangle {
652    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
653    my ( $x, $y, $x0, $y0, $x1, $y1, $x2, $y2 ) = @_;
654    return undef unless defined defined $x
655      && defined $y
656      && defined $x0
657      && defined $y0
658      && defined $x1
659      && defined $y1
660      && defined $x2
661      && defined $y2;
662    my $cw0 = Clockwise( $x0, $y0, $x1, $y1, $x, $y );
663    return 1 if abs( $cw0 ) < epsilon;
664    my $cw1 = Clockwise( $x1, $y1, $x2, $y2, $x, $y );
665    return 1 if abs( $cw1 ) < epsilon;
666    return 0 if ( $cw0 < 0 and $cw1 > 0 ) or ( $cw0 > 0 and $cw1 < 0 );
667    my $cw2 = Clockwise( $x2, $y2, $x0, $y0, $x, $y );
668    return 1 if abs( $cw2 ) < epsilon;
669    return 0 if ( $cw0 < 0 and $cw2 > 0 ) or ( $cw0 > 0 and $cw2 < 0 );
670    return 1;
671}
672
673sub PolygonArea {
674    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
675    my @xy = @_;
676    my $P;
677    if ( $xy[ -1 ] =~ /p/i ) {
678        ( $P = pop @xy ) =~ s/\D//g;
679    } else {
680        $P = $DECIMALS;
681    }
682    my $modP = modP( $P );
683    return undef unless @xy % 2 == 0 && @xy > 0;
684    my $A = Math::BigFloat->new();
685    for ( my ( $xa, $ya ) = @xy[ -2, -1 ] ; my ( $xb, $yb ) = splice @xy, 0, 2 ; ( $xa, $ya ) = ( $xb, $yb ) ) {
686        $A += Determinant( $xa, $ya, $xb, $yb, $modP );
687    }
688    $A->bdiv( 2, $modP );
689    $A->babs()->bfround( -$P );
690    return $A;
691}
692
693sub Determinant {
694    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
695    my $result = Math::BigFloat->new();
696    $result = $_[ 0 ] * $_[ 3 ] - $_[ 1 ] * $_[ 2 ];
697    return $result;
698}
699
700sub CircleArea {
701    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
702    my $radius = shift;
703    my $P      = shift || $DECIMALS;
704    my $modP   = modP( $P );
705    return undef unless defined $radius;
706    my $area = Math::BigFloat->new( 1 );
707    $area = $PI->copy()->bmul( ( $radius ** 2 ), $modP );
708    $area->bfround( -$P );
709    return $area;
710}
711
712sub Circumference {
713    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
714    my $diameter = shift;
715    my $P        = shift || $DECIMALS;
716    return undef unless defined $diameter;
717    my $modP = modP( $P );
718    my $circumference = $PI->copy()->bmul( $diameter, $modP );
719    $circumference->bfround( -$P );
720    return $circumference;
721}
722
723sub SphereVolume {
724    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
725    my $radius = shift;
726    my $P      = shift || $DECIMALS;
727    return undef unless defined $radius;
728    my $modP = modP( $P );
729    my $volume = $PI->copy()->bmul( ( 4 / 3 ), $modP )->bmul( ( $radius ** 3 ), $modP );
730    $volume->bfround( -$P );
731    return $volume;
732}
733
734sub SphereSurface {
735    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
736    my $radius = shift;
737    my $P      = shift || $DECIMALS;
738    return undef unless defined $radius;
739    my $modP = modP( $P );
740    my $surface = $PI->copy()->bmul( 4, $modP )->bmul( ( $radius ** 2 ), $modP );
741    $surface->bfround( -$P );
742    return $surface;
743}
744
745sub RuleOf72 {
746    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
747    my $pct  = shift;
748    my $P    = shift || $DECIMALS;
749    my $modP = modP( $P );
750    return undef unless defined $pct;
751    my $num = Math::BigFloat->new( 72 );
752    $num->bdiv( $pct, $modP )->bfround( -$P );
753    return $num;
754}
755
756sub CylinderVolume {
757    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
758    my $radius = shift;
759    my $height = shift;
760    my $P      = shift || $DECIMALS;
761    return undef unless defined $radius && defined $height;
762    my $modP = modP( $P );
763    my $volume = $PI->copy()->bmul( ( $radius ** 2 ), $modP )->bmul( $height, $modP );
764    $volume->bfround( -$P );
765    return $volume;
766}
767
768sub ConeVolume {
769    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
770    my $lowerbase = shift;
771    my $height    = shift;
772    my $P         = shift || $DECIMALS;
773    return undef unless defined $lowerbase && defined $height;
774    my $modP = modP( $P );
775    my $num = Math::BigFloat->new( $lowerbase );
776    $num->bmul( $height, $modP )->bdiv( 3, $modP );
777    $num->bfround( -$P );
778    return $num;
779}
780
781sub deg2rad {
782    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
783    my $degrees = shift;
784    my $P       = shift || $DECIMALS;
785    return undef unless defined $degrees;
786    my $modP = modP( $P );
787    my $radians = Math::BigFloat->new( $degrees );
788    $radians->bdiv( 180, $modP )->bmul( $PI->copy(), $modP );
789    $radians->bfround( -$P );
790    return $radians;
791}
792
793sub rad2deg {
794    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
795    my $radians = shift;
796    my $P       = shift || $DECIMALS;
797    return undef unless defined $radians;
798    my $modP = modP( $P );
799    my $degrees = Math::BigFloat->new( $radians );
800    $degrees->bdiv( $PI->copy(), $modP )->bmul( 180, $modP );
801    $degrees->bfround( -$P );
802    return $degrees;
803}
804
805sub C2F {
806    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
807    my $x = Math::BigFloat->new( shift );
808    my $P = shift || $DECIMALS;
809    return undef unless defined $x;
810    my $degrees = Math::BigFloat->new( $x );
811    my $modP = modP( $P );
812    $degrees->bmul( 1.8, $modP )->badd( 32 )->bfround( -$P );
813    return $degrees;
814}
815
816sub F2C {
817    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
818    my $x = Math::BigFloat->new( shift );
819    my $P = shift || $DECIMALS;
820    return undef unless defined $x;
821    my $degrees = Math::BigFloat->new( $x );
822    my $modP = modP( $P );
823    $degrees->bsub( 32 )->bdiv( 1.8, $modP )->bfround( -$P );
824    return $degrees;
825}
826
827sub cm2in {
828    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
829    my $x = Math::BigFloat->new( shift );
830    my $P  = shift || $DECIMALS;
831    return undef unless defined $x;
832    my $cm   = Math::BigFloat->new( $x );
833    my $modP = modP( $P );
834    $cm->bmul( 0.3937007874, $modP )->bfround( -$P );
835    return $cm;
836}
837
838sub in2cm {
839    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
840    my $x = Math::BigFloat->new( shift );
841    my $P = shift || $DECIMALS;
842    return undef unless defined $x;
843    my $inches = Math::BigFloat->new( $x );
844    my $modP   = modP ( $P );
845    $inches->bmul( 2.54, $modP )->bfround( -$P );
846    return $inches;
847}
848
849sub m2ft {
850    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
851    my $temp   = shift;
852    my $P      = shift || $DECIMALS;
853    return undef unless defined $temp;
854    my $modP = modP( $P );
855    my $meters = Math::BigFloat->new( $temp );
856    $meters->bmul( 3.280839895, $modP )->bfround( -$P );
857    return $meters;
858}
859
860sub ft2m {
861    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
862    my $temp = shift;
863    my $P    = shift || $DECIMALS;
864    return undef unless defined $temp;
865    my $modP = modP( $P );
866    my $feet = Math::BigFloat->new( $temp );
867    $feet->bmul( 0.3048, $modP )->bfround( -$P );
868    return $feet;
869}
870
871sub km2miles {
872    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
873    my $x = shift;
874    my $P = shift || $DECIMALS;
875    return undef unless defined $x;
876    my $modP = modP( $P );
877    my $miles = Math::BigFloat->new( $x );
878    $miles->bmul( 0.6213711922, $modP )->bfround( -$P );
879    return $miles;
880}
881
882sub miles2km {
883    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
884    my $x = shift;
885    my $P = shift || $DECIMALS;
886    return undef unless defined $x;
887    my $modP = modP( $P );
888    my $km = Math::BigFloat->new( $x );
889    $km->bmul( 1.6093440000966945374266097172311, $modP )->bfround( -$P );
890    return $km;
891}
892
893sub kg2lb {
894    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
895    my $temp = shift;
896    my $P    = shift || $DECIMALS;
897    return undef unless defined $temp;
898    my $modP = modP( $P );
899    my $kg   = Math::BigFloat->new( $temp );
900    $kg->bmul( 2.204622622, $modP )->bfround( -$P );
901    return $kg;
902}
903
904sub lb2kg {
905    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
906    my $temp = shift;
907    my $P    = shift || $DECIMALS;
908    return undef unless defined $temp;
909    my $modP = modP( $P );
910    my $lb   = Math::BigFloat->new( $temp );
911    $lb->bmul( 0.45359237, $modP )->bfround( -$P );
912    return $lb;
913}
914
915sub RelativeStride {
916    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
917    my $stride_length = shift;
918    my $leg_length    = shift;
919    my $P             = shift || $DECIMALS;
920    return undef unless defined $stride_length && defined $leg_length;
921    my $modP = modP( $P );
922    my $rs = Math::BigFloat->new( $stride_length );
923    $rs->bdiv( $leg_length, $modP )->bfround( -$P );
924    return $rs;
925}
926
927sub RelativeStride_2 {
928    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
929    my $temp = shift;
930    my $P    = shift || $DECIMALS;
931    return undef unless defined $temp;
932    my $modP = modP( $P );
933    my $DS   = Math::BigFloat->new( $temp );
934    $DS->bmul( 1.1, $modP )->badd( 1 )->bfround( -$P );
935    return $DS;
936}
937
938sub DimensionlessSpeed {
939    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
940    my $temp = shift;
941    my $P    = shift || $DECIMALS;
942    return undef unless defined $temp;
943    my $modP = modP( $P );
944    my $RSL  = Math::BigFloat->new( $temp );
945    $RSL->bsub( 1 )->bdiv( 1.1, $modP )->bfround( -$P );
946    return $RSL;
947}
948
949sub DimensionlessSpeed_2 {
950    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
951    my $speed     = shift;
952    my $legLength = shift;
953    my $P         = shift || $DECIMALS;
954    return undef unless defined $speed && defined $legLength;
955    my $modP = modP( $P );
956    my $DS   = Math::BigFloat->new( $speed );
957    my $root = Math::BigFloat->new( $legLength );
958    $root->bmul( 9.80665, $modP )->babs();
959    my $sqrt = SqrRoot( $root->bstr(), $modP );
960    $DS->bdiv( $sqrt->bstr() )->bfround( -$P );
961    return $DS;
962}
963
964sub ActualSpeed {
965    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
966    my $legLength = Math::BigFloat->new( shift );
967    my $AS        = Math::BigFloat->new( shift );
968    my $P         = shift || $DECIMALS;
969    my $modP      = modP( $P );
970    $legLength->bmul( 9.80665, $modP )->babs();
971    my $root = Root2( $legLength->bstr(), 2, $modP );
972    $AS->bmul( $root->bstr(), $modP );
973    return $AS->bfround( -$P );
974}
975
976sub Eccentricity {
977    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
978    my $a    = Math::BigFloat->new( shift );
979    my $b    = Math::BigFloat->new( shift );
980    my $P    = shift || $DECIMALS;
981    my $modP = modP( $P );
982    my $root = Math::BigFloat->new();
983    my $A    = $a->copy();
984    $a->bpow( 2 );
985    $b->bpow( 2 );
986    $a->bsub( $b->copy() )->babs();
987    $root = SqrRoot( $a->bstr(), $modP );
988    $root /= $A;
989    return $root->bfround( -$P );
990}
991
992sub LatusRectum {
993    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
994    my $a    = shift;
995    my $b    = shift;
996    my $P    = shift || $DECIMALS;
997    unless ( $a && $b ) { return undef }
998    my $modP = modP( $P );
999    my $result = Math::BigFloat->new( 2 );
1000    $result->bmul( $b ** 2, $modP )->bdiv( $a, $modP )->bfround( -$P );
1001    return $result;
1002}
1003
1004sub EllipseArea {
1005    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1006    my $a = shift;
1007    my $b = shift;
1008    my $P = shift || $DECIMALS;
1009    unless ( $a && $b ) { return undef }
1010    my $modP = modP( $P );
1011    my $area = Math::BigFloat->new( $a );
1012    $area->bmul( $b, $modP )->bmul( $PI->copy(), $modP )->bfround( -$P );
1013    return $area;
1014}
1015
1016sub OrbitalVelocity {
1017    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1018    unless ( scalar @_ >= 3 ) { return undef }
1019    my $r    = Math::BigFloat->new( shift );
1020    my $a    = Math::BigFloat->new( shift );
1021    my $M    = Math::BigFloat->new( shift );
1022    my $P    = shift || $DECIMALS;
1023    my $modP = modP( $P );
1024    my $num = Math::BigFloat->new( 2 );
1025    $num->bmul( $_g_, $modP )->bmul( $M, $modP )->bmul( ( ( 1 / $r ) - ( 1 / ( 2 * $a ) ) ), $modP );
1026    my $x    = Math::BigFloat->bstr( $num );
1027    my $v    = SqrRoot( $x, $modP );
1028    $v /= 1000000;
1029    return $v->bfround( -$P );
1030}
1031
1032sub sin {
1033    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1034    my $x = shift;
1035    my $P = shift || $DECIMALS;
1036    return undef unless defined $x;
1037    my $num;
1038    if ( ref $x eq "Math::BigFloat" ) {
1039        $num = $x->copy();
1040    } else {
1041        $num = Math::BigFloat->new( $x );
1042    }
1043    my $sign = -1;
1044    my $power = 3;
1045    my $sin = Math::BigFloat->new();
1046    my $current = $num->copy();
1047    my $modP = modP( $P );
1048    my $factorial = Math::BigFloat->new( 6 );
1049    while ( $sin ne $current ) {
1050        $sin = $current->copy();
1051        my $numerator = $num->copy()->bpow( $power );
1052        my $denominator = $factorial;
1053        my $fraction = $numerator->bdiv( $denominator, $modP )->bmul( $sign );
1054        $current->badd( $fraction );
1055        $factorial->bmul( $power + 1 )->bmul( $power + 2 );
1056        $sign *= -1;
1057        $power += 2;
1058    }
1059    return $current->bfround( -$P );
1060}
1061
1062sub cos {
1063    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1064    my $x = shift;
1065    my $P = shift || $DECIMALS;
1066    return undef unless defined $x;
1067    my $num;
1068    if ( ref $x eq "Math::BigFloat" ) {
1069        $num = $x->copy();
1070    } else {
1071        $num = Math::BigFloat->new( $x );
1072    }
1073    my $sign = -1;
1074    my $power = 2;
1075    my $cos = Math::BigFloat->new();
1076    my $current = Math::BigFloat->new( 1 );
1077    my $modP = modP( $P );
1078    my $factorial = Math::BigFloat->new( 2 );
1079    while ( $cos ne $current ) {
1080        $cos = $current->copy();
1081        my $numerator = $num->copy()->bpow( $power );
1082        my $denominator = $factorial;
1083        my $fraction = $numerator->bdiv( $denominator, $modP )->bmul( $sign );
1084        $current->badd( $fraction );
1085        $factorial->bmul( $power + 1 )->bmul( $power + 2 );
1086        $sign *= -1;
1087        $power += 2;
1088    }
1089    return $current->bfround( -$P );
1090}
1091
1092sub tan {
1093    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1094    my $num = shift;
1095    my $P   = shift || $DECIMALS;
1096    my $modP = modP( $P );
1097    my $sin = &sin( $num, $modP );
1098    my $tan = Math::BigFloat->new();
1099    my $cos = &cos( $num, $modP );
1100    $tan = $sin / $cos;
1101    return $tan->bfround( -$P );
1102}
1103
1104sub asin {
1105    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1106    my $x = Math::BigFloat->new( shift );
1107    my $P = shift || $DECIMALS;
1108    my $modP = modP( $P );
1109    my $asin = $x->copy();
1110    my $i = 3;
1111    my $last = 0;
1112    while ( $last ne $asin ) {
1113        $last = $asin->copy();
1114        my $frac = $x->copy()->bpow( $i )->bdiv( $i, $modP );
1115        my $z = $i;
1116        while ( $z > 1 ) {
1117            my $n1 = --$z;
1118            my $n2 = --$z;
1119            $frac->bmul( $n2 / $n1, $modP );
1120        }
1121        $asin->badd( $frac );
1122        $i += 2;
1123    }
1124    return $asin->bfround( -$P );
1125}
1126
1127sub acos {
1128    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1129    my $x    = shift;
1130    my $P    = shift || $DECIMALS;
1131    return undef unless defined $x;
1132    if ( $x == 0 ) {
1133        return $PI->copy()->bdiv( 2, $P );
1134    }
1135    my $num;
1136    if ( ref $x eq "Math::BigFloat" ) {
1137        $num = $x->copy();
1138    } else {
1139        $num = Math::BigFloat->new( $x );
1140    }
1141    my $modP = modP( $P );
1142    my $asin = asin( $num->bstr(), $modP );
1143    my $acos = $PI->copy()->bdiv( 2, $modP );
1144    $acos->bsub( $asin->bstr() )->bfround( -$P );
1145    return $acos;
1146}
1147
1148sub atan {
1149    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1150    my $x = shift;
1151    my $P = shift || $DECIMALS;
1152    return undef unless defined $x;
1153    my ( $num, $current, $numerator, $fraction, $atan );
1154    if ( ref( $x ) eq "Math::BigFloat" ) {
1155        $num = $x->copy();
1156        $current = $x->copy();
1157    } else {
1158        $num = Math::BigFloat->new( $x );
1159        $current = Math::BigFloat->new( $x );
1160    }
1161    my $sign = -1;
1162    my $power = 3;
1163    my $modP = modP( $P );
1164    my $atan2 = 0;
1165    while ( $atan ne $current ) {
1166        $atan = "$current";
1167        $numerator = $num->copy()->bpow( $power );
1168        $fraction = $numerator->bdiv( $power, $modP )->bmul( $sign );
1169        $current->badd( $fraction );
1170        $sign *= -1;
1171        $power += 2;
1172    }
1173    return $current->bfround( -$P );
1174}
1175
1176sub sec {
1177    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1178    my $num = shift;
1179    my $P   = shift || $DECIMALS;
1180    return undef unless defined $num;
1181    if ( ref $num eq "Math::BigFloat" ) {
1182        $num = $num->bstr();
1183    }
1184    my $modP = modP ( $P );
1185    my $cos = &cos( $num, $modP );
1186    return Inverse( $cos->bstr(), $P );
1187}
1188
1189sub asec {
1190    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1191    my $x    = shift;
1192    my $P    = shift || $DECIMALS;
1193    my $modP = modP( $P );
1194    return undef unless defined $x;
1195    my $num;
1196    if ( ref $x eq "Math::BigFloat" ) {
1197        $num = $x->copy();
1198    } else {
1199        $num = Math::BigFloat->new( $x );
1200    }
1201    my $inv = Inverse( $x, $modP );
1202    my $acos = acos( $inv->bstr(), $modP );
1203    return $acos->bfround( -$P );
1204}
1205
1206sub csc {
1207    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1208    my $num  = shift;
1209    my $P    = shift || $DECIMALS;
1210    my $modP = modP( $P );
1211    my $csc = Math::BigFloat->new( 1 );
1212    my $sin = &sin( $num, $modP );
1213    return Inverse( $sin->bstr(), $P );
1214}
1215
1216sub acsc {
1217    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1218    my $x = shift;
1219    my $P = shift || $DECIMALS;
1220    return undef unless defined $x;
1221    my $modP = modP( $P );
1222    my $inv = Inverse( $x, $modP );
1223    return asin( $inv->bstr(), $P );
1224}
1225
1226sub exsec {
1227    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1228    my $x = shift;
1229    return undef unless defined $x;
1230    my $P = shift || $DECIMALS;
1231    my $modP = modP( $P );
1232    my $cos = &cos( $x, $modP );
1233    my $inv = Inverse( $cos->bstr(), $P );
1234    $inv->bdec();
1235    return $inv;
1236}
1237
1238sub cot {
1239    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1240    my $num = shift;
1241    my $P   = shift || $DECIMALS;
1242    return undef unless defined $num;
1243    my $modP = modP ( $P );
1244    my $cos = &cos( $num, $modP );
1245    my $sin = &sin( $num, $modP );
1246    my $cot = Math::BigFloat->new();
1247    $cot = $cos / $sin;
1248    return $cot->bfround( -$P );
1249}
1250
1251sub acot {
1252    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1253    my $x = shift;
1254    my $P = shift || $DECIMALS;
1255    return undef unless defined $x;
1256    my $modP = modP( $P );
1257    my $num  = Math::BigFloat->new( $x );
1258    my $atan = atan( $num, $modP );
1259    my $acot = $PI->copy()->bdiv( 2, $modP )->bsub( $atan->bstr() );
1260    return $acot->bfround( -$P );
1261}
1262
1263sub vers {
1264    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1265    my $num = shift;
1266    my $P   = shift || $DECIMALS;
1267    return undef unless defined $num;
1268    my $modP = modP( $P );
1269    my $vers = Math::BigFloat->new( 1 );
1270    my $cos = &cos( $num, $modP );
1271    $vers->bsub( $cos->bstr() );
1272    return $vers->bfround( -$P );
1273}
1274
1275sub covers {
1276    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1277    my $num = shift;
1278    my $P   = shift || $DECIMALS;
1279    return undef unless defined $num;
1280    my $modP   = modP( $P );
1281    my $covers = Math::BigFloat->new( 1 );
1282    my $sin    = &sin( $num, $modP );
1283    $covers->bsub( $sin->bstr() );
1284    return $covers->bfround( -$P );
1285}
1286
1287sub hav {
1288    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1289    my $num = shift;
1290    my $P   = shift || $DECIMALS;
1291    return undef unless defined $num;
1292    my $modP = modP( $P );
1293    my $vers = vers( $num, $modP );
1294    my $hav = $vers / 2;
1295    return $hav->bfround( -$P );
1296}
1297
1298sub Commas {
1299    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1300    local $_ = shift;
1301    1 while s/^(-?\d+)(\d{3})/$1,$2/;
1302    return $_;
1303}
1304
1305sub SqrRoot {
1306    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1307    my $num = abs( shift );
1308    my $P   = shift || $DECIMALS;
1309    return Root( $num, 2, $P );
1310}
1311
1312sub sqrt {
1313    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1314    my $num = abs( shift );
1315    my $P   = shift || $DECIMALS;
1316    return Root( $num, 2, $P );
1317}
1318
1319sub Root {
1320    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1321    my $num  = shift;
1322    my $root = shift;
1323    my $P    = shift || $DECIMALS;
1324    my $modP = modP( $P );
1325    if ( $num < 0 ) { return undef }
1326    if ( $root == 0 ) { return 1 }
1327    my $Num;
1328    if ( ref( $num ) eq "Math::BigFloat" ) {
1329        $Num = $num->copy();
1330    } else {
1331        $Num = Math::BigFloat->new( $num );
1332    }
1333    my $Root    = Math::BigFloat->new( $root );
1334    my $current = Math::BigFloat->new()->bfround( -$modP );
1335    my $guess   = Math::BigFloat->new( $num ** ( 1 / $root ) )->bfround( -$modP );
1336    my $t       = Math::BigFloat->new( $guess ** ( $root - 1 ) );
1337    {
1338        $current = $guess - ( $guess * $t - $Num ) / ( $Root * $t );
1339        $guess =~ /^(.{$P})/;
1340        my $x = $1;
1341        $current =~ /^(.{$P})/;
1342        my $y = $1;
1343        last unless $x cmp $y;
1344        $t     = $current ** ( $root - 1 );
1345        $guess = $current->copy();
1346        redo;
1347    }
1348    return $current->bfround( -$P );
1349}
1350
1351sub Root2 {
1352    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1353    my ( $n, $r, $p ) = @_;
1354    $p++;
1355    my $log = Ln( $n, $p ) / $r;
1356    Exp( $log, $p )->bfround( 1 - $p );
1357}
1358
1359sub PICONST {
1360    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1361    my $max_pi = $PI->length() - 1;
1362    my $P      = shift || $DECIMALS;
1363    if ( $P <= $max_pi ) {
1364        return $max_piconst->copy()->bfround( -$P );
1365    }
1366    my $x = Root( 2, 2, $P );
1367    my $Pi = 2 + $x;
1368    my $y = Root2( $x, 2, $P );
1369    $x = $y;
1370    {
1371        $x = 0.5 * ( $x + 1 / $x );
1372        my $NewPi = $Pi * ( $x + 1 ) / ( $y + 1 );
1373        last unless $Pi cmp $NewPi;
1374        $Pi = $NewPi;
1375        $x = Root2( $x, 2, $P );
1376        $y = ( $y * $x + 1 / $x ) / ( $y + 1 );
1377        redo;
1378    }
1379    return $Pi;
1380}
1381
1382sub ECONST {
1383    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1384    my $max_p = $_e_->length() - 1;
1385    my $P     = shift || $DECIMALS;
1386    if ( $P <= $max_p ) {
1387        return $max_econst->copy()->bfround( -$P );
1388    }
1389    my $Eps = 0.5 * Math::BigFloat->new( "1" . "0" x $P );
1390    my $N   = Math::BigFloat->new( "1" )->bfround( -$P );
1391    my $D   = Math::BigFloat->new( "1" )->bfround( -$P );
1392    my $J   = Math::BigFloat->new( "1" )->bfround( -$P );
1393    {
1394        $N = $J * $N + 1;
1395        $D = $J * $D;
1396        if ( $D >= $Eps ) {
1397            $max_p = $P;
1398            return $max_econst = $N / $D;
1399        }
1400        $J++;
1401        redo;
1402    }
1403}
1404
1405sub exp {
1406    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1407    my $num = shift;
1408    my $P   = shift || $DECIMALS;
1409    return undef unless defined $num;
1410    return Exp( $num, $P );
1411}
1412
1413sub Exp {
1414    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1415    my $x = shift;
1416    my $P = shift || $DECIMALS;
1417    my $X;
1418    if ( ref( $x ) eq "Math::BigFloat" ) {
1419        $X = $x->copy();
1420    } else {
1421        $X = Math::BigFloat->new( $x );
1422    }
1423    my $modP = modP( $P );
1424    $X->bfround( -$modP );
1425    my $Y = $X->copy->bfround( 0 );
1426    $Y->bfround( -$modP );
1427    $Y += ( 0 cmp $X ) if abs( $X - $Y ) > 0.5;
1428    $X = $X - $Y;
1429    my $Sum  = Math::BigFloat->new( "1" )->bfround( -$modP );
1430    my $Term = Math::BigFloat->new( "1" )->bfround( -$modP );
1431    my $J    = Math::BigFloat->new( "1" )->bfround( -$modP );
1432    {
1433        $Term *= $X / $J;
1434        my $NewSum = $Sum + $Term;
1435        last unless $NewSum cmp $Sum;
1436        $Sum = $NewSum;
1437        $J++;
1438        redo;
1439    }
1440    return $Sum->bfround(-$P) unless $Y cmp 0;
1441    my $E   = ECONST( $modP );
1442    my $E_Y = 1;
1443    $E_Y *= $E for 1 .. $Y;
1444    $E_Y *= $Sum;
1445    return $E_Y->bfround(-$P);
1446}
1447
1448sub log {
1449    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1450    my $num = abs( shift );
1451    my $P   = shift || $DECIMALS;
1452    return undef unless defined $num;
1453    return Ln( $num, $P );
1454}
1455
1456sub Ln {
1457    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1458    my $X = shift;
1459    my $P = shift || $DECIMALS;
1460    $X = ref( $X ) ? $X->copy : Math::BigFloat->new( $X );
1461    my $modP = sprintf( "%.0f", $P * 1.1 );
1462    $X->bfround( -$modP );
1463    return -Ln( 1 / $X, $P ) if $X < 1;
1464    my $M = 0;
1465    ++$M until ( 2 ** $M ) > $X;
1466    $M--;
1467    my $Z        = $X / ( 2 ** $M );
1468    my $Zeta     = ( 1 - $Z ) / ( 1 + $Z );
1469    my $N        = $Zeta;
1470    my $Ln       = $Zeta;
1471    my $Zetasup2 = $Zeta * $Zeta;
1472    my $J        = 1;
1473    {
1474        $N = $N * $Zetasup2;
1475        my $NewLn = $Ln + $N / ( 2 * $J + 1 );
1476        unless ( $NewLn cmp $Ln ) {
1477            my $ans = $M * LN2P( $modP ) - 2 * $Ln;
1478            return $ans->bfround(-$P);
1479        }
1480        $Ln = $NewLn;
1481        $J++;
1482        redo;
1483    }
1484}
1485
1486sub LN2P {
1487    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1488    my $max_p = $max_ln2p->length() - 1;
1489    my $P     = shift || $DECIMALS;
1490    my $modP  = modP( $P );
1491    if ( $P <= $max_p ) {
1492        return $max_ln2p->copy->bfround( -$P );
1493    }
1494    my $one      = Math::BigFloat->new( "1" )->bfround( -$modP );
1495    my $N        = $one / 3;
1496    my $Ln       = $N->copy();
1497    my $Zetasup2 = $one / 9;
1498    my $J        = 1;
1499    {
1500        $N->bmul( $Zetasup2 );
1501        my $NewLn = $Ln + $N / ( 2 * $J + 1 );
1502        unless ( $NewLn cmp $Ln ) {
1503            $max_ln2p = $Ln * 2;
1504            return $max_ln2p->bfround(-$P);
1505        }
1506        $Ln = $NewLn;
1507        $J++;
1508        redo;
1509    }
1510}
1511
1512sub PythagTriples {
1513    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1514    my $s = shift;
1515    my $t = shift;
1516    my $P = shift || $DECIMALS;
1517    if ( $s <= 0 || $t <= 0 ) { return undef }
1518    my $x = Math::BigFloat->new( abs( $t ** 2 - $s ** 2 ) )->bfround( -$P );
1519    my $y = Math::BigFloat->new( 2 * $s * $t )->bfround( -$P );
1520    my $z = Math::BigFloat->new( $t ** 2 + $s ** 2 )->bfround( -$P );
1521    return $x, $y, $z;
1522}
1523
1524sub PythagTriplesSeq {
1525    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1526    my $s = shift;
1527    my $t = shift;
1528    my $P = shift || $DECIMALS;
1529    if ( $s <= 0 || $t <= 0 ) { return undef }
1530    my $x   = Math::BigFloat->new( $s ** 2 );
1531    my $y   = Math::BigFloat->new( $t ** 2 );
1532    my $sum = $x->copy()->badd( $y );
1533    return SqrRoot( Math::BigFloat->bstr( $sum ), $P );
1534}
1535
1536sub SIS {
1537    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1538    my $num  = shift || 1;
1539    my $nums = shift || 50;
1540    my $inc  = shift || 1;
1541    my @nums;
1542    push ( @nums, $num );
1543    my $next = Math::BigFloat->new( $num + 2 );
1544    my $sum  = Math::BigFloat->new( $num + $next );
1545    for ( 1 .. --$nums ) {
1546        $num = $next;
1547        push ( @nums, $next );
1548        $next = $sum + $inc;
1549        $sum += $next;
1550    }
1551    return @nums;
1552}
1553
1554sub Inverse {
1555    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1556    my $x = shift;
1557    my $P = shift || $DECIMALS;
1558    return undef unless defined $x;
1559    my $modP = modP( $P );
1560    my $num;
1561    if ( ref $x eq "Math::BigFloat" ) {
1562        $num = $x->copy();
1563    } else {
1564        $num = Math::BigFloat->new( $x );
1565    }
1566    my $inv = Math::BigFloat->new( 1 );
1567    $inv->bdiv( $num, $modP );
1568    return $inv->bfround( -$P );
1569}
1570
1571sub modP {
1572    my $num = shift;
1573    return $num + 5;
1574}
1575
1576sub CONSTANT {
1577    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1578    my @const = @_;
1579    my $prec  = 0;
1580    my $temp;
1581    my @array;
1582    if ( $const[-1] =~ /^\d*$/ ) {
1583        $prec = pop @const;
1584    }
1585    foreach $_ ( @const ) {
1586        if ( /_gm_/i || /all/i ) {
1587            $temp = new Math::BigFloat "1.61803398874989484820458683436563811772030917980576286213544862270526046281890244970720720418939113748475408807538689175212663386222353693179318006076672635443338908659593958290563832266131992829026788067520876689250171169620703222104321626954862629631361443814975870122034080588795445474924618569536486444924104432077134494704956584678850987433944221254487706647809158846074998871240076521705751797883416625624940758906970400028121042762177111777805315317141011704666599146697987317613560067087480710131795236894275219484353056783002287856997829778347845878228911097625003026961561700250464338243776486102838312683303724292675263116533924731671112115881863851331620384005222165791286675294654906811317159934323597349498509040947621322298101726107059611645629909816290555208524790352406020172799747175342777592778625619432082750513121815628551222480939471234145170223735805772786160086883829523045926478780178899219902707769038953219681986151437803149974110692608867429622675756052317277752035361393621076738937645560606059216589466759551900400555908950229530942312482355212212415444006470340565734797663972394949946584578873039623090375033993856210242369025138680414577995698122445747178034173126453220416397232134044449487302315417676893752103068737880344170093954409627955898678723209512426893557309704509595684401755519881921802064052905518934947592600734852282101088194644544222318891319294689622002301443770269923007803085261180754519288770502109684249362713592518760777884665836150238913493333122310533923213624319263728910670503399282265263556209029798642472759772565508615487543574826471814145127000602389016207773224499435308899909501680328112194320481964387675863314798571911397815397807476150772211750826945863932045652098969855567814106968372884058746103378105444390943683583581381131168993855576975484149144534150912954070050194775486163075422641729394680367319805861833918328599130396072014455950449779212076124785645916160837059498786006970189409886400764436170933417270919143365013716";
1588            if ( $prec ) {
1589                $temp->bfround( -$prec );
1590            }
1591            if ( wantarray ) {
1592                push ( @array, $temp );
1593            }
1594        }
1595        if ( /_catalan_/i || /all/i ) {
1596            $temp = new Math::BigFloat "0.91596559417721901505460351493238411077414937428167213426649811962176301977625476947935651292611510624857442261919619957903589880332585905943159473748115840699533202877331946051903872747816408786590902470648415216300022872764094238825995774150881639747025248201156070764488380787337048990086477511322599713434074854075532307685653357680958352602193823239508007206803557610482357339423191498298361899770690364041808621794110191753274314997823397610551224779530324875371878665828082360570225594194818097535097113157126158042427236364398500173828759779765306837009298087388749561089365977194096872684444166804621624339864838916280448281506273022742073884311722182721904722558705319086857354234985394983099191159673884645086151524996242370437451777372351775440708538464401321748392999947572446199754961975870640074748707014909376788730458699798606448749746438720623851371239273630499850353922392878797906336440323547845358519277777872709060830319943013323167124761587097924554791190921262018548039639342434956537596739494354730014385180705051250748861328564129344959502298722983162894816461622573989476231819542006607188142759497559958983637303767533853381354503127681724011814072153468831683568168639327293677586673925839540618033387830687064901433486017298106992179956530958187157911553956036689036990493966753843775810493189955385516262196253316804016273752130120940604538795076053827123197467900882369178615573389124417223833938148120775994298491724397668575632718068808279982979378849432724934657607490543874819526813074437046294635892810276531705076547974494839948959477092788591195848724127866084088554597823812492260505610094584486698958576871611171786662336847409949385541321093755281815525881591502228244454441718609946588151766496078223678970519269711312571375454370124329673057246845015819313016087766215650957554679666786617082347682558133518681937745650014565261704096074688953930234791980600084245562175108423471736387879369577878440933792219894575340961647424554622478788002922914803690712";
1597            if ( $prec ) {
1598                $temp->bfround( -$prec );
1599            }
1600            if ( wantarray ) {
1601                push ( @array, $temp );
1602            }
1603        }
1604        if ( /_apery_/i || /all/i ) {
1605            $temp = new Math::BigFloat "1.20205690315959428539973816151144999076498629234049888179227155534183820578631309018645587360933525814619915779526071941849199599867328321377639683720790016145394178294936006671919157552224249424396156390966410329115909578096551465127991840510571525598801543710978110203982753256678760352233698494166181105701471577863949973752378527793703095602570185318279000307654710756304884332086971157374238079344503160762531771453544441183117818224971852635709182448998796203508335756172022603393785870328131267807990054177348691152537065623705744096622171290262732073236149224291304052855537234103307757779806424202430488281521000914602653822069627155202082274335001015294801198690117625951676366998171835575234880703719555742347294083595208861666202572853755813079282586487282173705566196898952662018776810629200817792338135876828426412432431480282173674506720693507626895304345939375032966363775750624733239923482883107733905276802007579843567937115050900502736604711400853350343646722485653151811776618109222791910224883968002666065687051906275973877353574444787753791641427381322569573196020187488474710469933656614008069303256185371886007271853594828847886245041855546408571556300712509027138634689374168266546657729261117182460363056604653004752217032651363910586988578842450413400076174727913718427741087508679050188965396356958643081961372990232749349702416226454339239292672783678655715558177739663771912814182246641268663452811055140131673253668418279295372660503415185270488028902683158334795920387559849886178670059637310157271720001143347673515418825525246632629720253866142593759334901124954451888445879883653237605006862164259284618801137166666350356560100251312752001243465381788522516645056739550573863152637659543028146224230177475011676844571496704880344021307302412787315402904251150919940878348620142801404071621446547887481775826042066673402505321077025830183813299386697331994584062329039605703190927264068388085608407473895683350520941514917330483633047714345825539212218204516560042779";
1606            if ( $prec ) {
1607                $temp->bfround( -$prec );
1608            }
1609            if ( wantarray ) {
1610                push ( @array, $temp );
1611            }
1612        }
1613        if ( /_gamma_/i || /all/i ) {
1614            $temp = new Math::BigFloat "0.57721566490153286060651209008240243104215933593992359880579380974513817815465490800204253092981070122661042345860326928047135164791203948800536606461037345976626826019171409439518699073089342082681115272480646074000944609116161372351563626398584332141045445255931180817233263931404751492309635212264793910054305740237382398192590888596562246640403865322573375012737199489715017633198800855368701197527793558011235088254024640367896574303325878962716603487225275031179009858629238810808962050809520067455082680224906319303106435436009045606431673108916223923423190040278230666282475597561157456466827358396496682688009413716435206727056861387351803685766804293903992265911164881521862119020073349724767061284456884146036064744615108373824011763115970464920405436675762580329560583963444735408889698752182450095288981177450041810756648904797158973756550829956482367512558077410478243365464613480082832289424563397700319507555405116866793330362135883160170752605749442210232264689088158360185921103264127646556897969599278787214257345380787716603498201211893231736410397999823581853335889687416623200928470782748962584090972449329340007986308871514443489019148224047176657926174820342428101040873335655351128574656826084318283564695216594139535628369271429378337337433379686429368401688003555383437291075008327921702577375368037141188793180192584269933204984701562366914646472808066790853122048683498800599411781499587969727029583663306226025115000498496743205991135161767514469471756389224373936066230391750756913421783083438469788717316008107601214614336979590782255534871391388706661385738263517698540450095680405658465359138468527276721930455900430554051676747620967652827322185928446094750362369908488472423752534464421059769032210076173483219811625852956410683745061833659140846327500871722721667257557546723325255507410472504211384839027461378721484740738636819445610096494307626397264106526101612141759589572570377325729026542963991795591864852272592134741602693673095062700415514589022024205056";
1615            if ( $prec ) {
1616                $temp->bfround( -$prec );
1617            }
1618            if ( wantarray ) {
1619                push ( @array, $temp );
1620            }
1621        }
1622        if ( /_khintchine_/i || /all/i ) {
1623            $temp = new Math::BigFloat "2.68545200106530644530971483548179569382038229399446295305115234555721885953715200280114117493184769799515346590528809008289767771641096305179253348325966838185231542133211949962603932852204481940961806866416642893084778806203607370535010336726335772890499042707027234517026252370235458106863185010323746558037750264425248528694682341899491573066189872079941372355000579357366989339508790212446420752897414591476930184490506017934993852254704042033779856398310157090222339100002207725096513324604444391916914608596823482128324622829271012690697418234847767545734898625420339266235186208677813665096965831469952718374480540121953666660496482698908275481152547211773303196759473837193935781060592304018907113496246737068412217946810740608918276695667117166837405904739368809534504899970471763904513432323771510321965150382469888832487093539946960826478181205663494671257843666457974097784836620497777486827656970871631929385128993141995186116737926546205635059513857137616971268722998053276732787105137639563719023145289003058136910904799672757571385043565050641590820999623402779053834180985121278529455415101923273972716796875156245586879771758718269365955450251304196818650938031303858435298686363516207327699768066589087224447928200866225489913933297068356193419225233102502064171141533503247429916507784663858387916406031641983699585984006546488792260661608729245292500798046966253372984581955631152133430603971500916257685064406302588048429495308259622428230270523401977496361785043061562923047970063433750838442657870812777765055245804574026584834104343691031824570657550964215122825303405704026824079955791184653756937399558705407245503084329065125986768185602850521404175676201486584618506998513595754222651824007565424310843506639303023675791283864525364442268496173451042322285580615818211343332191663663651426916459806509174002240402290641026530659304514091234316906333539039761913861815498584776177231773745183069689797687961364339064741256229235947371959033823110561543639353830741004542531";
1624            if ( $prec ) {
1625                $temp->bfround( -$prec );
1626            }
1627            if ( wantarray ) {
1628                push ( @array, $temp );
1629            }
1630        }
1631        if ( /_landau_/i || /_ramanujan_/i || /_lr_/i || /all/i ) {
1632            $temp = new Math::BigFloat "0.76422365358922066299069873125009232811679054139340951472168667374961464165873285883840150501313123372193726912079259263418742064678084323063315434629380531605171169636177508819961243824994277683469051623513921871962056905329564467041917634977065956990571293866028938589982961051662960890991779298360729736972006403169851286365173473921065768550978681981674707359066921830288751501689624646710918081710618090086517493799082420450570666204898612757713333895484325083035682950407721597524121430942470953115765559404064229125772724071563491218723272555640889999512705135849728552347645942418505999635800934732669411548076911671455813028066898593167493626295259560163215843892463887558347193993864581698751045893518777945872755226448709943505595943671299977780669880564555921300690852242867691102264527531455816088116296997029876937094388422089495290791626363527791432286156863284215944899347183748322904155863814951281527102068249218645827978145098870379211809629840943604891233924014852514327407923660178532707078811584944045092539519718157085780907690772192962552262890529967200510669638584207655081660527132551761150093619010182152039541621744474356571314026496051480322439134457528009739604967190734667398621127034770623094786463721777245551191609693349580116501538146897732947400254272699518373881294004390465050310091210361980535760952228835847669743267507757379848939356645406017251962513826671863828822629657399438626453078913514555113206475947913245582423662405126070382560901984614575152951511943211356814416716008974384391847402590826495013602834007260634108659796382596784136373377680857831279147106417370573337040146024737648200768231118490558678994106995743922457089666910491534089500139419890965785853368531985664042350494746329804481593573838687414276915611134778612290893976432134279879206472381493290546264824907766030881348705331723336407298994245656611424036824812873959790915799781062723446426357233234127834780836022424212901203199698485951429216878840715626887034517436895639117073";
1633            if ( $prec ) {
1634                $temp->bfround( -$prec );
1635            }
1636            if ( wantarray ) {
1637                push ( @array, $temp );
1638            }
1639        }
1640        if ( /_sierpinski_/i || /all/i ) {
1641            $temp = new Math::BigFloat "2.5849817595792532170658935873831711600880516518526309173215449879719320440011571202111177245270642830313439675271584928060821923187897872299047382419051407714877706637709789429909195090538916174018888831788039474837941321858645411399643090785375197001772338459009456944053756794492374579585597585179674255629788277342303908244843480262486476177584815605758512072553110971230816071780849788881500097806833030977613625019242141160511588550990118827059119218516656116658218146525723473922090837893539685181949438885982596384972943476477372870062514541167941067265109122200895264359146879706881090332595323182513862833772115297045105929773651123477119322088137378396355981597690986669776205753580929932151604999024815911781450614362312734398131562638214155851981996626439748302199986452462457953855125313088680611106070915380230497079775582427852984168946632924118537731025052638887185848026204906420983889401089245443403362540129958473807838153959208083537401553970775472153329227768035176719828981294382970877918341439096416723782335709603207108445520524820489725953913958704083992677439523357417468867744155305621034247542053834049067381062651986166749473771332141368588137842121600803261051984625149759697256509652152510827678630218450123284235823215983358532186222762388342128395782679202593945375516184778819579281898007193722908635995686666472901738758059961533111087477134416091651746860298589917405489412578156737071439005466554727980569541597398523605309831919658877736650415899607456983010590346999086310281901369036616807688073485295874875200896192730898004336728426431758229969326891997729365727971979849493629766039927445616318901564033213832906150701871783495293008326554891249827853248516958090952881854747420204575686195783951958051727205868237631974666042360351571383265906332448659953586584188803129625125027777547500955928437209118570584106256159701299968613543163628083412230815949199322528466660157369540251616007609659383652145974190555337555264511348670765096655113592156822083955";
1642            if ( $prec ) {
1643                $temp->bfround( -$prec );
1644            }
1645            if ( wantarray ) {
1646                push ( @array, $temp );
1647            }
1648        }
1649        if ( /_wilbraham_/i || /_gibbs_/i || /all/i ) {
1650            $temp = new Math::BigFloat "1.8519370519824661703610533701579913633458097289811549098047837818769818901663483585327103365029547577016843616480071570093724507999019639342272322414165036365074788027757760407005425387045947037548070012549126196000327078575312602462781280151598692712625156658037819170657049819111714215383017286869095002766891969837835648786933759294319175361858839873281361537111741600533650285988928906414670095488877382247112955736673406636533206353917604135172039112403028911351451318386134929257744182407526476030905279207782148560221871814904254471501463635842777947117746613775605839980813601589774035700341407559120370214113987005974964457642432794571720297914619514587500552129836800839402275440787337189077600233378591748197346154415354013755202065349536370774797232235307627711101354680926841172462714308267187960091741576168508046447756294559627846381809450570206315108346086296761158384244642331395026518568824439952885040681806714182600926258083237153223244690004091924289785349238396416174935955727240494968269552946845809079249430987481637553106432748666108251208300982355329789484892215432004091308695606679907614565980645118075299785348720969353909946887244982337583293642008796934457410799803453503483850708482981091242218216243859534562275020112435860924659080184472456129920697558489165392134942490568469061049118625847278188759772510374005405482372866324572279022748900193843433072546840353737116696380141809972845395860023313759602450865523234036167185219941330755467527227101818708316522432071143111039033646408977079116322065177143217096902088039791122531643152182505080070064980216107233022936897562281947118346305483027762905687068773859998729098422829244431023193577879027813710683659782342659876790900757742205098230554650080486761301518847811425034427156656362916286133872434789956033111453170457363680159688876278436595524346040307734465776554486790217538281709048389723416938426939063740559103208012841498333933171221169378991591489862913304490585047568388032591042048255902683612038";
1651            if ( $prec ) {
1652                $temp->bfround( -$prec );
1653            }
1654            if ( wantarray ) {
1655                push ( @array, $temp );
1656            }
1657        }
1658        if ( /_sqrt2_/i || /all/i ) {
1659            $temp = new Math::BigFloat "1.41421356237309504880168872420969807856967187537694807317667973799073247846210703885038753432764157273501384623091229702492483605585073721264412149709993583141322266592750559275579995050115278206057147010955997160597027453459686201472851741864088919860955232923048430871432145083976260362799525140798968725339654633180882964062061525835239505474575028775996172983557522033753185701135437460340849884716038689997069900481503054402779031645424782306849293691862158057846311159666871301301561856898723723528850926486124949771542183342042856860601468247207714358548741556570696776537202264854470158588016207584749226572260020855844665214583988939443709265918003113882464681570826301005948587040031864803421948972782906410450726368813137398552561173220402450912277002269411275736272804957381089675040183698683684507257993647290607629969413804756548237289971803268024744206292691248590521810044598421505911202494413417285314781058036033710773091828693147101711116839165817268894197587165821521282295184884720896946338628915628827659526351405422676532396946175112916024087155101351504553812875600526314680171274026539694702403005174953188629256313851881634780015693691768818523786840522878376293892143006558695686859645951555016447245098368960368873231143894155766510408839142923381132060524336294853170499157717562285497414389991880217624309652065642118273167262575395947172559346372386322614827426222086711558395999265211762526989175409881593486400834570851814722318142040704265090565323333984364578657967965192672923998753666172159825788602633636178274959942194037777536814262177387991945513972312740668983299898953867288228563786977496625199665835257761989393228453447356947949629521688914854925389047558288345260965240965428893945386466257449275563819644103169798330618520193793849400571563337205480685405758679996701213722394758214263065851322174088323829472876173936474678374319600015921888073478576172522118674904249773669292073110963697216089337086611567345853348332952546758516447107578486024636008";
1660            if ( $prec ) {
1661                $temp->bfround( -$prec );
1662            }
1663            if ( wantarray ) {
1664                push ( @array, $temp );
1665            }
1666        }
1667        if ( /_sqrt3_/i || /all/i ) {
1668            $temp = new Math::BigFloat "1.73205080756887729352744634150587236694280525381038062805580697945193301690880003708114618675724857567562614141540670302996994509499895247881165551209437364852809323190230558206797482010108467492326501531234326690332288665067225466892183797122704713166036786158801904998653737985938946765034750657605075661834812960610094760218719032508314582952395983299778982450828871446383291734722416398458785539766795806381835366611084317378089437831610208830552490167002352071114428869599095636579708716849807289949329648428302078640860398873869753758231731783139599298300783870287705391336956331210370726401924910676823119928837564114142201674275210237299427083105989845947598766428889779614783795839022885485290357603385280806438197234466105968972287286526415382266469842002119548415527844118128653450703519165001668929441548084607127714399976292683462957743836189511012714863874697654598245178855097537901388066496191196222295711055524292372319219773826256163146884203285371668293864961191704973883639549593814575767185337363312591089965542462483478719760523599776919232357022030530284038591541497107242955920670620250952017596318587276635997528366343108015066585371064732853862592226058222051040368027029750479872807946165810041705268194001909573346217594389367024932042269103436981246372011118526108426891029972031120210006350717637458240520384755519727993379761490610789498554422332600401885130363156114488684728158928816324518726506664538487759916257664287211124084206801676351710010294318071551519096164246090703940812921690351749296136400413967043104125363232703092257732796029237659774553709546911574214042423078199232761740190642451245487751686269610533369421621360539460424565414012853300781363344985673640670397734222981196104292553450160140594047954715453454840727173765626236654916664023300601326574407010783685846845231316046775448050040224063991197036221860292023886715071101716940029686875966350004089531621423342522795683406701347018590202836071676214774349344956359580808213044258646946852261";
1669            if ( $prec ) {
1670                $temp->bfround( -$prec );
1671            }
1672            if ( wantarray ) {
1673                push ( @array, $temp );
1674            }
1675        }
1676        if ( /_sqrt5_/i || /all/i ) {
1677            $temp = new Math::BigFloat "2.23606797749978969640917366873127623544061835961152572427089724541052092563780489941441440837878227496950817615077378350425326772444707386358636012153345270886677817319187916581127664532263985658053576135041753378500342339241406444208643253909725259262722887629951740244068161177590890949849237139072972889848208864154268989409913169357701974867888442508975413295618317692149997742480153043411503595766833251249881517813940800056242085524354223555610630634282023409333198293395974635227120134174961420263590473788550438968706113566004575713995659556695691756457822195250006053923123400500928676487552972205676625366607448585350526233067849463342224231763727702663240768010444331582573350589309813622634319868647194698997018081895242644596203452214119223291259819632581110417049580704812040345599494350685555185557251238864165501026243631257102444961878942468290340447471611545572320173767659046091852957560357798439805415538077906439363972302875606299948221385217734859245351512104634555504070722787242153477875291121212118433178933519103800801111817900459061884624964710424424830888012940681131469595327944789899893169157746079246180750067987712420484738050277360829155991396244891494356068346252906440832794464268088898974604630835353787504206137475760688340187908819255911797357446419024853787114619409019191368803511039763843604128105811037869895185201469704564202176389289088444637782638589379244004602887540539846015606170522361509038577541004219368498725427185037521555769331672300477826986666244621067846427248638527457821341006798564530527112418059597284945519545131017230975087149652943628290254001204778032415546448998870617799819003360656224388640963928775351726629597143822795630795614952301544423501653891727864091304197939711135628213936745768117492206756210888781887367167162762262337987711153950968298289068301825908140100389550972326150845283458789360734639611723667836657198260792144028911900899558424152249571291832321674118997572013940378819772801528872341866834541838286730027432";
1678            if ( $prec ) {
1679                $temp->bfround( -$prec );
1680            }
1681            if ( wantarray ) {
1682                push ( @array, $temp );
1683            }
1684        }
1685    }
1686    return wantarray ? @array : $temp;
1687}
1688
1689sub Bernoulli {
1690    shift if UNIVERSAL::isa( $_[ 0 ], __PACKAGE__ );
1691    my $n = shift || return undef;
1692    if ( $n > 498 || $n < 2 || $n % 2 ) { return undef }
1693    my @bernoulli = (
1694        "1","6",
1695        "-1","30",
1696        "1","42",
1697        "-1","30",
1698        "5","66",
1699        "-691","2730",
1700        "7","6",
1701        "-3617","510",
1702        "43867","798",
1703        "-174611","330",
1704        "854513","138",
1705        "-236364091","2730",
1706        "8553103","6",
1707        "-23749461029","870",
1708        "8615841276005","14322",
1709        "-7709321041217","510",
1710        "2577687858367","6",
1711        "-26315271553053477373","1919190",
1712        "2929993913841559","6",
1713        "-261082718496449122051","13530",
1714        "1520097643918070802691","1806",
1715        "-27833269579301024235023","690",
1716        "596451111593912163277961","282",
1717        "-5609403368997817686249127547","46410",
1718        "495057205241079648212477525","66",
1719        "-801165718135489957347924991853","1590",
1720        "29149963634884862421418123812691","798",
1721        "-2479392929313226753685415739663229","870",
1722        "84483613348880041862046775994036021","354",
1723        "-1215233140483755572040304994079820246041491","56786730",
1724        "12300585434086858541953039857403386151","6",
1725        "-106783830147866529886385444979142647942017","510",
1726        "1472600022126335654051619428551932342241899101","64722",
1727        "-78773130858718728141909149208474606244347001","30",
1728        "1505381347333367003803076567377857208511438160235","4686",
1729        "-5827954961669944110438277244641067365282488301844260429","140100870",
1730        "34152417289221168014330073731472635186688307783087","6",
1731        "-24655088825935372707687196040585199904365267828865801","30",
1732        "414846365575400828295179035549542073492199375372400483487","3318",
1733        "-4603784299479457646935574969019046849794257872751288919656867","230010",
1734        "1677014149185145836823154509786269900207736027570253414881613","498",
1735        "-2024576195935290360231131160111731009989917391198090877281083932477","3404310",
1736        "660714619417678653573847847426261496277830686653388931761996983","6",
1737        "-1311426488674017507995511424019311843345750275572028644296919890574047","61410",
1738        "1179057279021082799884123351249215083775254949669647116231545215727922535","272118",
1739        "-1295585948207537527989427828538576749659341483719435143023316326829946247","1410",
1740        "1220813806579744469607301679413201203958508415202696621436215105284649447","6",
1741        "-211600449597266513097597728109824233673043954389060234150638733420050668349987259","4501770",
1742        "67908260672905495624051117546403605607342195728504487509073961249992947058239","6",
1743        "-94598037819122125295227433069493721872702841533066936133385696204311395415197247711","33330",
1744        "3204019410860907078243020782116241775491817197152717450679002501086861530836678158791","4326",
1745        "-319533631363830011287103352796174274671189606078272738327103470162849568365549721224053","1590",
1746        "36373903172617414408151820151593427169231298640581690038930816378281879873386202346572901","642",
1747        "-3469342247847828789552088659323852541399766785760491146870005891371501266319724897592306597338057","209191710",
1748        "7645992940484742892248134246724347500528752413412307906683593870759797606269585779977930217515","1518",
1749        "-2650879602155099713352597214685162014443151499192509896451788427680966756514875515366781203552600109","1671270",
1750        "21737832319369163333310761086652991475721156679090831360806110114933605484234593650904188618562649","42",
1751        "-309553916571842976912513458033841416869004128064329844245504045721008957524571968271388199595754752259","1770",
1752        "366963119969713111534947151585585006684606361080699204301059440676414485045806461889371776354517095799","6",
1753        "-51507486535079109061843996857849983274095170353262675213092869167199297474922985358811329367077682677803282070131","2328255930",
1754        "49633666079262581912532637475990757438722790311060139770309311793150683214100431329033113678098037968564431","6",
1755        "-95876775334247128750774903107542444620578830013297336819553512729358593354435944413631943610268472689094609001","30",
1756        "5556330281949274850616324408918951380525567307126747246796782304333594286400508981287241419934529638692081513802696639","4357878",
1757        "-267754707742548082886954405585282394779291459592551740629978686063357792734863530145362663093519862048495908453718017","510",
1758        "1928215175136130915645299522271596435307611010164728458783733020528548622403504078595174411693893882739334735142562418015","8646",
1759        "-410951945846993378209020486523571938123258077870477502433469747962650070754704863812646392801863686694106805747335370312946831","4206930",
1760        "264590171870717725633635737248879015151254525593168688411918554840667765591690540727987316391252434348664694639349484190167","6",
1761        "-84290226343367405131287578060366193649336612397547435767189206912230442242628212786558235455817749737691517685781164837036649737","4110",
1762        "2694866548990880936043851683724113040849078494664282483862150893060478501559546243423633375693325757795709438325907154973590288136429","274386",
1763        "-3289490986435898803930699548851884006880537476931130981307467085162504802973618096693859598125274741604181467826651144393874696601946049","679470",
1764        "14731853280888589565870080442453214239804217023990642676194878997407546061581643106569966189211748270209483494554402556608073385149191","6",
1765        "-3050244698373607565035155836901726357405007104256566761884191852434851033744761276392695669329626855965183503295793517411526056244431024612640493","2381714790",
1766        "4120570026280114871526113315907864026165545608808541153973817680034790262683524284855810008621905238290240143481403022987037271683989824863","6",
1767        "-1691737145614018979865561095112166189607682852147301400816480675916957871178648433284821493606361235973346584667336181793937950344828557898347149","4470",
1768        "463365579389162741443284425811806264982233725425295799852299807325379315501572305760030594769688296308375193913787703707693010224101613904227979066275","2162622",
1769        "-3737018141155108502105892888491282165837489531488932951768507127182409731328472084456653639812530140212355374618917309552824925858430886313795805601","30",
1770        "10259718682038021051027794238379184461025738652460569233992776489750881337506863808448685054322627708245455888249006715516690124228801409697850408284121","138",
1771        "-81718086083262628510756459753673452313595710396116467582152090596092548699138346942995509488284650803976836337164670494733866559829768848363506624334818961419869","1794590070",
1772        "171672676901153210072183083506103395137513922274029564150500135265308148197358551999205867870374013289728260984269623579880772408522396975250682773558018919","6",
1773        "-4240860794203310376065563492361156949989398087086373214710625778458441940477839981850928830420029285687066701804645453159767402961229305942765784122421197736180867","230010",
1774        "1584451495144416428390934243279426140836596476080786316960222380784239380974799880364363647978168634590418215854419793716549388865905348534375629928732008786233507729","130074",
1775        "-20538064609143216265571979586692646837805331023148645068133372383930344948316600591203926388540940814833173322793804325084945094828524860626092013547281335356200073083","2490",
1776        "5734032969370860921631095311392645731505222358555208498573088911303001784652122964703205752709194193095246308611264121678834250704468082648313788124754168671815815821441","1002",
1777        "-13844828515176396081238346585063517228531109156984345249260453934317772754836791258987516540324983611569758649525983347408589045734176589270143058509026392246407576578281097477","3404310",
1778        "195334207626637530414976779238462234481410337350988427215139995707346979124686918267688171536352650572535330369818176979951931477427594872783018749894699157917782460035894085","66",
1779        "-11443702211333328447187179942991846613008046506032421731755258148665287832264931024781365962633301701773088470841621804328201008020129996955549467573217659587609679405537739509973","5190",
1780        "4166161554662042831884959593250717297395614318182561412048180684077407803317591270831194619293832107482426945655143357909807251852859279483176373435697607639883085093246499347128331","2478",
1781        "-1369347910486705707645621362512824332220360774476594348356938715366608044588614657557436131706543948464159947970464346070253278291989696390096800799614617317655510118710460076077638883999","1043970",
1782        "1124251816617941290026484851206299982774720467712867275292043701618829826708395745459654170718363182143418314514085426692857018428614935412736063946853033094328968069656979232446257101741","1074",
1783        "-6173136454016248924640522272263470960199559328290655337530202055853397791747341312347030141906500993752700612233695954532816018207721731818225290076670213481102834647254685911917265818955932383093313","7225713885390",
1784        "4277269279349192541137304400628629348327468135828402291661683018622451659989595510712915810436238721139546963558655260384328988773219688091443529626531335687951612545946030357929306651006711","6",
1785        "-857321333523056180131194437347933216431403305730705359015465649285681432317514010686029079324479659634642384809061711319481020030715989009140595170556956196762318625529645723516532076273012244047","1410",
1786        "22258646098436968050639602221816385181596567918515338169946670500599612225742487595012775838387331550474751212260636163500086787417640903770807353228157478339547041472679880890292167353534100797481","42",
1787        "-14158277750623758793309386870401397333112823632717478051426522029712001260747920789473711562165031101665618225654329210473605281619696918061316240634857984019071572591940586875558943580878119388321001","30",
1788        "5411555842544259796131885546196787277987837486638756184149141588783989774511509608733429067517383750706299486822702171672522203106730993581242777825864203487238429479957280273093904025319950569633979493395","12606",
1789        "-346465752997582699690191405750952366871923192340955593486485715370392154894102000406980162521728492501917598012711402163530166516991115122131398542029056286959857727373568402417020319761912636411646719477318166587","868841610",
1790        "2269186825161532962833665086968359967389321429297588337232986752409765414223476696863199759981611817660735753831323900456495253961837175924312108872915089534970310604331636484174526399721365966337809334021247","6",
1791        "-62753135110461193672553106699893713603153054153311895305590639107017824640241378480484625554578576142115835788960865534532214560982925549798683762705231316611716668749347221458005671217067357943416524984438771831113","171390",
1792        "88527914861348004968400581010530565220544526400339548429439843908721196349579494069282285662653465989920237253162555666526385826449862863083834096823053048072002986184254693991336699593468906111158296442729034119206322233","244713882",
1793        "-498384049428333414764928632140399662108495887457206674968055822617263669621523687568865802302210999132601412697613279391058654527145340515840099290478026350382802884371712359337984274122861159800280019110197888555893671151","1366530",
1794        "2250525326187264545900714460628885135841050444551247116222631411681549780530233516069957534394574922579290608180427520318235621123686109474343887857944611842438698399885295153935574958275021715116120056995036417537079471","6",
1795        "-110636644250856903590976481422794879200517231299540994715372334521128669716264196333811025709747746193210786820114369025849897345722531098042760530922656878891556664782168465095563132092311332073097630676251482491663634626858373","281190",
1796        "2525292668891404920279427026668969389456388249389889339455604316691573384284678293620100066924361693666444722338743839198221347931651916807651198800935942493038194104759967208073711284671045255047521429204396148980705984836743","6",
1797        "-12407390668433023412711473483696990726334795896412761472587854072142800403373577087021298541061094633377354326966623278849423631924808044397822651135905640812063181221280972334965193338438214107578486417026806166184210160001817890901","27030",
1798        "4708181368529492614110644197951837317202610608341257204206693195241245204360822875910613010433572133227831741097261618833206537519198885812254347219150482005543422997225440204041473518187636442241332621804718967775203938403965710395632762155","9225988926",
1799        "-1856110669947388268389361040689764027464160460436671923253131176853224087741924378432403442710398247642246902212818749685974336641529240178398124235555437625251481044526024910356131819016670047949661636539964662370375622630863327168696307","3210",
1800        "4005748930070152861935826766476856180706477227448622268042052745245798242539770546339789899546160341590069109467023517085578618986055969187202731878271685432460708841118506310943865592568791360294244451765746911808994482063783730693703607","6",
1801        "-11993122770108617858536443322964878003618156069559794803117809279608039120818829088000103355036592864877954563564831932363414886811786054601318517206937549605059298307895591515771731031691422489377098686236263367916404512751010916862894129855138281961","15270994830",
1802        "5646413644023523531472659729552574911763686780871700375627426663366507837122353997075891736706811337698248660838754243486778190331522785903661556458651175061469825204821206611273990250663647381261360659950519735730925202117606150672170127523599","6",
1803        "-8717064809960074651332043679796544474420053189621803209941154764877242575579617540185166306094852502223738126111951612525627910517501081576202264770178546608710937474005742969950212404955732473198451623024108934373488641161751158901712323446768306053","7590",
1804        "13368053158552172665210852539359893340369870651951497976111882891296650008003955172160792457229376320993686817755409436399268291095350295968657381088168219133490277914269064723832062615431730061224649980566693258603099340996988542301914519271322675688591","9366",
1805        "-3018240015081392087620978688311925380399983229633120268872695911807562982111154053235820050168829922189401964755775948260724401542319799566237745986245598102255191922935742610508280966462644022540839619861091049093129359799053781543195492373882916779852781709","1671270",
1806        "3168221108903401670436878558215734893322849540781208738796672473984272484396317849596978630272031342024194689871467916186898192771267438982228710525079886956295106197431401217357893460897221381410667385636049264583380749631776691121592016493432807733153743581","1362",
1807        "-1906502909997888166123201923177385833567729039151413143876271870599057704445939156075718972624235764657102074902610737729027517674632609562387841658709266014329005407533521950744449109807215808770201247724932231495252981632908042371307965561986133763291349835081839","625170",
1808        "12620737044909818561832856090355555624016887319605261762997689571062646649745107532482632213152948299491122977690702642386377706799989565320538433072623252159464115918057294711396625436506736041542403370025258477205101808063594056759977272469883621048184279331863155","3102",
1809        "-2277782962749544647786193093844710000584561145527861043594866852071596440076587235747852022645695523676561694856919437156369375152041285055935622032497285897151359345040171955786899011718891877065760345722163063921177723102631428767500963315657655779698470644544064472947","412410",
1810        "480867224771001697116513683699011649496855159878152931805742068517626950204279499281932407966997815888727039144001177194783700618900740782637516562256421883686495287564849123342843286211825800198337962453448529082007644439295666002669973893196613894216505936316966183107269","63042",
1811        "-321467133590936589398380572003196190798000628347443663674019204361034039315014370869884972981404460888272855773233080186485230316544246541168364468343631969225480324799028067015621769939718443419712110857965409179947456994970687005134042835903494172569465751671057323145801","30",
1812        "21954828084819816230623427376391154946628510524415478471365831649487269438565442138452375719697629394886161086199009071745032148355699097673730377547354631757000924688032674454208386076360699273002129330161098228962653466112032300056534380609136268347425707537340996253935261001","1434",
1813        "-212191008579182060478563837945461872287372869333130175188325135660897759482730035448146388858205966593247205572842356279587190846758925659038643395344396334821348861518596112373605365460921821915409949436456394757252173113308288776951844012432992719925522001911419529928297005743854184672707","9538864545210",
1814        "4535049728315239205021018362829154800039522502549714840841943960634084990270225535043892135549408608572877513963384530325758104248133124392959997485849004663162061065909846598215984547677506961593292880902830325868627515047168286738527241360778218692535254144583771935549805772798793","138",
1815        "-1480677868678810347062135814574727890490996459903153909612611791768134015908900253197632543925157559965099005581639883558125985134242978146873558628010545299879178299856929100217178891524159543673803785481607540954533057560554704283718320006046424881681934129216249889269701182688055001","30",
1816        "262668605206102118430195944953058387315319589613263079853189811674338644589939356676047394737425651989092115424635212785341601958591806782599345416064945113338168378810138043832891536358769753916174695061208903056552776976154338169779827414504134808495078925108569042616724875466536400029203","3486",
1817        "-3507445546375253318342667741949651253516986113349672180095468758145505921133172244240023419466113277413785800736682458966212429373095894934752041434119711352215164597094886530278429206637066302031501412499166263642447542289785723225064437740655097573160922291075175522551466276032634749948001","30",
1818        "3053985414762198703102020975667535181880294373982570358384262712216099013630127381226064671248970168701164841859177635139561220373102316273710980757409295039231697360536736113380473515993580188204966808323201772081668695321316425977503713933894934399351348240918411922863059423248058823372368106375","16566",
1819        "-89774288871630307843644520580598652632613559409413344033893019938898862113876440183943307724290509165539815764354254584162176909463971675189351699717189454959415964995262823025743004001915975128242680213612581748367743077774468104545056300803334120339537905702856877606246855867891786763741486735368783092237","303940201110",
1820        "2884759276650094186013272224259584433367802183845217969121665253431996536437413850716488026971480747022580182931940812035472752343853999483816456012957551120812115224018480617695398681102870535795278832047330709610952952566962354787854087985127629191094880078635434933391568643930279007460403102076327","6",
1821        "-104203435297641800915507592803138769896333127491490532929531166632958450572751929026080533456531382568746630428805839898606131705503722712567379172835813497582706763790263297372586804372654090271412733654409742029405259156130663560412599913831518338400385290760913773448329773518728092595019422753511174189569","131070",
1822        "56081693586889462275051589447484617682656486262465693822181172785524218338951641908717980184709547885060218908068568000776767834360070429300833879662707069272458173745883316346369196270761641325176838738153606175963665383216626724969370509650128327157862121275587666040290047452303859942885949047756404247033","42",
1823        "-5218507479961513801890596392421261361036935624312258325065379143295948300812040703848766095836974598734762472300638625802884257082786883956679824964010841565051175167717451747328911935282639583972372470105587187736495055501208701522099921363239317373617854217050435670713936357978555246779460902210809009009539232173","2291190",
1824        "6225500408881102797510043328034969466304425964851480796588968345345616249056035479080489691323529597769377127237910326626106353639640430862662758276037155818313797361926063056784629357346246050085674910484712331211955483400507386944492614175497512823803191508029088188458205305675244351279445756172428826335261196513","1578",
1825        "-2597900408162896058710572658949818524468219266931291835479959418094890120803707468900281325124704535527442908101022615930505212028093980003061941163493997948367473301269937152541190254678821414979800759785215240308773060568770563333752614741579673909061322574414039421772773951907381124389083123642187877106727324831409104559","374416770",
1826        "74333782054653889423743469213152655799143173104421681719953140280099539295673341397636592724713785397334372991707306538931885461890982390991361955644690371434518008626267056713898856350449471518770634624568188698734199687848634136577815884007651660546155376602651878199847336975669075705029220024942988801596540479832623","6",
1827        "-181215287752963297591649028858266281596159320543301504003425598637697061538040522696126018237983659329700516696739148841351912957174525785315414144429136591554131215700560325949733256954674813149699783016318415338094395743152932542723256830883071638325957369194991055197776850603497527503600109515880553950890416293966947553269","8070",
1828        "3046520935696995573515824242272701117754774305790084108344788704634810391537752246673386049133783821395486096239148415641722199623177956366799133314531041402021345562799254459886856571341036470447536143349353694309187363616222605848012490292828963830868021774140963898080446961958927794750526681977409210580806334902707423302910890705","73743978",
1829        "-5389479287019828859012200442343326192355509012382167997444221939499202340581932428533725140323939479047999949436682248001244967154707104896269136682582392582110484972339408510822230035412430172890925658230152084954903826143390588854318458185729506889417244891371385482660658834892925089584501833253476537744400283864408173645832802729","69870",
1830        "876921640102601760383168273082801324788039055564836901409234298002253995782034121007322766508391364331313364222035613918206397636790569250311165509441010945170696883786500399759308075144882184008248023432975450281978336106603716700929424210717924897570827498252893067062775490058563356806930680935200945326351681714440032862515566687","6",
1831        "-1388182772753480224524259543538426023294781120119840766760993950456297256537536762358314128425435909635398077811882539199662711471185655335807263675649242817349083277475519035786895107126764049407549359856712249235123586468592541328115099760851178135792262280438168965402774928136789551166963401124472782522202058296880872854126156007513528607131","4940319930",
1832        "3288574272791325983707258196648395370596305758341699742238893225440550539879910098002426164478935667759666907669638986616709137799200377160559516125812795457124992507785560266265469172399329113247421069765604730822849303184464007727659407021486798188927253139176046166778181472669994161900278671855691024038916704559966546186751742476372279","6",
1833        "-8490228242595311199119920565849009760217791023572798918266374350550246366225435144540895862924269878096526376888988896473279190955432566511173411205746712397692187887537779437718173050158951105924254322875775498449791436769311967741327489636985694874281460516082075289295724305512800386456299241907049270636679469815489106102786158831873003971980149529","7828173870",
1834        "25875593499348128833220857256544133843109199942629172661571575831896009002717379314993110538912888759932261638883397790531299585084426973156913779915781848556311555818927755071422637872741239146696501724664788565419251002500248036396312656297993421996811407029386847063439782311617981818206546311041142991312322619717261878315979033851316268092742291","11886",
1835        "-132959963583526612558423182997702578242047892089625754241387587021938995020766563989990185987109982064786136276402406795608981540644133005568845598733701027267134930895872207756078937107416605123323092749911180763524271940975157312568091937739185681066178901105822385739566191012805835839265694429000135104051173197039224461050364993416415448157657001","30",
1836        "1262786340764743836543650366374180273258073609841543441787781675688541881984249268360325012510167918286412665710785332495002470256724535451201518584261363284544626088778140281353574806840036456927404244048497506339393258324386379574048562781825209799628261527492591555049229691481538713081331723845462120745928641410668923052992096071568118182416749081409","138",
1837        "-4426158122015997656918081228667589801499808940970897924307086870252646875948841652238923933674489065786806672342321537012751674084821846458685429618183214525289276658012594349154796284960120596423811831259041410965661717709958659062147534583171671141655182016580199364485449629732438469542501153689517957906763036169880300563526920600913900347892060168169080524563821","231026334630",
1838        "158378960459933745353732673983377704622086418682000111407542005580986544145352715435649528631197645211762019722767661122078958064217177336225880988737955740461066476486097551548991955116319546315427449493825430704787375985685635337262886136570504412827663590990198163434342208367501561777230870660164068668650416176368657029205310526817730948015526849382859712055","3894",
1839        "-769496271232217048566454403451119653773906354892060893481293713784107951847019857310739767770463076516560575249619621659718985646126964039659981825577646225301189082287264198134787013538174178383958494473734277012586809351930781277345290942985827690026019144748986116319143448962997698665018791994064217597704711709748695678908604180272686662449493086793580970828093","8790",
1840        "344977368185637380304333898730614153996525419281986770522385915720139215554847760949870834982202303432587625997234795118134923343663265548204594792841811439303451011251437993865656175552484850145014444202587617378736472160288579064961864004275382498994164923679644150818450558302750885234058860046684381556954878962153088786251490404112306387375170114820733240249707147","1806",
1841        "-1888575590158495996440108433926788958859450449505827069108529221911737605491094242128180304966356165162290173103044092702582217959772178183606221416867140355106476799094653551123469827906281690228047604677838808485752658196964935031379587556530427680306537879720377479310817318584744711522453387072627648438830200285728320936098194745411999104991803838665658699650727659349","4470",
1842        "5683175611528696607080062478072501820703250802520963611339740016328084145085898175251616359343168595858108624167762131679681416977895536025003761007598913177279789271071880600643885302907670572288368816385439218707793725555126269377090645484013184046265389329719449523663221529092764894703813815177897980268745011893371899061378940560801507189664711567121796987700152153839","6",
1843        "-1863878995204859011995045341848156066182191846635905937518715320655775958174360523134990756922303410810482600528769479642021001218415879006164302955370460829146434807964717737195356935144151583424833154250047747433575584999029126775186293388721514970183351129809976971603227633930434923843984829580311593372565398574762880028289167635570012415606941367995702212211519561707046505473575241","866054419230",
1844        "29732914652005326308069038299046877228165947489285775398232775554326821654997284254954541436759755619186929818878962419179787160694148360451817997811323189321730570436847238527392204080083670713746104123180927363450365552833743396887002966884296848913481132785586683120732369858597153076261127975076406564459518695821788635361711463230165952857117226399909961673058570176834883271","6",
1845        "-5896752302971586952817880137553548037522952101385717862219869597959540714708235406233632122536677446426246295692119684486036889235399687553979484886254426933587637080451003789769961323298135257953657891721726854085498886402779031921701364180427128346719739542203363106920603436185256970882544947794042043194815113693600814746189709447821507079960920412101317160689806936226622681106017","510",
1846        "689735739192864790938869842200448522491622174065530181755219067505817267965189978776569968680399983089600630968341102178380777713630583306140103562669027424552384203113812510243331805627699809267240054344772443480095752572066698167737729190346148297831547431732607091731942172905658944645576233950351675859756367833238423434368779109181349635099892383418900862720397266277623080279813067604999","25233558",
1847        "-1310028313878567907114692495908315703121191755911131111984360268010339291958527789632281425987720420051304810543373657881574797353695997625625576992435250954115647419075544944736054854015711074248873686155373327776946348644567049129984144807198780234107278733859628145334731952115798551768101817753443736023659080348249001450368140437083583099127029870537678245450148711191680759214684176445667","20010",
1848        "3260606235276851195099437047707916012280054709370077715052312896011476347066426313756906147929976324639617487079038182718315850586912998681890040154101217023601746358548442139356881761209491879344448144235471024519916730280314406401197057752850694124801781598205133338269426259564611182358074475139266646779460335469700559434884540406485299632013344690352700850048935613677493364152601347478976955","20526",
1849        "-219310221923496085124921823808404575155215719519880743817262187359846676808720933963496193287320441312135479981243085410780329920552595746251706979452929554249409283208988310305531179583355242953344731240610852110199480347842183280905881606393786883347661578747110380713468900860626627458517735797820729101052113795399477516376330508889102106163067720003973339359314302749857453315620772706780505888813811597","561706691910",
1850        "5831963211655227469096364047000879785191356939385161602680269994879257080042651495531334620369899431331189755049833242684170073167978643108706019166324688447112227778757784415956624614457808484212568132108937476630344240535741105415692911032767001030027558118855855057952089777910853586423077582781831060554157175555095233630076901455631073919888020711465283697536929732899568769525385745975473155407","6",
1851        "-23306762036180954707239902303733182709978468580325093149239040051634149745258544596562134540960231525935441270015799599621135682414015382311861534597843756730857227145200622246509066032795288017634942726754909458234645264430702930529842456653080378544365215066466596986953601385998432554535709196917229350478449221406496506432167664369239478872974247912601100426306339312275565475595669640092216631914518917","9510",
1852        "28122967090195287981055019523746421026174611154606607881410834553821571195600077156388139504050817072318352135336014937133504488344591033637378875175488860842142200215544952457462251551728305545862703146471647483394138505742989863336231595764600339418866218322681285266901817977284569573044604541872422517736195703069383015638214571378546092320911382385141634160156861783095114673655212399723765743717305197851","4494",
1853        "-3721822981147963543473628721718160389675395897437503779226807860482566691214705959016239951890193486758232537597591642279654459920857716018601998378015710967458002507660860624568924448328452560808376796224541899685756821793434321138838512925307216420554703983171171540795387100799809353390519107647717976213032752752815159815739322121041527686676774334782328158773314452206678687770745486976844581551304048486668867","230010",
1854        "11947011202513683767723104489271471097868501873399758887514336934592822829976646702161909309680237660201998669294617552009638464244432709441651377979627696796385888263305915811924616794303093449550416640107641304231571552322204363229110497299236034396583681906657985198177302972578027688117382447661848083476283296893356637990999976607059658190271101386708876477455968880271181353449021120174672279118531405266405457","282",
1855        "-3829393826694851442367963605345986414644123194149055126713238369521199883032584377556163984602775965899294620028160948209504229098195604243661949859717087550540749941206115020546490009175884700307199551574107407855402424207686775155726579245497021834551699678553804651177095786249593672236240744422195553711332954204942399883879329051596664648871055498926334411033125995310873151284771962309916237507813788879633942447325937291","34098248730",
1856        "1808383072515392862130370135635372000131626793093148880779287649928820293768508592707521055628027708276165829705670838064031068574527612930799078981578676791129015219438793213734036863426110360623900757855812426023541159774905365518756692026508213803668789819570234151094586211640196633499617063941661470860479560003474875427679061750173843509685528353151826146684405535206504936136026796030365288231301797295663647612503","6",
1857        "-2038920955732166704563544675871539298534664548832344594940419211031004686244887980129990722833481906996574995006272283772989552251711645339529824674672295721387444989765197094356625686277258759887529063884063441729464466695304817154879753936013360244814470171255312860420052459456595966677396165132301906628818307744454944935388834925108370678592735494494812411438800013100442226704887804736057100084506757161238437561170752083","2490",
1858        "16450740261761201488286514086666341385844412580813949724000155038308870188720719658119876451207931032748030811189911338403495201879917853940007778344522354341186748832271809821828409298207968739139227014609598106703237472043228932298563800984396865943755795758369082954888437302866500408097109163018555655697695594114688192421141599313689189778427690359665256283948257291245464556741556264558183330312680675299274129472749543281530376155","7305236862",
1859        "-31404740573139387013607323899060788461129734403267877740839079468522117216337083625061725107066347029224470970177830290797457104821060170534363148276651160571699094881049180803660482733644591278636224534484481850814754420273709821201199908209760469191567840252247706876075856090781209456764232767086688491357087469900906692342923336064928179724546370791109956974800789303951709599892246386582961297095349501749155074025993332911366367","5010",
1860        "105959450712139358561614326750678748291627849550191512662380413294068107407247886566713525350838198788534238121678429166419664041236233803606051817627731575939553790154021616196778501189528768078400811107009317850171384096998455082824610437919969197421385189475358775596537026092374519645052975059595487888618614275545703577191947383711761849232715855501325445057452043875112980841528855594013764288125280661970576229405994478822053767","6",
1861        "-110968358518556750910569006059104612436322304793290197896476208195051113971187304256671882613201063543710887491006989203213468862095728086421327649757926706375373160032479060881414285316805050049998250096629298691714712337731033133063244894833650574221421335477496494331790029997128928883270452912343595756835759921082499027715633187483508227173954421007445739076452215234463953952029377612275096135644427285378787750547213980148328191792731369493029","2203871994870",
1862        "871667613876289947152259213998641324359438367113245978510277002562210922288209180150151793852558736029189876363096736208315035600167649273464678134646229708161579214485318416128870291425640987624153304026736765522468963515580422620761261195314156252082303726799452238319359220112879725078451979791840912313059346932550771416018721541922207119333620480831014488368439858212426770461809324187079997350701358985200574078187666941445255865364959","6",
1863        "-139969199374294532188729693702577000000465664934371179035889932037927222630615161871008568258760858093560212970955512057983700284904139662988796471588843837914664871339017278834524915750717069973487760152722966363103180527139199261072959205960570607504345672675033426810959288741490491658133099197637151741354653719448795255505120643394167214213806837785207626077991226383344774422283805854628774864645574223075309564129905544107459811288704862611","330",
1864        "999866869350557363722458479004395830847953169240152118030482972984856254139343285563059710426142001310818604300278555911104883557186325333153917124330048368249458087059040520529706556642732809436335612822769043297715578151918469336944179483860946506334004186007324467864842978559250174108263899703085692274142353761805113570952308838426945093917953996642662776561930695988488695006166233701104386823791719249365863288962981240951731129450241042813203","798",
1865        "-19435667946412236600006355752680318177378018998686558693047046847185922842569405850478643326624807107024135785299799674685846584030230598008329461654964640062557602833521189158599030670653044302429200521250792329717629531968392422269289900710918148435060503372045597753741275339144955223007618411222546697775810780747446478661790743019087588298611993809224843559412437328287256779856482337754827908028601591041648043989896358637724179061858684905231780173","5190",
1866        "23574815095404701575006138342677814737611117002189168884514625588124591088321558700476292561495266316421750730681754426899890383276767887403267480686627676268685966493677023072554534770488586483642062259843015737466820417963854802166458232968082382836216212063844897312505105552626378830504014452088433926239997258280111620289793360632141212087877945966075033938687737132829040482438463044913103827680386355585682702291192692673999705957317036687127526488661","2082",
1867        "-1946958244926408519819928540035312319794025231154312283846941445533833885684884185938663359749293528922631351379885288320544746558558870084671027543542465268801245629638273090846490014696481293058172945225573249920181321736453062629022981688433233577329127459928765113606300176648832109875139888454987326886636880988465878410619882237596781152380489370724295257792104332503696582825373077048086525583941198173749084587321884232599450628448598711165788748886648206781","56213430",
1868        "502172089109637973079905157308517324166048351810351159422375337574332573514046742139861365117790762795398322755462389540876498686325519933386024542998456467504352201616875674226702066988781736776780983675270711591178536297327588046501114273115924942240089195877277170211145955479537604561876530663522061491492548171664129861255728492613537859772665546205365872403842851465588616941966018356565181922169705164734493393477247978508055348618032704534473742084629023675","4686",
1869        "-123595612458163473336102111542838674764831835495331502431774135351050074096606458390189947949855899933224956382090454659506523733852842909886351324966274684823091833542355636103776512914842749831251392589846063826895163366115470020631076054800071496634622909233303706997995125130120515859762839469193960001284685791277467069230539953616534532549898448827387071275667857164787765007253139087944551495091418632200084873393873944046407490699550148256008464127109223283763137247","368521410",
1870        "44586958800126786790452414833451737276485269217604041493957832153847748776691019563443875873891685730191225792159808508144405955278309411839337401759278396174609624192685884123993968501104767184346477939067903846512439522927159649898761965816525286791728854346485686727420275292569511311545587705052888988969140031330003151200601432709893613802729701685072968681090703554561525613786761794587069881004247080811531298610122589066840923344874780872870506826198179146216889","42",
1871        "-18249520604687390491188873019557839254274215818941916557220148305863171323999328631073000715354712349442824804877473011208435647924577393976121564194070391790525169749128982828132077146902780721841387850162775862442199934474879939605150425883830706636702361084052335568869843628198626221242494076538767614229186803516353130263179834649242288787960314531149564506544752097614513847899020629776701115098262552562151629993531946710430672386743677617308205579019072302696542508379","5370",
1872        "23698146493089350432280213985162669959262841111102470546871456237323693618435911411642184657663645971277792033537906347762183635901158641375259341789919760310568146035244593237879245436647800609988327386298127942914981131212727307691555784531721623262685912374782653802898934200911324353096490548706043962852465907958751649549222286885339380268684947173053181005751121488921877337975992137735046528729900037937548187674111044469106308867226010497700021930277136809473041010432321","2154",
1873        "-778920925563520089023277887716919839298513486135455047496265119151473505464734957690937635048394537294415972122421086877478533967411348994474171388256074145035237797802154822761617528172700015731068657547092679836095543042179681364991968425138238530710252585175316218097390905698054739393609117605669741881868029805237983772228322424502491583995084065878787104609542382045050185577023575049548616444750140539582676309391395859436260664778641733708499657047550739325300732722680932200039559718809","21626561658972270",
1874        "715341102258098940319912343996707811685084551356244645540933308810413928160891724298276851619649551829973697033587181069271380692407132317161739109284011945227649892560360729074672778640720680874112214109997197387525352062913504682830983772725393570468326491083187300331970668449041594296073065276692887678729691114916079642579900288026881432235979137544823640729868705093774645241454622836223176904124630526882732099120398430089179339282190314646212515622936459121100998906169917551","6",
1875        "-18399470428456050908014333117886494595730651014055636288722717429950917809500136272557117402977052839042105693597404693257617634201600598895359708213302738206862242849606975947076539483777341046657860376203296993121539300815161693660027527850273081527494673550440357634430596208643360358204371232957151165657607509017253243152916116719712350843286172274487990507840124030754833011930721808729659207188165402273267771908929719004881279785191361391914450528631012111465559730900716078151518537","46110",
1876        "20813243679670197402657864594700362248333105628735574505407748752194329210363364029904045498683920274785904391711886368296491226342074760936891789529952828353712861964615536938434038079631985205838319830159148784466890765577478710750411239086707279314113340423502528151089538256774657292895067372560433868261285932669643784461912885034984352745012134324851266880718011466493208024572308839966018431480701067038072465848677106468318507151450224167148317101357641051334844151154035628560478914687","15414",
1877        "-110725230879881321879568054811797526267115166186538309029826133137308115277324207270571040263441920262386242076342056949047866353534514266980024412370434220187359560062050493440457221035440089226046696014690219122751069156630764313360960761587856716000690131098571105908045239988081704908481660126863809763243757816063349431021332536632587298653399772920535661074013239670953426088742833252653499763654446938447050846790045870317112050534444840968778815676707175954631059728590843972403931253118799","23970",
1878        "1054364681019267121830159662373420279818225480653781839149178920561555806813202030269003305294989784104213489890426332509828119920399712689187184649711689415835280384720994979874232796860338473523737136345243289020296312370876276856303158404225225651223912805235434755784151235675589206423964860803286777387807480810963355146229723274524146948489417884118537852386976826754447246469153480716621041563108596300892747446848622658459755002905969181236655822204297783975692352381312064835994355860201685","66",
1879        "-56868988747972506254659596935683883463716573524427146597376366171981802233487630537261600452515928180502785214820124586232687460333822056190451678107485446400177877446417556715711032055230455854240595868388654451080202675047744244991416192115524025084813812227930163123287145588705727429456779082968390877093897979160965323905993195706278985485896440595636775123847679345910845046714163301325224520710766325936895600819434580574518875433796519293151389916795779147962584026929392606854485528272956190655743","1018290",
1880        "27233521998844117466571164441681703778822231688822656418920849017523981696494911618651082582707143082377020546888216407884466480834780583096823087696141619926246646927066679156637706667759198043548965512948036645725109831717082230384665542660596536212413188407637782216353393129995452568219713750918323446571634462816666446106121343241772625570723451925602179037604045813434065564118195563426039220880907100156027093090547687363620867038635244235643402355449232121698090257432456131056504444320853157345201","138",
1881        "-21144886325968847421136862640933605790730878065158231117894812330764874940817232709655287449203776974910988576503561661800249278991319113309582335281592647442495030553208258548341339318764197461147432848991686338242711098654931144304725148541437215370995347051713517773502906169174089944259231510484935144428438091067266131242394201122670888252772064856423958372374082615970685597386991737179844349495505700142239123902982043714578671787314293339023048726217119025452344630237832665698519639606247957083240801","30",
1882        "4202152324841490851582198491734154054898910989229835692416027795481485467411812516892366095348619673337490550280425989691409092856177027846893427878738640633197809375983848427036487478678839362355759610808529212110936501477008097429319817529927993776215011685533261432410077791728458394351343228912812829117361997697793433448908490531697660725893829519644454139967318758247135377025444352170089166553939801451419735179550891309557305319624844923751198154628001739638691998206529289924753014269328705445340396548203417533","1651635762",
1883        "-585016197049150516429606755440285431916223316191798033687181061470667150509536113671258201054717866226570169976401850360230659699699301963057743166370714339259000981058060988677132844370590538962044247092297537672905410985648312911796535010450651733583075122028847986584025427113069993778348629306554503407744980050377057969061170508810412927919930789637438916383763636475443274955653725558963454824998976936559504337387796596709727547538729277633022545136032386860825147973153846428111580196229784765036789760566304701","63030",
1884        "78631979609289345055469501119395008379118096513723589634568877004799927202495863945986283823608714469880713486237817037414979036031006235383025032476262083964387973243898060840722936009590301165543812264089314197515688838727901018644687053389435289951515626740339208758682375061416597195951759290458697298602770109818453665947288873748972700580769864649544946817207071904078063454972954577208463927530056857030200986057709786810533019838606174872050746460837475834260566055441410395805129523481902672214740177851314801113","2298",
1885        "-110754162756402172153014550168566972984876925083693264831143127544432717811904992984077240231351850620045442158975752476209437945003921292257289655906112174208294630503691959241009346814626779961417257948400146389011140929210993075671069409245612465372014808443841918963732454979635619634196765514065569765406021995951240364020270070411437377187081353045931837302854850452592893523482298799136698133585855070465304605234654821626817847160151252277371211732021978740981118935331815732788765628347320001270460720535171080710693555387","868841610",
1886        "2879114883186609620826723532068276104139845449825894019948530676086618349883233245139084663416929458760597036667413717978530021472148264070597286070616857893183807259158771980877711836973364470227289919568822358886590073234578910075775184060367188577311832614887005612354400063672439629984465735358226535865532910705094242132700635264096296205052370626838252677418749801518077664289390262828965585623070897002640768910479324438878660938953441067263764560452279108930457475258025425810584201571774449957241464340109525100976383","6",
1887        "-21299114788094757397153706724714536413292077775222009197296487570276368549767668483739905654759629794289722059088574391474394384116388842393472198779890845792835783368209072403070166818305644992935634095949389005892015797599060514475400351300665686095027668428158520105490245612578226308332064019874832238296603858849991741989885454092253889660892754436759635904520725519771622791567806625636568537908996210911982475409155309551762403125143065067389753960600797968698009862634474184891040512288957749239742206998790145616588343479389","11670",
1888        "1039554412204304060999050471842935871664238236129525031050621197904673190252061816306076783004652192616133527774723493749394461545816137751676691165637363944930483826319036927422553141330104483296290342476507053344206817648719129001417301179730182909494636606795995576156560517060235847592695182729113023238975218843912996538779543747649999358252604590089345118921257104639445913236211149024538216358616556115899687337155136781007563112551684977465108319504450250735437471011340987371661862741413337000516697089619657504784045957077408711185","148218378",
1889        "-4666956320257517816926480178489878702719387148900232522527912339796188714409776663856362141719980502875970806945762712937180144475887677992558175799481301586019202060028185307056311432687319799854135340241062231395259919629305490635354239890090217072557890791528689479267987496699263598146171495893203013485575227547175970131771602894568516495986613037151863666577390670765165933422564746168996994413125484969251071189411262200344536746879649016328422523901374916569404613335757056186544409616787583076594896302119414868190918059336879970513","171390",
1890        "640808912350356174378073333715330795309912947431801296035514496305956924194050478724997308276302531933458844716623312895516343306208024932183790655351844102967791249508631724180677304236114904329953319103713608769803340518450516179464461284806009110502342220743568034682133273443656465906102058462073733441986152637932219455262690498783033366796807895595451968974170018589790265173871732540219203496470134837096254721957294282692027321105998231580155273075570903000709400332769684568638209079266772432973508791091418552966827434770217964847","6",
1891        "-98872162196656114876533219362632891391774839342005000436956938726481191695417122725991435991920485003784886516420949151408274403750897124053086048470809705387715020563676318491207223028465413318847625119271235350352217955302858427960756286710348612322463571556662841990100314844270405671459605112346068263402145070889841697007530455932282600691314028702416310723147417755942945877131986599840610369192441659461647425950295710718456957848137284733782503430592302360661270817820862416426647830439342025098395116603056716015807642935747540566951812504354652379","233649143825370",
1892        "10161900313215566319953421591388637402328180425349211048816698980319072816903521145345726221163084105611450343443221137527690722856383209708142471452342875346539666863558482835019156804751669116932322755085721460877938790181414341174421583630169802248409852214724801939766701626707662229904797008779706437578874336997303238218116573968772216970036489561373747407630110362020821819680854208929607112013855668348426244596961480597381595137016528270815603652018778457314919825182032768360703403384404569260876760458965758437801728169587522360995301639","6",
1893        "-63783639210289998710088973624588233946566075008131539459746377225328551512763402330278905606150965942861218654753440979855135225080751599455762303374322801376784474875903535720147084146967619965166480436166324915138652313161158459147438440579563255414402868679468643443725452316589751257887343519347376962412796543353378510329227640325960354800434372777469890448460017461964152127411007649663198873691188302850776549851095445319674329831032713432828366945325817352958745749412125525231627662590308189489017131904402298383855303749523779093820302379681027746367","9315635010",
1894        "1174239835539994552190497513182360535032387583605057260007722004883059102557958144704940674287492304226247212849059989833910250445125922032127665436388239538661698487981932513000579577766999325929606223291472664224576178291713742872723127632853994169543584403270186239889883934311766077803787140536031309605772620994720158177432481254737877501095184402972656917745175226661913130878960138618730840741958163788195850901560843232573640710143138001159729682511610291187234554610346621711370121346135534129198492308718185558865166313048808914617330314355609497","42",
1895        "-3459038918426949807320150484933956945853243159595350413053075601122470521916955664973829953506019998422380872167968566067266857401026670310184013796611289015884412134623896705526558633772056365694507697389600886696154015160862337140032194111117938579681157914085424486545326355341017764587852271884318010152210498883089816668030790992239717499457275164513815414692072356653908445653642502944394623033131400812242246305136709927393331111317205972671896631185417486747798871333208226992838344653421293924893416433794558513589011611309370041563779286876725823001","30",
1896        "170003857441113621055924711218201581479060969973774823543550655221015298985029713903581691625806778175064153013683952617568923683195904825491671562597778394848940698100060988046881674679194048284424046096930240621809667552628402886187489064670995742803216624905991737595831860741767548013712391467620520309900241319253287943675129083092747241525094828376937402590205346662298390862562643436019735148452642958172300132516340448730506554851256643796686115524433824609798825320819953662000059899761347431172299049404919145053054118348793890650129877368554102264664237","354",
1897        "-31826882167847269036795264135997035261326106742393531631262513546951621739538129483912871462979086464612269631159787510831358937051013784806855551647371523568233230776308094926266913151377871925647720334390123342313810138336493293789276924639127443823386014514502682882716709333886588047556737429630427737978891250291212164857441066067952971934975332936367526583068834875982474415309122221371543750940258579267452360489050995048217629796997409692793673935051432353130414145045187876172429836566644615416398029958336369046123993446514393441794181039041943485420490171529217309","15755919270",
1898        "47002375233508680615492575022874518250580216595322188118729388581779057270357890909757660925961204865010764157357677339640679419801075529375841405834200075859710151539426985792722765113963198773076233038710305332376160833224527455197843754378242330663125902354233318353062799203163603734366135826055833545210943663972100422012871119726280972153082552774282015316310107675884741591238987329522432070124995685585469794166878230612009713252291615918394199230644202167695522574730295764175556463982949470912949002495605091357997599244798045099081157657094553333043809915382015","5478",
1899        "-1104074382679040449322457640474906559564674033012018040552761693975196623398701263395843134574919063111979334979001689856745801791431844835192880349897217288293513020820761816536292340772895413243616457931467439624320619545324679396693426479472670415005491693729301629890910865952928640236591593400805258179861429283108782677380853799729252550169975507441841656878779871922009807575659278905966164200037561186284335874588669047079319694680280463911923325143499841983743646211629813016265241858462427586701167961402272877528256110590588675820001821779929699856406082793746601","30",
1900        "830966056816864165219563221691255172504816245427194942775778460297239574008254673723276280782798993878568376502765275874049468633347470545478154168815827272595416765071482194437630416545142455883532025709379623773918422191234979559992742090527602639659473260751398640580101548381014113971979542275363074811349527026580264739170503267212543977947471780118824749687917395812884246155157646166836076504356150965742303158590029916282708649948973023160673477590749360021799996465506161348017383225242179003450331850903647117458837892565999818223864269663453594385683125383821474811083023","5213334",
1901        "-18840632174653835864968232583541556084424607654449238363350816356454780409921315489669153042640423571904583615446707687427327682018540206290243741793046153140931788677169249364676259458302231885326120895149522010431725107071591063425627272025446059167464033814661996655342902404678020463898666363623366708637084399729725592271785009331862696316729612579245598615561813162507665448387066401343329999818243511962109583265019621658883154231991893552479038765074765515515428784392050898556152995633956321371743766583910838111997813677461374152418833832478425757962198139801618258600291701","27030",
1902        "177948829071244383748461224300047737954439245366393998814954157666837290036442548237545354559654248140762601235025124016956688571054331709927511568833209554328300942088687575803158770352778902946800416194092364129368960809507136675944946524355529271391513095831833834095214141413201872443533465259338559071060276534081662109185490926535115165429017697889522123426359409741271701359896235732490827288153220352311977810488328925043422296933754175073177948425972950246944493271200423473793223213272148300689418662045121140562185310116950518967585206995379939580805302928813353440147414740763","57822",
1903        "-6126914449391766463123273102284652981370509248391971126816109275821466096976801669348238710361647712578179757979414761992440067232917335164055107920780363262942782633482006999971350591213671109122331751122081373688391335194625200099100643843244514246383231743873282258287892336340426940736781587230644995704054335003998425945469627545685623200550970324534080404065484535299953386337281009271537528001339179644699195404880031646796095372570028624110131487675467842335793522121252328612723178264941963550114480615755660790108038424472918326815297861578110970433546154844146996383907171442617576638221916277","446617991732222310",
1904        "370417641649773203811105792421530317349989434553176611030631322780037740921920903496806496686319103807921979015270437429241936687918394087764863038019578527638851481686403847480960464641545840936004289242306394237365086284210791291343989458767793541220855587660095196948191905734357956461535121861354233393019487098593606068830314737953863466706245485662337370451037774820208643752612698330292374600196793018724057118532300846184682702881492970790966220893763118615367891028920559173128847556946881708070136133779182213573320504698972245409107308503781210793753700634159761806826499289415831","6",
1905        "-900309704889057549412703055487959300250691255243632314995356843469409197586598181836315775665695203810817228257767800995208832223406044287944134381117528642679552487333128252142131375359121793415508211496537834167452784284536412724329064314223405159310288664698829339128075019109308684058290006604682848685269913509631592635991061291910848588896849609801431590823098722376633326783362640890294817117428008250878447534568070706275064939211743430644145531251375121140903646902567865610094036531542192311798685900713287395549217317410664958760801352234360088468504559813801913565013476710034147368107","3210",
1906        "54022537803660305299034619951080701815883244369748020425073564798052876966542443431399087275857901045210724481101594546860129641620399468102807793533371090463155019439376300426001111110082554857406340378848339696012824680258786824575893953186081095194738031315027987901599699063076400850373460929533644671669856383541718682959908309284048987253429549866158169260473899283324734840555471844451471581906629081405048357875350755464695286397907639138344100989311969291336190793916170606327076279174475858929332637620247498179721571404122810497608597392310341741689230018260320576738940580419471886026121","42",
1907        "-178631832601908536171004558032761096296603555642069337022466657948224618220000462860819307447255038425012723042570990636159393162920489265496945881843415223197089039558399046230412722134470564185454526546692168928856644009854678892565378884116563496393911506002791803756925372034984042436244890477408434546164564440123159575118382731762120541349243900708792397436983071602078285813659425066659156365690242457074462226267478374244243395380526750495191218544888861880143020865916520055167383158922400486128739977249726236354738468036532667283720083951543434880992597829033564904575040629182418739065389001","30",
1908        "791452428000395052990463674367408774641299277169449382275797418384963666743647706199968141726577699871535620647877271680371011363641725875016441804003162706216814077280766553803448788712737779201122868425785256450531121303233116813021906603098871081632396261587905577208955217566658493203885877312418398060846091120975165767009344894679683148265158035296420969642632188332444968324025169811975004853721487927435823390089679208323886268818820308870377103259885100640557295098484358112295185600733503211183865432570364543715288518332255545140994972346903247999234036372815320143622542006284391778075614230370915","28446",
1909        "-14750577990707069090892640643627890396790588057597932870361312797605339877218555318290483697927533647647152620207786705312787195868985834595937926758107692208298403338603695989231215167692521223914659443461688651180043473122600527834618929052569980773595567060954830287077779819370475900148020136696922033144977602316267162377565896641961582226753571000379259049438007003714263998786954349001656043288757745179230747108219268158244213583535560775121257993743283335146172958812409641851314398652207671071800026121000251318559834563893238437683153099891646143428578083417759291169439869395337214083816200847999672131739013721","112409792943630",
1910        "3747779487230132491900784839670360702621969953197210049501251358419118141986190115493564102700854690955620245536804192644210548283959203482329576239229460642292110251158496706142540965642461705544532623844725391336993674360586090434899264849404330015359243509766947985382840962057755044491819343032537950612194324352600097426145059000885091707883730246793387324549656578506520410600783165826659781478577755966537904539285939024535631431583555662911156572916830417427298791392192762504813321994018595884981406242206002171583618182559028568019214315784771237935918859841167713230012475726056827372638626381963585567","6",
1911        "-90024360230387232926513776967323326794697521699014742623495978698606015115117275163938036449925194916848722477911818639649193586017040660336424190717814222779195576374240102808595115759449145109628358599211511306846046027111627976699604296466038648254628913968475783753631135765609682725844055075467512860033319632741902401109514097069147699670253945763189080474357381924555174022637039432337932007567967360013189225484557699891500723386807516891972255376800483482791579348941380904127063436404359870716285253826928676282091132256974170075250937201408937099856598602948280585945432968199829410867368026769199252595801","30",
1912        "268255361246337126367778228815259545045868431029437849645589178925638795843497414272945160749514714084165426561127725005493039810079480452517290207498185200717509668055892683546318815665341535697032808106797631781867079778636900617376382054640817598394038676153517033958597909974830544311249490853122270098045783563026899566057862076965478508567587449877041567729076139718636854981066664507897053160317558190734639340443419436918626125596334898307900637096349388220502296543130848656587202710882968136604232104512201075293126007959863944130551986019010434312311477565523149656237135607816533477227160969517648113475415388407","18438",
1913        "-1971549561515498045945281707781358725202052418820469210287100375124408635052889958528421462206148998040363604314916692108315970766317682360194568697651597939992942608607863240436440670587263898825900455009148047595679837461799849465343687138591048044960789240368282898014290879082621858961237901928382366671020037189395193173411370422088058940123720113875828797411348581899556685938975754080936775135639140747303520886284176117348920517204041675362641162399925576576653291783247744987573436710429594592096815313807926352292662762861675065257056115521036077036706726455044401225578559120921795944351508096992367201543994966078858397","27695910",
1914        "934217844394639862121876161182125782629662966798874841047916685310521671553962326103991809647198275117953002131456853984058360472753503867675200497447364861473654065939291320525986944570166180778722315422080724702421323386482784126733369118721815203984881266971497573448614086191554625693617723971683850141789064565402372607797348631247017594179350444633051496340472506374723019936564961282558140233726105870589191017887960191293607258806456016263231512800235871497605920995641449875958591869460832774676734017514779283735936104646380107986104389604327378999649850427118391694149262577932479289635308490619974233553612397882181813","2658",
1915        "-158845135682583637310050414378214595799174539534413035558457927180260655307343938945950776141757261524377755785578389794390713400424204955224064287421714210308705949074744929596932632241358688262813726174958748320025797188191224303535908875129465586234092629175800592249211106320935394054787272761311671502199224520443098879151573888801771623192883526328927129807622792528335792700318705689540933290598347707055522491251314640804320693900281215224262309217387438431268666875059083616022905739177343532424279695920452691897028237708410245697642781733124304062104070106669202262888345577753573178251769219880004565578924500260455594131928657","90709710",
1916        "52820988550912804071451884491873707077611397679084704441786148606552305037050330902476831330763618274617948135285949957728028695870688501651855365293214261665763024004555565085202897776922988749231163862901818153108264790734253251559229193150038194654616836673701295399639842711610882623238385854157477150420529686140635722118527587063478286467506193498289348590140209058848431175257465051090606564081140432276211107220257560552836024246626312132223879207171064190421730933829702220468746118116687590172977326777807226148473634828964585294322560361223896420954749054906526495733099883401725824247665659378641424764700866014787892724263","6",
1917        "-33509969560310951077673300385182939173120218233464878294252193880216174489593173505823838246617056145585882555706860451887313675357842060706217782059762885287365030658978615168275599290138352468354789757259008089648541137547348110439089815672768327313180473534796888815851856313394644427062261422524265399511451721069073409293158382781768172256354047772889044225230507510811482737634673210381142106400974943160605739571810507050337114409751599167906471400307782233953490003352622899793151422802066630790773271959078676470883432223098011999197297967048507289842238199051754189231014498571762573655266437187969662451821328060350075608863875769571341","750400230",
1918        "9391055572765375337084396117314557817718180166181509161271561059129443589397485081513018424705071474206174166639082356181454829444072491877870552062694059378849770958442322223463829270989901750468144609171516219648902784638712130372569713998758894226403154437364546237046936234486878313444740508912905759996491085406313529674400225188471243945850641568659710670402270969430805001894842922440442618558016542494360661268867597331176547554831945263940482855409773467788716571946235493331362953065136196925332392537268524281812660120686916907551456071718059004900562940126515655638458838017179782433551203385736031994621640626329584486625252298313721425","41089818",
1919        "-8036788600329789637554905408606434751151763343827636999239373595458189088012226087443084910869936431803106244357312755090962116368482475276694252905815888538084063223859971526076599239358507163699508542178943276854150562290682856423508900798777112787185623445870747150615150782172422604213843730688721433197193927654307261557625517124764287228841421171512634534585391620401335799765458353101463743603983062786131145017655611914534755704275915462765633857066185411573999789535405190641054970837982446446314245641839258145246102021303507933166787126508119006984137728812640057612823780541560848078401500756860998000273851310662643557278203530509156427","6810",
1920        "36887651095190359395978631707450897701552180104326675120919609694820560648030474979171133539537862210662044736026621541679008373532300695155313150561362479551904410240432070856446263923979299500092444516028968803538104878672539814981236198854450478175391010338000353791751383268673072000938944246515277386132140583257856132696558562847161814996976075174769694748695438632575090586824991284249542651404132205213670261425227222486264304182560679963038720949995579422933796694986011135385677383318542530506241731018476230796808445574699332521366132935744942092063737117719972419446305379041821531028902499051591918182645467840089773846410376436748099927","6",
1921        "-9231251497205337786805280106627035365614237424730160108604924651736976180075798896215248792654225922031690989989936814428616639799004585592248516552105228745115600772090469618649949110864874437364870818222898493364585027928374556652768749981360929058298243205794495656723748164403781748730579086335638752417431277593473211309077070359896017274789066550552469547698979828077256389389589765286385150263946322137435703489322701057413339409715907420401015745559280769997899350755949293570194811595617539775432300877948131084500469766803745905331871839118178605120693482483577762003279580322707699760471279269534859347690610372106577357727353920929624706796074041823","285702690",
1922        "1027825635381261025084137509518564962560728450132223539521526948817353636181702656928583828828717439941582449877622669707429906934288552278291471518913824570952553234368789724603442838567318141441615655257676812187239952554217908948438076178377910225560489116748839886969453337695309781684455928141552611581080216941930192968523416897195063762538557409287539980661256894009392408089358915496977576627225148859236468505369342003480428366798569334698706219583396397757047422466674537196192952267180837549901710490935337309113676890484857602521544773316821275084647927509602731787989854783170081148513873085444058741939360926584684355544622280387129474006982319","6",
1923        "-6550760154113509955322519520324689762818345309330319372044343946911206397745316500483812442237014438941335293944750087907405849513426436624910542114024406950056919447434524096723352448017446297358873710248617421245570015225121262135612027972284197922153718331776532945085169610851088437877134878008919341980631531214431991715959299563428081327762284421994333768251921757188656927704213872171465925501861603885110199062495929817244864046765240300372782070942306808931353298237329193104023648523727135025613272860312485759776097808925561726646874622086120677449738708016595010877049123286575216517173385795826972669505322466175106064791167244415577645190187783109558537","7150110",
1924        "6368886588855890529015573157270243329857473927510055032448029799766753973277733490743802481131249684574480516536707396064341401134831688716527812955468642497995612469570965195817160487421985680788793411096140356661976146745150228102075529102340373723482996983265766606390143860966670113764902836475227042342140853703902690953653738071098617987052339494008437294258743953709135498602090607507693690330811502343804351805921729587694688949530767925472547494085687987805410410505963373499481785590861551304643751957929639863906451583045448988420712969430083368919548554903489673904877295745230080254098114031395040280378645921227902389659674535661931328259201472684713522055433","1288550298",
1925        "-188572988301799749241127705615416162812212213841895020766690762649760470928083270444552850572643956114497215752611047787304109375175847805695151784532184851620320647579855705925496701761448465808354795813973779680808888167818764964700925538842134054390992956001669225520942778629366261309269472085949866705644243058822229069911401050274835800403163105897188918692416516840686738234446171587614205164199212128303162013809506690718962024295099161224082167003915292076143373222944136866167556404170101140178710147326875675233494871699877286279336655831877302931691870727029196060937873178589575033953909097962613832380524097498719417990532304114093172036385627503846689476887699","7010970",
1926        "413664867987222956244209355804553132010737311619873937552617773013514076367023256312176160835062661892330568050774492900833764272232882273128949537336490168744752502048943049241434525025420643944559491826884487174091266835019079457750298308269731643882842579624683236689298288625957015934020852065186333120537819299541510992711232749590370342660947186143661434743740645683801183663840994358198551336011498396150912151423938206307866725487635171065947029616755552973030152102744795495331947296938243287814790516243968203740094323906300807997945263580249967044922095488324059855868304323339162224087731767838889133184638682384777584272685488996638922752221799803151101166568141","2802",
1927        "-1031107841715829258343463911912030698033011526937248355743056689674697602723540112990172956648590063031143260600724004347065027947507100488450593390705595096537384272550134562580809036689346664720973019306107636742232165503223693884139986374204369981544077656811163988149479258369895200746359684803254154401124622678491322587940166253688363451761279210145812935556959005617096357884154904315169217040557989251752098016262897990226072225936970292748900280469462419239068829802849814616525273947160256662458312243507612184481586372693712226732955571148746164223940290646509705984881548422992991644255618458312931517172972824955609103863870050004830358982153156128215328921678827090764016307","1261596819210",
1928        "301188512670574409973522039731799774332284237966284386701245748565406164428681036128926063056095855532173991462477799681429260898037617758622985941710976235577943559100259786932475830383970536995316845412258295246742979426181833225367398491810295210078233611825547244385629725708901995508766557190851612705431746456650510895626334976810141560416672624945285953742157497174172068239724208754377218814505115022247180573164262593883452362474468782982223281207823275802362572430210812017251300349771480681113489830696638136842926883664618708287688430079713258060558978096084995295625323965331459365135957743350303225791768643943305879602716567061428469743627314569410061103860282955485","66",
1929        "-770937004570847542710916650565931789612450235859854061364164699012203056638918055135259895416107273426948797500232897203853071790393885116465171393875976673693385000879337198213556447794398434782352654738563122022624642967864571454383814267648876048449590603480115693689162580755992686116534563663816203506520969500954741353683081691871170841572853932684457019190727452579388269518761070066012225637906766598808889134766906870857170491069727522854987604110140202165701028056317715546755482888281141979064206161119893213557441124682697844058662867837317718568074632672139404775850673801367503450134341628737627508058275331031355246875710092188894581462212753275513865164976897399277601","30",
1930        "6129522921698903924225094098894493975355342495202583967070268649960491527425425546286128994652359747795614576729300901193288468364098351921425898605865990737737763196859623939630382031884908027628137877044490771874628629780057338048896492947577968007519815003168246395643465584343043193237880881635998419368585273823574720837496255282601496997088803126305575254814207239481124819310777601851780008004534530982225452330277362311532174615219685651873091360796309808900938980799843045197724124205459703692404692815206856959614667498156439854860378606039596585529438689777885628862162187008631192421450307071680585728728072572580947780239520470275779631590328798136283209653623721740978652009","42",
1931        "-173794233423007364057178441776478220477329295247992746246194028033098955394598434863688615219916214324301023001093357893746330083903744204991166650812061841423746487587175286790822799578241212757499312100323647899801856816587814694250535160608651214203178715194330829602196319299830498255481024809653516681647938371069642056989453238492194899464534194575984259211778777474726825907276342679607943595029925764974686120493830750858815144111018156639748863065679944594300279690320295481014774768770077734050367942856894507914571711299927762683035011101547081082113463840150439878440505062095907817860452253434600485043521832982039217192780017205043143542812887267054217968944656636429743815938366731","207930",
1932        "13873674832923745170636547493407974778812031787162484396922203584245099363875266410392579249977654901546226138255725682890165166445534326211967788770965050606335620329629392266509006667636060346080828403980900870728472844701200849051031032721110383896992993960600720447264177011379489696922012639234054945046493515830646916075528778995875263891724422324224871104321292779816866047869594198817075501831319907776354543266924566231340617726210269713505336228589002649382821152876769886796237625064939931597087898226387660507299139385252112368241337865982079462586705656758449262979906152409703115811425535504703782935179149721682887228861005266932595149810184020272544138033841220106232249961984874041","2874",
1933        "-26012984231082938122321205711106137025736836308811099420978321017146930273005671218213801512995150233049010223325808387207760652312185884489040204754381130425223206104908949348483694397053700779323806598166397483624967026385155584314813089275200872985931085543928822164664368156111249818285706698582816847860316768251385568355911376789419906820360322920270272727732554866797105317523759568495736021299387762599572352567156747373223964524139683333489545030164669822955224449969772265055374313158826483073058617899550929715104043848380098396583498089878468223019258239295983312553586622706514962679548326500707178484227310027019877548445744372693423428185302230083093969160616600533454208539007576331211435405256579","925269860885370",
1934        "990616118004405209547336756803741296643221987382321931721822820900299097210222317336001831102939236490395767769383589450743538317573906480157831910388399026009549363482047723522295682385225255084608513720106406633464966275709378145660935302980557499077874664887504148114522460257788204486783699727500508138689641080038916313291730518751063292404119140988768624151887967557709255327422581104429283482510069570158085040668138860662778558327704232203358594262804010229103345601462990376149249652299540260177492690219461257167020609270211173637016502365724336529333125021542322561542486553234490872949324335738034611680591470710376104053043734142034650737431499785882884589522913671849942865299633710178111","6",
1935        "-674583921974243806727583626605726893305001263343470277847644141353711283719144748396214643464815527720051176636837830641202715031912558545751272265963459158832678711777574436725109883048688054377892408030247686556398952252471378342025741011438575099484896341336771053946168596117603779433009583486087778261092733102095362575294914815257723299413552987428930625441349976483182063565793384408266465380474529520534826114954047630754878668040735849610788709880820541989250906334800119386014176221640220454426941644202670624599594997416757514254570930031264863092952803531973860879930032899029342578124854325385199562765121526103213588359312368828453688417972133504757523704888418645741605357575301183566273541023","690",
1936        "369763997577304898527792056092078784633064579850215089906718927953975408532819407204558577594984146089620623998022051043813802512840569387856574528038951679028036813020229062675066916989178254799161955591299056486170972788431152654760791766791190599466061280142672740585701636947897534361260451582259758691868004875864277790983654375073056163053627813804545640618169241729192679530964042280652574338325670662894789508966565337142721844903092886592945025900894162911247219665577745221033424450158824779569415157792860757313903143859765449917326253963605938268777389669269349544451932422774330893531316131993509969738757701435869993709233337662014445697606454995935095546520785205839498933104702647973917548450654482159","63346038",
1937        "-1054181687381409638131616910853911354455980657913781692616470855145305376442167447677966359041858772154880528891461759759408670524036867330393458847554077805022113371551512761301353668641881943315640673862411986402084705908271833037248441148559758203373649679290148162009307104085713861437407027770249084665069566424802818695978402723455508336525572007924797995088668996393412803006736971627889522618048536287908252622198057721526239638561752413227961375709262141133204448616567780414625262935283998346756442958214463663418419136319815181399963230777047415998079220041299699187067412714951204651656933821520847527577289722382000754063804433988042815372396860021515945936788059415457004917777955610189853432234056001","30",
1938        "490708060345718168482949551065008296976352341812156850912460754793773228598743021163789697895339848508047563198613297970895771272423084451586762411542433689884882582217920150615320261296412657655420294950882472519077107436411891314989684272827405143453664676299901282427746214002776016429289992734596941647026794542862610180682813869753438462413802302318288921935542654412271101075502309306996459972517715826132089736597584320511669263794701105417099176436055541622722611853665479429929243382275785136760788118477444370781356616739167758192241507716005871662971992278037186350173271291067722673038829083878521474773044108789762491726385320922397469619970046438378452390321367551511627206765616414040686450916033159995734945","2300826",
1939        "-295710682036368660462889828594612232955924955002201524705248520206592222985898208468564409024643552663319217301279279456599515167727625150606450018918397365866826625510003837807762234721907229268939728260753209578951485832358379661753320905096937054954098245328282838820612002927999827875337943367582502020545668956501113411001279555402079055392887648432396357270387066730514118593546561439603919778742559512015901475146476123141564243291327909428505227081149290897384800809797729594081307795911125761034063373560371910341119307145118793762233074626775908914657621875311827273319354160489425009975655550242275976199551186361927525119846700395907422799296606519248369841030425783198567324724905036753320437659605858401058641353","226590",
1940        "48304952056070024985817756788239978908682438355340656791822950062590796355241246689448963611728938505601369463215719576289607035969877364225474402150107267021033472617633002498666739540719726067208908400868629471652429458078768484626439283884479517280978247192859399400878296914805548929586054952054550941690046726387385343252474507896205536559285343999286125673465838911018803676968634234058996143076453319007358347695097769648224143191413651220988411012365741972645136470905240747776225223601701498704103140362659865756058524036705273180849600017203422636276522059566343973347678362590339657974928808786150211385800610012387411013195229675295096785833478392423958039918096194458403331849842918089944679987935690252652536647","6",
1941        "-25535109222241698972694173193973518774716242178265877234279078154218306288190453574391345349607020347701104389306082727640662249328994752154372973422361999092715922759633521270378541507790769763357456276378726595626459724623679603915577643947700472106673593583118436960062894166402312867190376592931480464022802240967436286930003098509086131781553808335933575272888361402715107788113154287996385861486207257531899222760762747086758473894396360217594584387656887231071811145643960111628458232775253525983817907496407600357182919744162523813712195436135703038816616947504736640465073624938220282345566631808140813969037396579578863663701714202683178828730441933582773514024218332915875259194767588507001149870106803558286434937883617","510",
1942        "1098651228481515846116735605942901657544418942763685605483819128454976385464899434737242367240952385067444749008365642568063670579728409673461751272924549125285602795045413203741844230426358616759853983294227411821897092348333446120995756339616373490620446536298890167517709746328742870299157366295356390032576984786436553789735950197578753614711411671045054914305425293539519785975913292188458435472840034103326532984267720997838197057411785852422605414587329459898626238314431604983174665957920693815755227736013545929781603072985638923688324355299338791368472239159091263112230452362400314432752766640687084000308229806355172210181524971644942859893996940410875091940649202106109066788663807627887433701124817799811432100659519242101309","3499986"
1943    );
1944    my $P = shift || $DECIMALS;
1945    my $x = Math::BigFloat->new( $bernoulli[ $n - 2 ] );
1946    my $y = Math::BigFloat->new( $bernoulli[ $n - 1 ] );
1947    return wantarray ? ( $x, $y ) : $x->bdiv( $y )->bfround( -$P );
1948}
1949
19501;
1951__END__
1952
1953=head1 NAME
1954
1955Math::NumberCruncher - Collection of useful math-related functions.
1956
1957=head1 AUTHOR
1958
1959Kurt Kincaid <sifukurt@yahoo.com>
1960