Tue Apr 15 00:34:45 2008 run on win32
%
% Testing file for REDUCE Special Functions Package
%
% Chris Cannam, ZIB Berlin
% October 1992 -> Feb 1993
% (only some of the time, of course)
%
% Corrections and comments to neun@sc.zib-berlin.de
%
on savesfs;
% just in case it's off for some reason
off bfspace;
% to provide more similarity between runs
% with old & new bigfloats
let {sinh (~x) => (exp(x) - exp (-x))/2 };
% this will improve some results
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
% 1. Bernoulli numbers
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
off rounded;
procedure do!*one!*bern(x);
write ("Bernoulli ", x, " is ", bernoulli x);
do*one*bern
do!*one!*bern(1);
- 1
Bernoulli 1 is ------
2
do!*one!*bern(2);
1
Bernoulli 2 is ---
6
do!*one!*bern(3);
Bernoulli 3 is 0
do!*one!*bern(13);
Bernoulli 13 is 0
do!*one!*bern(14);
7
Bernoulli 14 is ---
6
do!*one!*bern(300);
Bernoulli 300 is ( - 186387899520485901199504534184815606618219184663590593751\
8715320655775958174360523134990756922303410810482600528769\
4796420210012184158790061643029553704608291464348079647177\
3719535693514415158342483315425004774743357558499902912677\
5186293388721514970183351129809976971603227633930434923843\
9848295803115933725653985747628800282891676355700124156069\
41367995702212211519561707046505473575241)/866054419230
do!*one!*bern(-2);
Bernoulli -2 is bernoulli(-2)
do!*one!*bern(0);
Bernoulli 0 is 1
for n := 2 step 2 until 100 do do!*one!*bern n;
1
Bernoulli 2 is ---
6
- 1
Bernoulli 4 is ------
30
1
Bernoulli 6 is ----
42
- 1
Bernoulli 8 is ------
30
5
Bernoulli 10 is ----
66
- 691
Bernoulli 12 is --------
2730
7
Bernoulli 14 is ---
6
- 3617
Bernoulli 16 is ---------
510
43867
Bernoulli 18 is -------
798
- 174611
Bernoulli 20 is -----------
330
854513
Bernoulli 22 is --------
138
- 236364091
Bernoulli 24 is --------------
2730
8553103
Bernoulli 26 is ---------
6
- 23749461029
Bernoulli 28 is ----------------
870
8615841276005
Bernoulli 30 is ---------------
14322
- 7709321041217
Bernoulli 32 is ------------------
510
2577687858367
Bernoulli 34 is ---------------
6
- 26315271553053477373
Bernoulli 36 is -------------------------
1919190
2929993913841559
Bernoulli 38 is ------------------
6
- 261082718496449122051
Bernoulli 40 is --------------------------
13530
1520097643918070802691
Bernoulli 42 is ------------------------
1806
- 27833269579301024235023
Bernoulli 44 is ----------------------------
690
596451111593912163277961
Bernoulli 46 is --------------------------
282
- 5609403368997817686249127547
Bernoulli 48 is ---------------------------------
46410
495057205241079648212477525
Bernoulli 50 is -----------------------------
66
- 801165718135489957347924991853
Bernoulli 52 is -----------------------------------
1590
29149963634884862421418123812691
Bernoulli 54 is ----------------------------------
798
- 2479392929313226753685415739663229
Bernoulli 56 is ---------------------------------------
870
84483613348880041862046775994036021
Bernoulli 58 is -------------------------------------
354
- 1215233140483755572040304994079820246041491
Bernoulli 60 is ------------------------------------------------
56786730
12300585434086858541953039857403386151
Bernoulli 62 is ----------------------------------------
6
- 106783830147866529886385444979142647942017
Bernoulli 64 is -----------------------------------------------
510
1472600022126335654051619428551932342241899101
Bernoulli 66 is ------------------------------------------------
64722
- 78773130858718728141909149208474606244347001
Bernoulli 68 is -------------------------------------------------
30
1505381347333367003803076567377857208511438160235
Bernoulli 70 is ---------------------------------------------------
4686
- 5827954961669944110438277244641067365282488301844260429
Bernoulli 72 is ------------------------------------------------------------
140100870
34152417289221168014330073731472635186688307783087
Bernoulli 74 is ----------------------------------------------------
6
- 24655088825935372707687196040585199904365267828865801
Bernoulli 76 is ----------------------------------------------------------
30
414846365575400828295179035549542073492199375372400483487
Bernoulli 78 is -----------------------------------------------------------
3318
Bernoulli 80 is ( - 4603784299479457646935574969019046849794257872751288919656\
867)/230010
1677014149185145836823154509786269900207736027570253414881613
Bernoulli 82 is ---------------------------------------------------------------
498
Bernoulli 84 is ( - 2024576195935290360231131160111731009989917391198090877281\
083932477)/3404310
Bernoulli 86 is 660714619417678653573847847426261496277830686653388931761996983/
6
Bernoulli 88 is ( - 1311426488674017507995511424019311843345750275572028644296\
919890574047)/61410
Bernoulli 90 is 11790572790210827998841233512492150837752549496696471162315452\
15727922535/272118
Bernoulli 92 is ( - 1295585948207537527989427828538576749659341483719435143023\
316326829946247)/1410
Bernoulli 94 is 12208138065797444696073016794132012039585084152026966214362151\
05284649447/6
Bernoulli 96 is ( - 2116004495972665130975977281098242336730439543890602341506\
38733420050668349987259)/4501770
Bernoulli 98 is 67908260672905495624051117546403605607342195728504487509073961\
249992947058239/6
Bernoulli 100 is ( - 945980378191221252952274330694937218727028415330669361333\
85696204311395415197247711)/33330
on rounded;
precision 100;
12
do!*one!*bern(1);
Bernoulli 1 is - 0.5
do!*one!*bern(2);
Bernoulli 2 is
0.1666666666666666666666666666666666666666666666666666666666666\
666666666666666666666666666666666666667
do!*one!*bern(3);
Bernoulli 3 is 0
do!*one!*bern(13);
Bernoulli 13 is 0
do!*one!*bern(14);
Bernoulli 14 is
1.166666666666666666666666666666666666666666666666666666666666\
666666666666666666666666666666666666667
do!*one!*bern(300);
Bernoulli 300 is -
2.15214997327998682971981737675608819857345657278380670184934\
3925846796995824343474766557085976828531e+375
do!*one!*bern(-2);
Bernoulli -2 is bernoulli(-2)
do!*one!*bern(0);
Bernoulli 0 is 1
do!*one!*bern(38);
Bernoulli 38 is
4.883323189735931666666666666666666666666666666666666666666666\
666666666666666666666666666666666666667e+14
do!*one!*bern(400);
Bernoulli 400 is -
6.84694485580645336061625858231088359767823009718062574241478\
1150311357197834589950328990573681886442e+549
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
% 2. Gamma function
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
on rounded;
off complex;
precision 40;
100
algebraic procedure wg(x);
write ("gamma (", x, ") ==> ", gamma x);
wg
algebraic procedure wp(x);
write ("-- precision ", x, ", from ", precision(x));
wp
wg (1/2);
gamma (0.5) ==> 1.772453850905516027298167483341145182798
wg (3/2);
gamma (1.5) ==> 0.8862269254527580136490837416705725913988
write ("sqrt(pi)/2 ==> ", sqrt(pi)/2);
sqrt(pi)/2 ==> 0.8862269254527580136490837416705725913988
wp(10);
-- precision 10, from 40
for x := 0 step 5 until 100 do
<< wg (1 + x/1000);
wg (-1 - x/13);
wp (8+floor(x/4)) >>;
gamma (1) ==> 1
gamma (-1) ==> gamma(-1)
-- precision 8, from 10
gamma (1.005) ==> 0.99713854
gamma ( - 1.3846154) ==> 2.7320314
-- precision 9, from 8
gamma (1.01) ==> 0.994325851
gamma ( - 1.76923077) ==> 2.89933597
-- precision 10, from 9
gamma (1.015) ==> 0.9915612888
gamma ( - 2.153846154) ==> - 2.919307224
-- precision 11, from 10
gamma (1.02) ==> 0.98884420326
gamma ( - 2.5384615385) ==> - 0.91247160689
-- precision 13, from 11
gamma (1.025) ==> 0.9861739631483
gamma ( - 2.923076923077) ==> - 2.407817725014
-- precision 14, from 13
gamma (1.03) ==> 0.98354995055382
gamma ( - 3.3076923076923) ==> 0.42665848359037
-- precision 15, from 14
gamma (1.035) ==> 0.980971560550586
gamma ( - 3.69230769230769) ==> 0.250121998146955
-- precision 16, from 15
gamma (1.04) ==> 0.9784382009142447
gamma ( - 4.076923076923077) ==> - 0.4868211588210416
-- precision 18, from 16
gamma (1.045) ==> 0.975949291822951489
gamma ( - 4.46153846153846154) ==> - 0.064315864992688343
-- precision 19, from 18
gamma (1.05) ==> 0.9735042655627756432
gamma ( - 4.846153846153846154) ==> - 0.07308480130893001268
-- precision 20, from 19
gamma (1.055) ==> 0.9711025662416699039
gamma ( - 5.2307692307692307692) ==> 0.026504298643014994546
-- precision 21, from 20
gamma (1.06) ==> 0.968743649511638364209
gamma ( - 5.61538461538461538462) ==> 0.00947958151841097406813
-- precision 23, from 21
gamma (1.065) ==> 0.96642698229883993296188
gamma (-6) ==> gamma(-6)
-- precision 24, from 23
gamma (1.07) ==> 0.964152042541366448869499
gamma ( - 6.38461538461538461538462) ==> - 0.00224562428754660955672785
-- precision 25, from 24
gamma (1.075) ==> 0.9619183189344448686422338
gamma ( - 6.769230769230769230769231) ==> - 0.0014913552210721080976799
-- precision 26, from 25
gamma (1.08) ==> 0.95972531068282223532653464
gamma ( - 7.1538461538461538461538462) ==> 0.0009821350464726696000373901
-- precision 28, from 26
gamma (1.085) ==> 0.9575725272601010524841387298
gamma ( - 7.538461538461538461538461538) ==> 0.0002081367552774423831456972591
-- precision 29, from 28
gamma (1.09) ==> 0.95545948817480124076971838245
gamma ( - 7.9230769230769230769230769231) ==>
0.00038372470359721945415080451977
-- precision 30, from 29
gamma (1.095) ==> 0.95338572274293305256571929565
gamma ( - 8.30769230769230769230769230769) ==>
- 0.0000487302698231202992921248405302
-- precision 31, from 30
gamma (1.1) ==> 0.9513507698668731836292487177265
gamma ( - 8.692307692307692307692307692308) ==>
- 0.00002092711689267207910068115762837
-- precision 33, from 31
wg(1/1000000003);
gamma (0.000000000999999997000000008999999973) ==>
1.00000000242278433608752313084681e+9
off rounded;
gamma(17/2);
2027025*sqrt(pi)
------------------
256
gamma(-17/2);
- 512*pi
-------------------
34459425*sqrt(pi)
gamma(4);
6
gamma(0);
gamma(0)
gamma(-4);
gamma(-4)
gamma(-17/3);
6*pi
------------------------
17
17*sqrt(3)*gamma(----)
3
p := gamma(x**2) * gamma(x-y**gamma(y)) - (1/(gamma(4*(x-y))));
2 gamma(y)
gamma(x )*gamma( - y + x)*gamma(4*x - 4*y) - 1
p := --------------------------------------------------------
gamma(4*x - 4*y)
y := 1/4;
1
y := ---
4
p;
2 1 gamma(1/4)
gamma(x )*gamma( - (---) + x)*gamma(4*x - 1) - 1
4
------------------------------------------------------------
gamma(4*x - 1)
x := 3;
x := 3
p;
1 gamma(1/4)
146313216000*gamma( - (---) + 3) - 1
4
------------------------------------------------
3628800
y := -3/8;
- 3
y := ------
8
p;
(128*(2490343877896875*sqrt(pi)
- 3 (8*pi)/(3*gamma(3/8)*sin((3*pi)/8))
3*(------) - 1
8
*gamma(---------------------------------------------------) - 64))/(
- 3 (8*pi)/(3*gamma(3/8)*sin((3*pi)/8))
(------)
8
7905853580625*sqrt(pi))
on rounded, complex;
*** Domain mode rounded changed to complex-rounded
precision 50;
33
p;
- 0.00000000058461000084165968732153392127582134179078414159599
+ 3.7721251013859508830301986850709684723938237902095e-60*i
off rounded, complex;
*** Domain mode complex-rounded changed to complex
clear y;
p;
gamma(y)
40320*gamma( - y + 3)*gamma( - 4*y + 12) - 1
------------------------------------------------------
gamma( - 4*y + 12)
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
% 3. Beta function. Not very interesting
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
algebraic procedure do!*one!*beta(x,y);
write ("Beta (", x, ",", y, ") = ", beta(x,y));
do*one*beta
do!*one!*beta(0,1);
Beta (0,1) = beta(0,1)
do!*one!*beta(2,-3);
Beta (2,-3) = beta(2,-3)
do!*one!*beta(3,2);
1
Beta (3,2) = ----
12
do!*one!*beta(a+b,(c+d)**(b-a));
b
(c + d)
gamma(a + b)*gamma(----------)
b a
(c + d) (c + d)
Beta (a + b,----------) = ---------------------------------------------
a a a b
(c + d) (c + d) *a + (c + d) *b + (c + d)
gamma(------------------------------------)
a
(c + d)
do!*one!*beta(-3,4);
Beta (-3,4) = beta(-3,4)
do!*one!*beta(-3,2);
Beta (-3,2) = beta(-3,2)
do!*one!*beta(-3,-7.5);
- 15 - 15
Beta (-3,-------) = beta(-3,-------)
2 2
do!*one!*beta((pi * 10), exp(5));
5
5 gamma(e )*gamma(10*pi)
Beta (10*pi,e ) = ------------------------
5
gamma(e + 10*pi)
on rounded;
precision 30;
50
do!*one!*beta(0,1);
Beta (0,1) = beta(0,1)
do!*one!*beta(2,-3);
Beta (2,-3) = beta(2,-3)
do!*one!*beta(3,2);
Beta (3,2) = 0.0833333333333333333333333333333
do!*one!*beta(a+b,(c+d)**(b-a));
b
(c + d)
gamma(a + b)*gamma(----------)
b a
(c + d) (c + d)
Beta (a + b,----------) = ---------------------------------------------
a a a b
(c + d) (c + d) *a + (c + d) *b + (c + d)
gamma(------------------------------------)
a
(c + d)
do!*one!*beta(-3,4);
Beta (-3,4) = beta(-3,4)
do!*one!*beta(-3,2);
Beta (-3,2) = beta(-3,2)
do!*one!*beta(-3,-7.5);
Beta (-3, - 7.5) = beta(-3, - 7.5)
do!*one!*beta((pi * 10), exp(5));
Beta (31.4159265358979323846264338328,148.413159102576603421115580041) =
3.26162024071771351768890966259e-37
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
% 4. Pochhammer notation
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
off rounded;
pochhammer(4,5);
6720
pochhammer(-4,5);
- 24
---------------
factorial(-1)
pochhammer(4,-5);
factorial(-2)
---------------
6
pochhammer(-4,-5);
- 1
-------
15120
pochhammer(17/2,12);
157783444591397625
--------------------
4096
pochhammer(-17/2,12);
- 516891375
--------------
4096
pochhammer(1/3,14)*pochhammer(2/3,15);
148260333813553014031851192320000000
--------------------------------------
68630377364883
q := pochhammer(1/5,11)*pochhammer(2/5,11)*pochhammer(3/5,11)*
pochhammer(1-1/5,11)*pochhammer(1,11)*pochhammer(6/5,11)*
pochhammer(70/50,11)*pochhammer(8/5,11)*pochhammer(9/5,11);
q := 3083594385972516709231715156294765932567301688416624945842705314534752626\
0421323040989302952930905614937458388095527696384131072/12924697071141057\
41986576081359316958696581423282623291015625
on complex;
pochhammer(a+b*i,c)*pochhammer(a-b*i,c);
pochhammer(a - i*b,c)*pochhammer(a + i*b,c)
a := 2;
a := 2
b := 3;
b := 3
c := 5;
c := 5
pochhammer(a+b*i,c)*pochhammer(a-b*i,c);
8950500
off complex;
on rounded;
pochhammer(1/5,11)*pochhammer(2/5,11)*pochhammer(3/5,11)*
pochhammer(1-1/5,11)*pochhammer(1,11)*pochhammer(6/5,11)*
pochhammer(70/50,11)*pochhammer(8/5,11)*pochhammer(9/5,11);
2.38581559706937212593793381562e+67
q;
2.38581559706937212593793381562e+67
pochhammer(pi,floor (pi**8));
2.47253079057195612235973919163e+33625
pochhammer(-pi,floor (pi**7));
5.91750008101140123889628058136e+9185
pochhammer(1.5,floor (pi**8));
1.88808885937650373473836451368e+33619
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
% 5. Digamma function
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
procedure do!*one!*psi(x);
<< precision (precision(0) + 4)$
write("Psi of ", x, " is ", psi(x) ) >> ;
do*one*psi
clear x, y;
z := x * ((x+y)**2 + (x**y));
y 2 2
z := x*(x + x + 2*x*y + y )
off rounded;
do!*one!*psi(3);
1
2*log(2) + psi(---) + psi(1) + 3
2
Psi of 3 is ----------------------------------
2
do!*one!*psi(pi);
Psi of pi is psi(pi)
do!*one!*psi(1.005);
201 1
Psi of ----- is psi(-----) + 200
200 200
do!*one!*psi(1.995);
pi 1
199*cot(-----)*pi + 199*psi(-----) + 200
399 200 200
Psi of ----- is ------------------------------------------
200 199
do!*one!*psi(74);
Psi of 74 is (269426837164032294756360054754050*log(2)
1
+ 134713418582016147378180027377025*psi(---)
2
+ 2138308231460573767907619482175*psi(1)
+ 667084944417653637854891458725877)/13685172681347672114608764\
6859200
do!*one!*psi(-1/2);
- 1 1
Psi of ------ is psi(---) + 2
2 2
do!*one!*psi(-3);
Psi of -3 is infinity
do!*one!*psi(z);
y 2 2 y 3 2 2
Psi of x*(x + x + 2*x*y + y ) is psi(x *x + x + 2*x *y + x*y )
on rounded;
precision 100;
62
do!*one!*psi(3);
Psi of 3 is
0.9227843350984671393934879099175975689578406640600764011942327651\
1513227322233532906305293670825325048537
do!*one!*psi(pi);
Psi of
3.141592653589793238462643383279502884197169399375105820974944592307816\
40628620899862803482534211706798214809 is
0.9772133079420067332920694864061823436408346099943256380095232865318105924777\
14131730207565436292873435576949
do!*one!*psi(1.005);
Psi of 1.005 is -
0.569020911344382831688295690950028809267968870263469094699243\
7771372094405771031990077045450115105840960011306543
do!*one!*psi(1.995);
Psi of 1.995 is
0.419554603028108628651620272287243954692129976503405961275768\
00237211077024989654159358528162438486780461247996589393
do!*one!*psi(74);
Psi of 74 is
4.297293118804667639106350362560367137919643935000983615964648570\
79142072638532815275329153012799334519059885201047316254
do!*one!*psi(-1/2);
Psi of - 0.5 is
0.03648997397857652055902366700124443280684039533956589295287\
2746128345029282945897851326282715415875401365590709051546051\
66846
do!*one!*psi(-3);
Psi of -3 is infinity
do!*one!*psi(z);
y 2 2 y 3 2 2
Psi of x*(x + x + 2*x*y + y ) is psi(x *x + x + 2*x *y + x*y )
precision 15;
132
x := 8/3;
x := 2.66666666666667
y := 7/1000;
y := 0.007
do!*one!*psi(z);
Psi of 21.74768766103287773 is 3.056340330052438423
off rounded;
clear x, y;
df(psi(z), x);
y 3 2 2 y y 2 2
polygamma(1,x *x + x + 2*x *y + x*y )*(x *y + x + 3*x + 4*x*y + y )
df(df(psi(z), y),x);
2*y y 3 2 2
x *log(x)*polygamma(2,x *x + x + 2*x *y + x*y )*x*y
2*y y 3 2 2
+ x *log(x)*polygamma(2,x *x + x + 2*x *y + x*y )*x
y y 3 2 2 3
+ 3*x *log(x)*polygamma(2,x *x + x + 2*x *y + x*y )*x
y y 3 2 2 2
+ 4*x *log(x)*polygamma(2,x *x + x + 2*x *y + x*y )*x *y
y y 3 2 2 2
+ x *log(x)*polygamma(2,x *x + x + 2*x *y + x*y )*x*y
y y 3 2 2
+ x *log(x)*polygamma(1,x *x + x + 2*x *y + x*y )*y
y y 3 2 2
+ x *log(x)*polygamma(1,x *x + x + 2*x *y + x*y )
y y 3 2 2 2
+ 2*x *polygamma(2,x *x + x + 2*x *y + x*y )*x *y
y y 3 2 2 2
+ 2*x *polygamma(2,x *x + x + 2*x *y + x*y )*x
y y 3 2 2 2
+ 2*x *polygamma(2,x *x + x + 2*x *y + x*y )*x*y
y y 3 2 2
+ 2*x *polygamma(2,x *x + x + 2*x *y + x*y )*x*y
y y 3 2 2
+ x *polygamma(1,x *x + x + 2*x *y + x*y )
y 3 2 2 4
+ 6*polygamma(2,x *x + x + 2*x *y + x*y )*x
y 3 2 2 3
+ 14*polygamma(2,x *x + x + 2*x *y + x*y )*x *y
y 3 2 2 2 2
+ 10*polygamma(2,x *x + x + 2*x *y + x*y )*x *y
y 3 2 2 3
+ 2*polygamma(2,x *x + x + 2*x *y + x*y )*x*y
y 3 2 2
+ 4*polygamma(1,x *x + x + 2*x *y + x*y )*x
y 3 2 2
+ 2*polygamma(1,x *x + x + 2*x *y + x*y )*y
int(psi(z), z);
y 3 2 2
log(gamma(x *x + x + 2*x *y + x*y ))
on rounded;
for k := 1 step 0.1 until 2 do do!*one!*psi(k);
Psi of 1 is - 0.57721566490153286060651
Psi of 1.09999999999999999999999999 is - 0.423754940411076795168216226
Psi of 1.19999999999999999999999999483 is - 0.2890398965921882955472079690017
Psi of 1.2999999999999999999999999948304368 is
- 0.16919088886679965563116117990224309
Psi of 1.39999999999999999999999999483043679294 is
- 0.0613845445851161457306754873481741603465
Psi of 1.499999999999999999999999994830436792938162 is
0.03648997397857652055902366216872537099062413
Psi of 1.5999999999999999999999999948304367929381615219 is
0.12604745277347625190600271342271679469086378102
Psi of 1.69999999999999999999999999483043679293816152191799 is
0.20854787487349395667996417990930337500109911872963
Psi of 1.799999999999999999999999994830436792938161521917991479 is
0.2849914332938615406087023630333301945809176983421812009
Psi of 1.8999999999999999999999999948304367929381615219179914789724 is
0.35618416116405971922472708037210055174041765918986424485621
Psi of 1.99999999999999999999999999483043679293816152191799147897240477 is
0.422784335098467139393487906583570145998489083710421411058968004
off rounded;
% PSI_SIMP.TST F.J.Wright, 2 July 1993
on evallhseqp;
factor psi;
on rat, intstr, div;
% for neater output
% Do not try using "off mcd"!
psi(x+m) - psi(x+m-1) = 1/(x+m-1);
1 1
-----------=-----------
m + x - 1 m + x - 1
psi(x+2) - psi(x+1) + 2*psi(x) = 1/(x+1) + 2*psi(x);
1 1
2*psi(x) + -------=2*psi(x) + -------
x + 1 x + 1
psi(x+2) + 3*psi(x) = 4*psi(x) + 1/x + 1/(x+1);
-1 -1
x + 2 x + 2
4*psi(x) + ---------=4*psi(x) + ---------
x + 1 x + 1
psi(x + 1) = psi(x) + 1/x;
-1 -1
psi(x) + x =psi(x) + x
psi(x + 3/2) = psi(x + 1/2) + 1/(x + 1/2);
1 2 1 2
psi(x + ---) + ---------=psi(x + ---) + ---------
2 2*x + 1 2 2*x + 1
psi(x - 1/2) = psi(x + 1/2) - 1/(x - 1/2);
1 2 1 2
psi(x + ---) - ---------=psi(x + ---) - ---------
2 2*x - 1 2 2*x - 1
psi((x + 3a)/a);
-1
1 2*(3*x + 8*x + 12)
psi(---*x) + ----------------------
2 2
x + 6*x + 8
psi(x/y + 3);
-1 2
-1 y*(3*x + 2*x *y + 6*y)
psi(x*y ) + --------------------------
2 2
x + 3*x*y + 2*y
off rat, intstr, div;
on rational;
psi(x+m) - psi(x+m-1) = 1/(x+m-1);
1 1
-----------=-----------
m + x - 1 m + x - 1
psi(x+2) - psi(x+1) + 2*psi(x) = 1/(x+1) + 2*psi(x);
2*psi(x)*x + 2*psi(x) + 1 2*psi(x)*x + 2*psi(x) + 1
---------------------------=---------------------------
x + 1 x + 1
psi(x+2) + 3*psi(x) = 4*psi(x) + 1/x + 1/(x+1);
2 2
4*psi(x)*x + 4*psi(x)*x + 2*x + 1 4*psi(x)*x + 4*psi(x)*x + 2*x + 1
------------------------------------=------------------------------------
2 2
x + x x + x
psi(x + 1) = psi(x) + 1/x;
psi(x)*x + 1 psi(x)*x + 1
--------------=--------------
x x
psi(x + 3/2) = psi(x + 1/2) + 1/(x + 1/2);
1 1 1 1 1 1
psi(x + ---)*x + ---*psi(x + ---) + 1 psi(x + ---)*x + ---*psi(x + ---) + 1
2 2 2 2 2 2
---------------------------------------=---------------------------------------
1 1
x + --- x + ---
2 2
psi(x - 1/2) = psi(x + 1/2) - 1/(x - 1/2);
1 1 1 1 1 1
psi(x + ---)*x - ---*psi(x + ---) - 1 psi(x + ---)*x - ---*psi(x + ---) - 1
2 2 2 2 2 2
---------------------------------------=---------------------------------------
1 1
x - --- x - ---
2 2
psi((x + 3a)/a);
1 2 2 8
psi(---*x)*x*(x + 6*x + 8) + 6*(x + 4*x + ---)
2 3
--------------------------------------------------
2
x*(x + 6*x + 8)
psi(x/y + 3);
x 2 2 2 2 2
psi(---)*x*(x + 3*x*y + 2*y ) + 3*y*(x + 2*x*y + ---*y )
y 3
------------------------------------------------------------
2 2
x*(x + 3*x*y + 2*y )
off rational;
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
% 6. Polygamma functions
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
procedure do!*one!*pg(n,x);
write ("Polygamma (", n, ") of ", x, " is ", polygamma(n,x));
do*one*pg
off rounded;
do!*one!*pg(1,1/2);
1 1
Polygamma (1) of --- is polygamma(1,---)
2 2
do!*one!*pg(1,1);
2
pi
Polygamma (1) of 1 is -----
6
do!*one!*pg(1,3/2);
2
3 pi - 8
Polygamma (1) of --- is ---------
2 2
do!*one!*pg(1,1.005);
201 201
Polygamma (1) of ----- is polygamma(1,-----)
200 200
do!*one!*pg(1,1.995);
399 399
Polygamma (1) of ----- is polygamma(1,-----)
200 200
do!*one!*pg(1,1e-10);
1 1
Polygamma (1) of ------------- is polygamma(1,-------------)
10000000000 10000000000
do!*one!*pg(2,1.45);
29 29
Polygamma (2) of ---- is polygamma(2,----)
20 20
do!*one!*pg(3,1.99);
199 199
Polygamma (3) of ----- is polygamma(3,-----)
100 100
do!*one!*pg(4,-8.2);
- 41 - 41
Polygamma (4) of ------- is polygamma(4,-------)
5 5
do!*one!*pg(5,0);
Polygamma (5) of 0 is infinity
do!*one!*pg(6,-5);
Polygamma (6) of -5 is infinity
do!*one!*pg(7,200);
Polygamma (7) of 200 is (17726903568079021516229662669344581999513453462651624\
712682433450930396569413638060869523209196416456149981841411350319061723276\
968944831274405326841106087088051862580482604049179034792880997036306450715\
491134920388121743532178874504889199026667348488846246287054278592284199023\
548572837863642592371308428279823642858774260678182876022899303888558378920\
013392172441409952343098219076353830605697314645020537255667718925997023946\
826863459069360054494011515187931095594830347819766125471751832158677864778\
297413490571816668099535649328289203181141151677544880676008798702466775519\
536091686976761627732998715124407667063278175501392084867643186474774933465\
8
396962161891279302009565617687762042880000000000000000000000*pi - 16820227\
432460991991188660317955383248583163197077694391321472430675801147397894113\
241309474829652161707416711766105349898340694530697667469860241947942732920\
236500114448913470876512035623237448423938602084906142995199152247445505854\
809437866123469569241946248640222536901257402027830808909074139755120553020\
556624279638002199727273508802505439369471450464398635245844612858514905584\
191236791869017030402868269742260341406479951357060228410068248515181595142\
787674466216948716753957289698597496283251496652438121703517099395342367021\
250291313681953557646793694608947926057983579614775135882107369446401248554\
435651401852361884349691269283596626501141334161214962487130496595028983095\
5139972163546512706990894178649281)/332379441901481653429306175050210912490\
877252424717963362795627204944935676505713641303560172432808552812159526462818\
482407311443167715586395099878270739132900972423384048825922106902366518694430\
745950915458779757277282691228353896966672481750012784165867117882267723605328\
731691535740709943298606962033030246693303602017387715928925429361947910469604\
750251103233276436606433091607681634323856824649594135073543769729862444199003\
003689857550501021762715909773708042403069021620614852595346852975209964593076\
502948221562526866293424905422559646396593953966512675164975671252040991301719\
130814280519993725908582643757436465790651101591268309746402030002476193040535\
46148691267935533164553830400000000000000000000000
on rounded;
precision 100;
63
do!*one!*pg(1,1/2);
Polygamma (1) of 0.5 is
4.9348022005446793094172454999380755676568497036203953132066746881100224112096\
02621500886701859276116
do!*one!*pg(1,1);
Polygamma (1) of 1 is
1.6449340668482264364724151666460251892189499012067984377355582293700074704032\
00873833628900619758705
do!*one!*pg(1,3/2);
Polygamma (1) of 1.5 is
0.9348022005446793094172454999380755676568497036203953132066746881100224112096\
026215008867018592761159
do!*one!*pg(1,1.005);
Polygamma (1) of 1.005 is
1.6329941567556809752535869786110269899606654266292810174237683064726679207401\
36690588658095268582677
do!*one!*pg(1,1.995);
Polygamma (1) of 1.995 is
0.6469608286405823512476399271078276952436864784225077824256922477813603546774\
633500406650324443323374
do!*one!*pg(1,1e-10);
Polygamma (1) of 0.0000000001 is
1.0000000000000000000164493406660781505587296600657326399164207520577564247801\
80942170336642993019363e+20
do!*one!*pg(2,1.45);
Polygamma (2) of 1.45 is -
0.9038374030762576882323095385298671928600903081114466582453675225591714826271\
66928803699348863802546
do!*one!*pg(3,1.99);
Polygamma (3) of 1.99 is
0.5029071324168653201805109908625701971861516895613080444480934082942919777908\
933506993776280887706183
do!*one!*pg(4,-8.2);
Polygamma (4) of - 8.2 is
74935.512595774270527120292307813135522959054934662779176176297186107113915939\
37121099755096024434196
do!*one!*pg(5,0);
Polygamma (5) of 0 is infinity
do!*one!*pg(6,-5);
Polygamma (6) of -5 is infinity
do!*one!*pg(7,200);
Polygamma (7) of 200 is
5.7240937253925583738370296815447272453468073681959561547435678473859266441297\
53317064645661270233713e-14
off rounded;
clear x;
% Polygamma differentiation has already
% been tried a bit in the psi section
df(int(int(int(polygamma(3,x),x),x),x),x);
polygamma(1,x)
clear w, y, z;
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
% 7. Zeta function
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
procedure do!*one!*zeta(n);
write ("Zeta of ", n, " is ", zeta n);
do*one*zeta
off rounded;
clear x, y, z;
z := x * ((x+y)**5 + (x**y));
y 5 4 3 2 2 3 4 5
z := x*(x + x + 5*x *y + 10*x *y + 10*x *y + 5*x*y + y )
do!*one!*zeta(0);
- 1
Zeta of 0 is ------
2
for k := 4 step 2 until 35 do
do!*one!*zeta(k);
4
pi
Zeta of 4 is -----
90
6
pi
Zeta of 6 is -----
945
8
pi
Zeta of 8 is ------
9450
10
pi
Zeta of 10 is -------
93555
12
691*pi
Zeta of 12 is -----------
638512875
14
2*pi
Zeta of 14 is ----------
18243225
16
3617*pi
Zeta of 16 is --------------
325641566250
18
43867*pi
Zeta of 18 is ----------------
38979295480125
20
174611*pi
Zeta of 20 is ------------------
1531329465290625
22
155366*pi
Zeta of 22 is -------------------
13447856940643125
24
236364091*pi
Zeta of 24 is -----------------------
201919571963756521875
26
1315862*pi
Zeta of 26 is ----------------------
11094481976030578125
28
6785560294*pi
Zeta of 28 is --------------------------
564653660170076273671875
30
6892673020804*pi
Zeta of 30 is ------------------------------
5660878804669082674070015625
Zeta of 32 is zeta(32)
Zeta of 34 is zeta(34)
do!*one!*zeta(-17/3);
- 17 - 17
Zeta of ------- is zeta(-------)
3 3
do!*one!*zeta(190);
Zeta of 190 is zeta(190)
do!*one!*zeta(300);
Zeta of 300 is zeta(300)
do!*one!*zeta(0);
- 1
Zeta of 0 is ------
2
do!*one!*zeta(-44);
Zeta of -44 is 0
on rounded;
clear x, y;
for k := 3 step 3 until 36 do <<
precision (31+k*3);
do!*one!*zeta(k) >>;
Zeta of 3 is 1.202056903159594285399738161511449990765
Zeta of 6 is 1.017343061984449139714517929790920527901817490033
Zeta of 9 is 1.002008392826082214417852769232412060485605851394888756549
Zeta of 12 is
1.00024608655330804829863799804773967096041608845800340453304095\
2133
Zeta of 15 is
1.00003058823630702049355172851064506258762794870685817750656993\
2893332267156
Zeta of 18 is
1.00000381729326499983985646164462193973045469721895333114317442\
9987630039542650045638
Zeta of 21 is
1.00000047693298678780646311671960437304596644669478493760020748\
7376596839087898159833876638564
Zeta of 24 is
1.00000005960818905125947961244020793580122750391883730279586424\
6972321724495355468544848206832825003614
Zeta of 27 is
1.00000000745071178983542949198100417060411945471903188256582999\
3239578352147606271570867900837100313523764933952
Zeta of 30 is
1.00000000093132743241966818287176473502121981356795513681618500\
8613360441960672940496363503624604027929086312123388047291
Zeta of 33 is
1.00000000011641550172700519775929738354563095165224717276359325\
6517739947029124624567548673934974376008810870912845774213829513\
369
Zeta of 36 is
1.00000000001455192189104198423592963224531842098380889412403806\
9139542218571745865030220152998942329578185363084791339999779092\
891491916899
precision 20;
139
do!*one!*zeta(-17/3);
Zeta of - 5.6666666666666666667 is - 0.0018766468228592287697
do!*one!*zeta(z);
y 5 4 3 2 2 3 4 5
Zeta of x*(x + x + 5*x *y + 10*x *y + 10*x *y + 5*x*y + y ) is zeta(
y 6 5 4 2 3 3 2 4 5
x *x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + x*y )
y := 3;
y := 3
x := pi;
x := 3.1415926535897932385
do!*one!*zeta(z);
Zeta of 27548.203250393209469 is 1
do!*one!*zeta(190);
Zeta of 190 is 1.0
do!*one!*zeta(300);
Zeta of 300 is 1
do!*one!*zeta(0);
Zeta of 0 is - 0.5
do!*one!*zeta(-44);
Zeta of -44 is 0
off rounded;
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
% 8. Kummer functions
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
off rounded;
t!*kummer!*a := { 1, 2.4, -1397/10 }$
t!*kummer!*b := { 0, 1, pi, -pi, 26 }$
for each a in t!*kummer!*a do
for each b in t!*kummer!*a do
for each z in t!*kummer!*a do
<< write "KummerM(", a, ",", b, ",", z, ") = ",
kummerm(a,b,z);
write "KummerU(", a, ",", b, ",", z, ") = ",
kummeru(a,b,z) >>;
KummerM(1,1,1) = e
KummerU(1,1,1) = kummeru(1,1,1)
12 2/5 2
KummerM(1,1,----) = e *e
5
12 12
KummerU(1,1,----) = kummeru(1,1,----)
5 5
- 1397 1
KummerM(1,1,---------) = ------------
10 7/10 139
e *e
- 1397 - 1397
KummerU(1,1,---------) = kummeru(1,1,---------)
10 10
12 12
KummerM(1,----,1) = kummerm(1,----,1)
5 5
2 12 12
2*gamma(---)*(gamma(----)*e - kummerm(1,----,1))
12 5 5 5
KummerU(1,----,1) = --------------------------------------------------
5 12
5*gamma(----)
5
12 12 12 12
KummerM(1,----,----) = kummerm(1,----,----)
5 5 5 5
12 12
KummerU(1,----,----) =
5 5
2 2/5 2/5 12 2 4/5 2/5 12 12
gamma(---)*(5*e *5 *gamma(----)*e - 12*2 *3 *kummerm(1,----,----))
5 5 5 5
-----------------------------------------------------------------------------
4/5 2/5 12
30*2 *3 *gamma(----)
5
12 - 1397 12 - 1397
KummerM(1,----,---------) = kummerm(1,----,---------)
5 10 5 10
12 - 1397 2
KummerU(1,----,---------) = (2*gamma(---)*(
5 10 5
7/10 2/5 12 - 1397 139
- 1397*e *1397 *kummerm(1,----,---------)*e
5 10
2/5 12 7/10 2/5 12 139
- 10*10 *gamma(----)))/(6985*e *1397 *gamma(----)*e )
5 5
- 1397 - 1397
KummerM(1,---------,1) = kummerm(1,---------,1)
10 10
- 1397
KummerU(1,---------,1) =
10
1397 - 1397 3*pi
1397*gamma(------)*kummerm(1,---------,1)*sin(------) - 10*e*pi
10 10 10
-----------------------------------------------------------------
1417 3*pi
10*gamma(------)*sin(------)
10 10
- 1397 12 - 1397 12
KummerM(1,---------,----) = kummerm(1,---------,----)
10 5 10 5
- 1397 12
KummerU(1,---------,----) = ( - 4868940170769807140163734629614870951763530546\
10 5
631554678436906268603895538747067202821442117249242625797089485978560806001\
2/5 2/5 7/10 2
0798787187515790195560655355904*e *2 *3 *e *pi + 2004596689573651914\
834381033324557143448090555188728353940095378016295057932438794523477554321\
7/10 1397 - 1397 12 3*pi
2890625*5 *gamma(------)*kummerm(1,---------,----)*sin(------))/(1434929\
10 10 5 10
627468612680625899093288874118430988228481552150279237922703146068670321255\
7/10 1417 3*pi
92231750488281250*5 *gamma(------)*sin(------))
10 10
- 1397 - 1397 - 1397 - 1397
KummerM(1,---------,---------) = kummerm(1,---------,---------)
10 10 10 10
- 1397 - 1397
KummerU(1,---------,---------) = (1397*( - 15215998257931517360033101903862815\
10 10
223888543388957222106115555299202666279048867916168002482075600873075777\
814303704269392393284644411305647344643215427879407971076364582034085385\
160268928702508481521778353248095651747223572373118125240197785695254459\
776199438476999663624567921028135476855359579600382069871166959862260275\
271370315144729355759279591209107609434598035344309681771044690476129730\
7/10
6235868715700448137787799653833690839289133*( - 1397) *pi + 100000000\
000000000000000000000000000000000000000000000000000000000000000000000000\
7/10 7/10
00000000000000000000000000000000000000000000000000000000000*e *10
1397 - 1397 - 1397 3*pi 139
*gamma(------)*kummerm(1,---------,---------)*sin(------)*e ))/(100000\
10 10 10 10
000000000000000000000000000000000000000000000000000000000000000000000000000\
7/10 7/10
000000000000000000000000000000000000000000000000000000000000*e *10
1417 3*pi 139
*gamma(------)*sin(------)*e )
10 10
12 12
KummerM(----,1,1) = kummerm(----,1,1)
5 5
12 12
KummerU(----,1,1) = kummeru(----,1,1)
5 5
12 12 12 12
KummerM(----,1,----) = kummerm(----,1,----)
5 5 5 5
12 12 12 12
KummerU(----,1,----) = kummeru(----,1,----)
5 5 5 5
12 - 1397 12 - 1397
KummerM(----,1,---------) = kummerm(----,1,---------)
5 10 5 10
12 - 1397 12 - 1397
KummerU(----,1,---------) = kummeru(----,1,---------)
5 10 5 10
12 12
KummerM(----,----,1) = e
5 5
2 - 2 2*pi
2*gamma(---)*kummerm(1,------,1)*sin(------) + 5*e*pi
12 12 5 5 5
KummerU(----,----,1) = -------------------------------------------------------
5 5 12 2*pi
5*gamma(----)*sin(------)
5 5
12 12 12 2/5 2
KummerM(----,----,----) = e *e
5 5 5
12 12 12
KummerU(----,----,----) =
5 5 5
2/5 4/5 2/5 2 2/5 2 - 2 12 2*pi
6*e *2 *3 *e *pi + 5 *gamma(---)*kummerm(1,------,----)*sin(------)
5 5 5 5
-----------------------------------------------------------------------------
4/5 2/5 12 2*pi
6*2 *3 *gamma(----)*sin(------)
5 5
12 12 - 1397 1
KummerM(----,----,---------) = ------------
5 5 10 7/10 139
e *e
12 12 - 1397
KummerU(----,----,---------) = (
5 5 10
7/10 2/5 2 - 2 - 1397 2*pi 139
- 4*e *10 *gamma(---)*kummerm(1,------,---------)*sin(------)*e
5 5 10 5
2/5 7/10 2/5 12 2*pi 139
+ 1397*1397 *pi)/(1397*e *1397 *gamma(----)*sin(------)*e )
5 5
12 - 1397 12 - 1397
KummerM(----,---------,1) = kummerm(----,---------,1)
5 10 5 10
12 - 1397 1431 1431 1417
KummerU(----,---------,1) = ( - 10*gamma(------)*kummerm(------,------,1)*pi +
5 10 10 10 10
1417 1397 12 12 - 1397
1397*gamma(------)*gamma(------)*gamma(----)*kummerm(----,---------,1)
10 10 5 5 10
3*pi 1431 1417 12 3*pi
*sin(------))/(10*gamma(------)*gamma(------)*gamma(----)*sin(------))
10 10 10 5 10
12 - 1397 12 12 - 1397 12
KummerM(----,---------,----) = kummerm(----,---------,----)
5 10 5 5 10 5
12 - 1397 12
KummerU(----,---------,----) = (2004596689573651914834381033324557143448090555\
5 10 5
7/10 1417
1887283539400953780162950579324387945234775543212890625*5 *gamma(------)
10
1397 12 12 - 1397 12 3*pi
*gamma(------)*gamma(----)*kummerm(----,---------,----)*sin(------) - 48689\
10 5 5 10 5 10
401707698071401637346296148709517635305466315546784369062686038955387470672\
2/5
028214421172492426257970894859785608060010798787187515790195560655355904*2
7/10 1431 1431 1417 12
*3 *gamma(------)*kummerm(------,------,----)*pi)/(143492962746861268062\
10 10 10 5
589909328887411843098822848155215027923792270314606867032125592231750488281\
7/10 1431 1417 12 3*pi
250*5 *gamma(------)*gamma(------)*gamma(----)*sin(------))
10 10 5 10
12 - 1397 - 1397 12 - 1397 - 1397
KummerM(----,---------,---------) = kummerm(----,---------,---------)
5 10 10 5 10 10
12 - 1397 - 1397
KummerU(----,---------,---------) = (1397*( - 15215998257931517360033101903862\
5 10 10
815223888543388957222106115555299202666279048867916168002482075600873075\
777814303704269392393284644411305647344643215427879407971076364582034085\
385160268928702508481521778353248095651747223572373118125240197785695254\
459776199438476999663624567921028135476855359579600382069871166959862260\
275271370315144729355759279591209107609434598035344309681771044690476129\
7/10 1431
7306235868715700448137787799653833690839289133*( - 1397) *gamma(------)
10
1431 1417 - 1397
*kummerm(------,------,---------)*pi + 100000000000000000000000000000000\
10 10 10
000000000000000000000000000000000000000000000000000000000000000000000000\
7/10 1417 1397
00000000000000000000000000000000000*10 *gamma(------)*gamma(------)
10 10
12 12 - 1397 - 1397 3*pi
*gamma(----)*kummerm(----,---------,---------)*sin(------)))/(1000000000\
5 5 10 10 10
000000000000000000000000000000000000000000000000000000000000000000000000000\
7/10 1431
00000000000000000000000000000000000000000000000000000000*10 *gamma(------)
10
1417 12 3*pi
*gamma(------)*gamma(----)*sin(------))
10 5 10
- 1397 - 1397
KummerM(---------,1,1) = kummerm(---------,1,1)
10 10
- 1397 - 1397
KummerU(---------,1,1) = kummeru(---------,1,1)
10 10
- 1397 12 - 1397 12
KummerM(---------,1,----) = kummerm(---------,1,----)
10 5 10 5
- 1397 12 - 1397 12
KummerU(---------,1,----) = kummeru(---------,1,----)
10 5 10 5
- 1397 - 1397 - 1397 - 1397
KummerM(---------,1,---------) = kummerm(---------,1,---------)
10 10 10 10
- 1397 - 1397 - 1397 - 1397
KummerU(---------,1,---------) = kummeru(---------,1,---------)
10 10 10 10
- 1397 12 - 1397 12
KummerM(---------,----,1) = kummerm(---------,----,1)
10 5 10 5
- 1397 12
KummerU(---------,----,1) = (
10 5
1411 - 1397 12 pi
7055*gamma(------)*kummerm(---------,----,1)*sin(----)*pi + 2794
10 10 5 10
1397 12 2 - 1411 - 2 3*pi
*gamma(------)*gamma(----)*gamma(---)*kummerm(---------,------,1)*sin(------)
10 5 5 10 5 10
2*pi 12 2*pi
*sin(------))/(50*gamma(----)*sin(------)*pi)
5 5 5
- 1397 12 12 - 1397 12 12
KummerM(---------,----,----) = kummerm(---------,----,----)
10 5 5 10 5 5
- 1397 12 12 2/5 1397 12 2
KummerU(---------,----,----) = (1397*5 *gamma(------)*gamma(----)*gamma(---)
10 5 5 10 5 5
- 1411 - 2 12 3*pi 2*pi
*kummerm(---------,------,----)*sin(------)*sin(------)
10 5 5 10 5
4/5 2/5 1411 - 1397 12 12 pi
+ 8466*2 *3 *gamma(------)*kummerm(---------,----,----)*sin(----)*pi)/(
10 10 5 5 10
4/5 2/5 12 2*pi
60*2 *3 *gamma(----)*sin(------)*pi)
5 5
- 1397 12 - 1397 - 1397 12 - 1397
KummerM(---------,----,---------) = kummerm(---------,----,---------)
10 5 10 10 5 10
- 1397 12 - 1397
KummerU(---------,----,---------) = (
10 5 10
2/5 1411 - 1397 12 - 1397 pi
1411*1397 *gamma(------)*kummerm(---------,----,---------)*sin(----)*pi - 4
10 10 5 10 10
2/5 1397 12 2
*10 *gamma(------)*gamma(----)*gamma(---)
10 5 5
- 1411 - 2 - 1397 3*pi 2*pi 2/5
*kummerm(---------,------,---------)*sin(------)*sin(------))/(10*1397
10 5 10 10 5
12 2*pi
*gamma(----)*sin(------)*pi)
5 5
- 1397 - 1397
KummerM(---------,---------,1) = e
10 10
- 1397 - 1397
KummerU(---------,---------,1) =
10 10
1397 1417 1417
1397*gamma(------)*(gamma(------)*e - kummerm(1,------,1))
10 10 10
------------------------------------------------------------
1417
10*gamma(------)
10
- 1397 - 1397 12 2/5 2
KummerM(---------,---------,----) = e *e
10 10 5
- 1397 - 1397 12 1397
KummerU(---------,---------,----) = (1397*gamma(------)*(717464813734306340312\
10 10 5 10
949546644437059215494114240776075139618961351573034335160627961158752441\
2/5 7/10 1417 2
40625*e *5 *gamma(------)*e - 243447008538490357008186731480743547\
10
588176527331577733921845313430194776937353360141072105862462131289854474\
2/5 7/10
29892804030005399393593757895097780327677952*2 *3
1417 12
*kummerm(1,------,----)))/(717464813734306340312949546644437059215494114\
10 5
7/10 1417
240776075139618961351573034335160627961158752441406250*5 *gamma(------))
10
- 1397 - 1397 - 1397 1
KummerM(---------,---------,---------) = ------------
10 10 10 7/10 139
e *e
- 1397 - 1397 - 1397 1397
KummerU(---------,---------,---------) = (1397*gamma(------)*( - 2125674956633\
10 10 10 10
032975196624335969635286777229511437323928224343075298612479183126847888\
669946745961441968686160658227486434117341864824259398934046657195274753\
293559368132110161728306889569339740434868592435948758962549087133060524\
602096055630661627048030735061555236853008352138567630526116693733270173\
375161002024292757760455410433025718690999571358891912333038013345537600\
7/10
06254341494325951532336811508595833526048489556116405666102486918801*e
7/10 1417 - 1397 139
*( - 1397) *kummerm(1,------,---------)*e + 1000000000000000000000\
10 10
000000000000000000000000000000000000000000000000000000000000000000000000\
7/10 1417
00000000000000000000000000000000000000000000000*10 *gamma(------)))/(
10
100000000000000000000000000000000000000000000000000000000000000000000000000\
7/10
0000000000000000000000000000000000000000000000000000000000000000000*e
7/10 1417 139
*10 *gamma(------)*e )
10
on rounded;
precision 30;
20
t!*k!*c := 7;
t*k*c := 7
% To test each and every possible combination of
% three arguments from t!*kummer!*b would take too
% long, but we want the possibility of trying most
% special cases. Compromise: test every seventh
% possibility.
for each a in t!*kummer!*b do
for each b in t!*kummer!*b do
for each z in t!*kummer!*b do
<< if t!*k!*c = 7
then << write "KummerM(", a, ",", b, ",", z, ") = ",
kummerm(a,b,z);
write "KummerU(", a, ",", b, ",", z, ") = ",
kummeru(a,b,z);
t!*k!*c := 0 >>;
t!*k!*c := t!*k!*c + 1 >>;
KummerM(0,0,0) = kummerm(0,0,0)
KummerU(0,0,0) = kummeru(0,0,0)
KummerM(0,1,3.14159265358979323846264338328) = 1
KummerU(0,1,3.14159265358979323846264338328) = kummeru(0,1,
3.14159265358979323846264338328)
KummerM(0,3.14159265358979323846264338328,26) = 1
KummerU(0,3.14159265358979323846264338328,26) = kummeru(0,
3.14159265358979323846264338328,26)
KummerM(0,26,1) = 1
KummerU(0,26,1) = kummeru(0,26,1)
KummerM(1,0, - 3.14159265358979323846264338328) = kummerm(1,0,
- 3.14159265358979323846264338328)
KummerU(1,0, - 3.14159265358979323846264338328) = kummeru(1,0,
- 3.14159265358979323846264338328)
KummerM(1,3.14159265358979323846264338328,0) = 1
KummerU(1,3.14159265358979323846264338328,0) = kummeru(1,
3.14159265358979323846264338328,0)
KummerM(1, - 3.14159265358979323846264338328,3.14159265358979323846264338328) =
2692.89480079631357528203659153
KummerU(1, - 3.14159265358979323846264338328,3.14159265358979323846264338328) =
0.129554194296952806409644905659
KummerM(1,26,26) = 7.74565667206271943920803547594
KummerU(1,26,26) = kummeru(1,26,26)
KummerM(3.14159265358979323846264338328,1,1) = 10.2259987795328570162092950355
KummerU(3.14159265358979323846264338328,1,1) = kummeru(
3.14159265358979323846264338328,1,1)
KummerM(3.14159265358979323846264338328,3.14159265358979323846264338328,
- 3.14159265358979323846264338328) = 0.0432139182637722497744177371717
KummerU(3.14159265358979323846264338328,3.14159265358979323846264338328,
- 3.14159265358979323846264338328) = ( - 0.137891580772667438638357954721
2.14159265358979323846264338328
*( - 3.14159265358979323846264338328)
+ 1.22938452238615186004844111454)/
2.14159265358979323846264338328
( - 3.14159265358979323846264338328)
KummerM(3.14159265358979323846264338328,26,0) = 1
KummerU(3.14159265358979323846264338328,26,0) = kummeru(
3.14159265358979323846264338328,26,0)
KummerM( - 3.14159265358979323846264338328,0,3.14159265358979323846264338328) =
kummerm( - 3.14159265358979323846264338328,0,3.14159265358979323846264338328)
KummerU( - 3.14159265358979323846264338328,0,3.14159265358979323846264338328) =
kummeru( - 3.14159265358979323846264338328,0,3.14159265358979323846264338328)
KummerM( - 3.14159265358979323846264338328,1,26) =
6.18522226564722800173513559462e+5
KummerU( - 3.14159265358979323846264338328,1,26) = kummeru(
- 3.14159265358979323846264338328,1,26)
KummerM( - 3.14159265358979323846264338328, - 3.14159265358979323846264338328,1
) = 2.71828182845904523536028747135
KummerU( - 3.14159265358979323846264338328, - 3.14159265358979323846264338328,1
) = 19.2419564406028465661346237348
KummerM( - 3.14159265358979323846264338328,26, - 3.14159265358979323846264338328
) = 1.42892253084220157246185820464
KummerU( - 3.14159265358979323846264338328,26, - 3.14159265358979323846264338328
) = kummeru( - 3.14159265358979323846264338328,26,
- 3.14159265358979323846264338328)
KummerM(26,1,0) = 1
KummerU(26,1,0) = kummeru(26,1,0)
KummerM(26,3.14159265358979323846264338328,3.14159265358979323846264338328) =
3.91638029828702661403357541917e+5
KummerU(26,3.14159265358979323846264338328,3.14159265358979323846264338328) =
4.18995913372628050180833640475e-32
KummerM(26, - 3.14159265358979323846264338328,26) =
2.19471945265322419268333614674e+34
KummerU(26, - 3.14159265358979323846264338328,26) = 0
off rounded;
clear x, y, z, t!*k!*c;
df(df(kummerM(a,b,z),z),z);
2
(kummerm(2,3,z)*z - 2*kummerm(2,3,z)*z + 2*kummerm(2,3,z) + 2*kummerm(1,3,z)*z
2
- 4*kummerm(1,3,z) + 2*kummerm(0,3,z))/z
df(kummerU(a,b,z),z);
kummeru(2,3,z)*z - kummeru(2,3,z) - kummeru(1,3,z)
----------------------------------------------------
z
z := ((x^2 + y)^5) + (x^(x+y));
x + y 10 8 6 2 4 3 2 4 5
z := x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y
df(df(kummerM(a,b,z),y),x);
4*x + 4*y
(x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 4*x + 4*y
*log(x) *x + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4*x + 4*y
*x + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3*x + 3*y
*y + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 11 3*x + 3*y
*log(x) *x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 9 3*x + 3*y
*log(x) *x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 7 2 3*x + 3*y
*log(x) *x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 5 3 3*x + 3*y
*log(x) *x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 3 4 3*x + 3*y
*log(x) *x *y + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 5 3*x + 3*y
*log(x) *x*y - x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 3*x + 3*y
*log(x) *x + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 3*x + 3*y
*x + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 3*x + 3*y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 3*x + 3*y
*x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 3*x + 3*y
*x *y + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 3*x + 3*y
*x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 2 3*x + 3*y
*x *y + 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 3*x + 3*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 2 3*x + 3*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 3*x + 3*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 3 3*x + 3*y
*x *y + 60*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 2 3*x + 3*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 3 3*x + 3*y
*x *y + 30*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 2 3*x + 3*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 4 3*x + 3*y
*x *y + 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 3 3*x + 3*y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 4 3*x + 3*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 3 3*x + 3*y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 5 3*x + 3*y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 4 3*x + 3*y
*x *y + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 3*x + 3*y
*x*y + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 3*x + 3*y
*x*y - x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3*x + 3*y
*x + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 3*x + 3*y
*y - x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3*x + 3*y
*y + 5*x
x + y 10 8 6 2 4 3 2 4 5 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x + 5
3*x + 3*y
*x
x + y 10 8 6 2 4 3 2 4 5 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
3*x + 3*y
20*x
x + y 10 8 6 2 4 3 2 4 5 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
3*x + 3*y
20*x
x + y 10 8 6 2 4 3 2 4 5 6 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
3*x + 3*y
+ 30*x
x + y 10 8 6 2 4 3 2 4 5 5 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
3*x + 3*y
+ 30*x
x + y 10 8 6 2 4 3 2 4 5 4 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
3*x + 3*y
+ 20*x
x + y 10 8 6 2 4 3 2 4 5 3 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
3*x + 3*y
+ 20*x
x + y 10 8 6 2 4 3 2 4 5 2 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
3*x + 3*y
+ 5*x
x + y 10 8 6 2 4 3 2 4 5 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y +
3*x + 3*y
5*x
x + y 10 8 6 2 4 3 2 4 5 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y +
3*x + 3*y
x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y ) + 2
3*x + 3*y
*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 3*x + 3*y
*log(x) *x + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3*x + 3*y
*x + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2*x + 2*y
*y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 21 2*x + 2*y
*log(x) *x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 19 2*x + 2*y
*log(x) *x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 17 2 2*x + 2*y
*log(x) *x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 15 3 2*x + 2*y
*log(x) *x *y + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 13 4 2*x + 2*y
*log(x) *x *y + 252*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 11 5 2*x + 2*y
*log(x) *x *y + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 9 6 2*x + 2*y
*log(x) *x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 7 7 2*x + 2*y
*log(x) *x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 5 8 2*x + 2*y
*log(x) *x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 3 9 2*x + 2*y
*log(x) *x *y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 10 2*x + 2*y
*log(x) *x*y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 2*x + 2*y
*log(x) *x + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
21 2*x + 2*y
*x + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
20 2*x + 2*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
20 2*x + 2*y
*x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
19 2*x + 2*y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
19 2*x + 2*y
*x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
18 2 2*x + 2*y
*x *y + 180*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
18 2*x + 2*y
*x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
17 2 2*x + 2*y
*x *y + 90*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
17 2*x + 2*y
*x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
16 3 2*x + 2*y
*x *y + 720*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
16 2 2*x + 2*y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
15 3 2*x + 2*y
*x *y + 360*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
15 2 2*x + 2*y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
14 4 2*x + 2*y
*x *y + 1680*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
14 3 2*x + 2*y
*x *y + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
13 4 2*x + 2*y
*x *y + 840*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
13 3 2*x + 2*y
*x *y + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
12 5 2*x + 2*y
*x *y + 2520*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
12 4 2*x + 2*y
*x *y + 252*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 5 2*x + 2*y
*x *y + 1260*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 4 2*x + 2*y
*x *y + 252*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 6 2*x + 2*y
*x *y + 2520*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 5 2*x + 2*y
*x *y - 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 2*x + 2*y
*x + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 6 2*x + 2*y
*x *y + 1260*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 5 2*x + 2*y
*x *y - 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 2*x + 2*y
*x + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 7 2*x + 2*y
*x *y + 1680*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 6 2*x + 2*y
*x *y - 80*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 2*x + 2*y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 7 2*x + 2*y
*x *y + 840*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 6 2*x + 2*y
*x *y - 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 2*x + 2*y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 8 2*x + 2*y
*x *y + 720*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 7 2*x + 2*y
*x *y - 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 2 2*x + 2*y
*x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 8 2*x + 2*y
*x *y + 360*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 7 2*x + 2*y
*x *y - 60*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 2 2*x + 2*y
*x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 9 2*x + 2*y
*x *y + 180*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 8 2*x + 2*y
*x *y - 80*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 3 2*x + 2*y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 9 2*x + 2*y
*x *y + 90*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 8 2*x + 2*y
*x *y - 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 3 2*x + 2*y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 10 2*x + 2*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 9 2*x + 2*y
*x *y - 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 4 2*x + 2*y
*x *y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 2*x + 2*y
*x*y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 2*x + 2*y
*x*y - 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 2*x + 2*y
*x*y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2*x + 2*y
*x + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 2*x + 2*y
*y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2*x + 2*y
*y + 10*x
x + y 10 8 6 2 4 3 2 4 5 19
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
2*x + 2*y
10*x
x + y 10 8 6 2 4 3 2 4 5 18
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 50*x
x + y 10 8 6 2 4 3 2 4 5 18
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
2*x + 2*y
90*x
x + y 10 8 6 2 4 3 2 4 5 17
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 90*x
x + y 10 8 6 2 4 3 2 4 5 16 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 400*x
x + y 10 8 6 2 4 3 2 4 5 16
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 360*x
x + y 10 8 6 2 4 3 2 4 5 15 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 360*x
x + y 10 8 6 2 4 3 2 4 5 14 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 1400*x
x + y 10 8 6 2 4 3 2 4 5 14 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 840*x
x + y 10 8 6 2 4 3 2 4 5 13 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 840*x
x + y 10 8 6 2 4 3 2 4 5 12 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 2800*x
x + y 10 8 6 2 4 3 2 4 5 12 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 1260*x
x + y 10 8 6 2 4 3 2 4 5 11 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 1260*x
x + y 10 8 6 2 4 3 2 4 5 10 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 3500*x
x + y 10 8 6 2 4 3 2 4 5 10 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 2*x
x + y 10 8 6 2 4 3 2 4 5 10
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
2*x + 2*y
1260*x
x + y 10 8 6 2 4 3 2 4 5 9 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
- 10*x
x + y 10 8 6 2 4 3 2 4 5 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
2*x + 2*y
1260*x
x + y 10 8 6 2 4 3 2 4 5 8 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 2800*x
x + y 10 8 6 2 4 3 2 4 5 8 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
2*x + 2*y
840*x
x + y 10 8 6 2 4 3 2 4 5 7 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
- 40*x
x + y 10 8 6 2 4 3 2 4 5 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
2*x + 2*y
840*x
x + y 10 8 6 2 4 3 2 4 5 6 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 1400*x
x + y 10 8 6 2 4 3 2 4 5 6 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
- 20*x
x + y 10 8 6 2 4 3 2 4 5 6 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 120*x
x + y 10 8 6 2 4 3 2 4 5 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
2*x + 2*y
360*x
x + y 10 8 6 2 4 3 2 4 5 5 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
- 60*x
x + y 10 8 6 2 4 3 2 4 5 5 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 360*x
x + y 10 8 6 2 4 3 2 4 5 4 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 400*x
x + y 10 8 6 2 4 3 2 4 5 4 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
- 40*x
x + y 10 8 6 2 4 3 2 4 5 4 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 120*x
x + y 10 8 6 2 4 3 2 4 5 4 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 90*x
x + y 10 8 6 2 4 3 2 4 5 3 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
- 40*x
x + y 10 8 6 2 4 3 2 4 5 3 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 90*x
x + y 10 8 6 2 4 3 2 4 5 2 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 50*x
x + y 10 8 6 2 4 3 2 4 5 2 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
- 30*x
x + y 10 8 6 2 4 3 2 4 5 2 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 2 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 10*x
x + y 10 8 6 2 4 3 2 4 5 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y -
2*x + 2*y
10*x
x + y 10 8 6 2 4 3 2 4 5 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y +
2*x + 2*y
10*x
x + y 10 8 6 2 4 3 2 4 5 10
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y -
2*x + 2*y
8*x
x + y 10 8 6 2 4 3 2 4 5 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y -
2*x + 2*y
x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y ) + 2
2*x + 2*y
*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 11 2*x + 2*y
*log(x) *x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 9 2*x + 2*y
*log(x) *x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 7 2 2*x + 2*y
*log(x) *x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 5 3 2*x + 2*y
*log(x) *x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 3 4 2*x + 2*y
*log(x) *x *y + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 5 2*x + 2*y
*log(x) *x*y - 3*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 2*x + 2*y
*log(x) *x + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 2*x + 2*y
*x + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 2*x + 2*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 2*x + 2*y
*x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 2*x + 2*y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 2*x + 2*y
*x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 2 2*x + 2*y
*x *y + 80*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 2*x + 2*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 2 2*x + 2*y
*x *y + 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 2*x + 2*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 3 2*x + 2*y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 2 2*x + 2*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 3 2*x + 2*y
*x *y + 60*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 2 2*x + 2*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 4 2*x + 2*y
*x *y + 80*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 3 2*x + 2*y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 4 2*x + 2*y
*x *y + 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 3 2*x + 2*y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 5 2*x + 2*y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 4 2*x + 2*y
*x *y + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 2*x + 2*y
*x*y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 2*x + 2*y
*x*y - 3*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2*x + 2*y
*x + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 2*x + 2*y
*y - 3*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2*x + 2*y
*y + 10*x
x + y 10 8 6 2 4 3 2 4 5 9
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
2*x + 2*y
10*x
x + y 10 8 6 2 4 3 2 4 5 8
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
2*x + 2*y
40*x
x + y 10 8 6 2 4 3 2 4 5 7
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
2*x + 2*y
40*x
x + y 10 8 6 2 4 3 2 4 5 6 2
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 60*x
x + y 10 8 6 2 4 3 2 4 5 5 2
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 60*x
x + y 10 8 6 2 4 3 2 4 5 4 3
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 3 3
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 2 4
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
2*x + 2*y
+ 10*x
x + y 10 8 6 2 4 3 2 4 5 4
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y +
2*x + 2*y
10*x
x + y 10 8 6 2 4 3 2 4 5 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y +
2*x + 2*y
x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y ) + 2
2*x + 2*y
*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 2*x + 2*y
*log(x) *x + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2*x + 2*y
*x + 2*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
x + y
*y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 21 x + y
*log(x) *x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 19 x + y
*log(x) *x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 17 2 x + y
*log(x) *x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 15 3 x + y
*log(x) *x *y + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 13 4 x + y
*log(x) *x *y + 252*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 11 5 x + y
*log(x) *x *y - x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 11 x + y
*log(x) *x + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 9 6 x + y
*log(x) *x *y - 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 9 x + y
*log(x) *x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 7 7 x + y
*log(x) *x *y - 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 7 2 x + y
*log(x) *x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 5 8 x + y
*log(x) *x *y - 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 5 3 x + y
*log(x) *x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 3 9 x + y
*log(x) *x *y - 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 3 4 x + y
*log(x) *x *y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 10 x + y
*log(x) *x*y - x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 5 x + y
*log(x) *x*y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
30 x + y
*x + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
29 x + y
*x + 140*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
28 x + y
*x *y + 70*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
27 x + y
*x *y + 910*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
26 2 x + y
*x *y + 455*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
25 2 x + y
*x *y + 3640*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
24 3 x + y
*x *y + 1820*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
23 3 x + y
*x *y + 10010*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
22 4 x + y
*x *y + 5005*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
21 4 x + y
*x *y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
21 x + y
*x + 20020*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
20 5 x + y
*x *y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
20 x + y
*x *y - 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
20 x + y
*x + 10010*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
19 5 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
19 x + y
*x *y - 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
19 x + y
*x + 30030*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
18 6 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
18 2 x + y
*x *y - 180*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
18 x + y
*x *y + 15015*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
17 6 x + y
*x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
17 2 x + y
*x *y - 90*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
17 x + y
*x *y + 34320*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
16 7 x + y
*x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
16 3 x + y
*x *y - 720*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
16 2 x + y
*x *y + 17160*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
15 7 x + y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
15 3 x + y
*x *y - 360*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
15 2 x + y
*x *y + 30030*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
14 8 x + y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
14 4 x + y
*x *y - 1680*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
14 3 x + y
*x *y + 15015*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
13 8 x + y
*x *y + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
13 4 x + y
*x *y - 840*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
13 3 x + y
*x *y + 20020*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
12 9 x + y
*x *y + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
12 5 x + y
*x *y - 2520*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
12 4 x + y
*x *y + 10010*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 9 x + y
*x *y + 252*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 5 x + y
*x *y - 1260*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 4 x + y
*x *y - x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 x + y
*x + 10010*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 10 x + y
*x *y + 252*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 6 x + y
*x *y - 2520*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 5 x + y
*x *y - x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 x + y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 x + y
*x + 5005*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 10 x + y
*x *y + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 6 x + y
*x *y - 1260*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 5 x + y
*x *y - 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 x + y
*x + 3640*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 11 x + y
*x *y + 210*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 7 x + y
*x *y - 1680*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 6 x + y
*x *y - 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 2 x + y
*x *y + 80*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 x + y
*x *y + 1820*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 11 x + y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 7 x + y
*x *y - 840*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 6 x + y
*x *y - 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 2 x + y
*x *y + 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 x + y
*x *y + 910*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 12 x + y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 8 x + y
*x *y - 720*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 7 x + y
*x *y - 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 3 x + y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 2 x + y
*x *y + 455*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 12 x + y
*x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 8 x + y
*x *y - 360*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 7 x + y
*x *y - 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 3 x + y
*x *y + 60*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 2 x + y
*x *y + 140*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 13 x + y
*x *y + 45*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 9 x + y
*x *y - 180*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 8 x + y
*x *y - 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 4 x + y
*x *y + 80*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 3 x + y
*x *y + 70*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 13 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 9 x + y
*x *y - 90*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 8 x + y
*x *y - 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 4 x + y
*x *y + 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 3 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 14 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 10 x + y
*x *y - 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 9 x + y
*x *y - 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 5 x + y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 4 x + y
*x *y + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
14 x + y
*x*y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 x + y
*x*y - 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 x + y
*x*y - x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 x + y
*x*y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 x + y
*x*y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 x + y
*y - x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 x + y
*y + 5*x
x + y 10 8 6 2 4 3 2 4 5 29
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
5*x
x + y 10 8 6 2 4 3 2 4 5 28
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 100*x
x + y 10 8 6 2 4 3 2 4 5 28
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
70*x
x + y 10 8 6 2 4 3 2 4 5 27
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 70*x
x + y 10 8 6 2 4 3 2 4 5 26 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 1300*x
x + y 10 8 6 2 4 3 2 4 5 26
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 455*x
x + y 10 8 6 2 4 3 2 4 5 25 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 455*x
x + y 10 8 6 2 4 3 2 4 5 24 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 7800*x
x + y 10 8 6 2 4 3 2 4 5 24 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 1820*x
x + y 10 8 6 2 4 3 2 4 5 23 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 1820*x
x + y 10 8 6 2 4 3 2 4 5 22 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 28600*x
x + y 10 8 6 2 4 3 2 4 5 22 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 5005*x
x + y 10 8 6 2 4 3 2 4 5 21 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 5005*x
x + y 10 8 6 2 4 3 2 4 5 20 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 71500*x
x + y 10 8 6 2 4 3 2 4 5 20 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ x
x + y 10 8 6 2 4 3 2 4 5 20
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
10010*x
x + y 10 8 6 2 4 3 2 4 5 19 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 10*x
x + y 10 8 6 2 4 3 2 4 5 19
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
10010*x
x + y 10 8 6 2 4 3 2 4 5 18 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 128700*x
x + y 10 8 6 2 4 3 2 4 5 18 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 20*x
x + y 10 8 6 2 4 3 2 4 5 18
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
15015*x
x + y 10 8 6 2 4 3 2 4 5 17 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 90*x
x + y 10 8 6 2 4 3 2 4 5 17
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 15015*x
x + y 10 8 6 2 4 3 2 4 5 16 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 171600*x
x + y 10 8 6 2 4 3 2 4 5 16 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 45*x
x + y 10 8 6 2 4 3 2 4 5 16 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 160*x
x + y 10 8 6 2 4 3 2 4 5 16
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 17160*x
x + y 10 8 6 2 4 3 2 4 5 15 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 360*x
x + y 10 8 6 2 4 3 2 4 5 15 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 17160*x
x + y 10 8 6 2 4 3 2 4 5 14 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 171600*x
x + y 10 8 6 2 4 3 2 4 5 14 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 240*x
x + y 10 8 6 2 4 3 2 4 5 14 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 560*x
x + y 10 8 6 2 4 3 2 4 5 14 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 15015*x
x + y 10 8 6 2 4 3 2 4 5 13 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 840*x
x + y 10 8 6 2 4 3 2 4 5 13 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 15015*x
x + y 10 8 6 2 4 3 2 4 5 12 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 128700*x
x + y 10 8 6 2 4 3 2 4 5 12 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 630*x
x + y 10 8 6 2 4 3 2 4 5 12 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 1120*x
x + y 10 8 6 2 4 3 2 4 5 12 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 10010*x
x + y 10 8 6 2 4 3 2 4 5 11 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 1260*x
x + y 10 8 6 2 4 3 2 4 5 11 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 10010*x
x + y 10 8 6 2 4 3 2 4 5 10
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
10 x + y
*y + 71500*x
x + y 10 8 6 2 4 3 2 4 5 10 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 1008*x
x + y 10 8 6 2 4 3 2 4 5 10 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 1400*x
x + y 10 8 6 2 4 3 2 4 5 10 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- x
x + y 10 8 6 2 4 3 2 4 5 10
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
5005*x
x + y 10 8 6 2 4 3 2 4 5 9 10
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 1260*x
x + y 10 8 6 2 4 3 2 4 5 9 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 10*x
x + y 10 8 6 2 4 3 2 4 5 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
5005*x
x + y 10 8 6 2 4 3 2 4 5 8 11
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 28600*x
x + y 10 8 6 2 4 3 2 4 5 8 10
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 1050*x
x + y 10 8 6 2 4 3 2 4 5 8 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 1120*x
x + y 10 8 6 2 4 3 2 4 5 8 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 5*x
x + y 10 8 6 2 4 3 2 4 5 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y -
x + y
40*x
x + y 10 8 6 2 4 3 2 4 5 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
1820*x
x + y 10 8 6 2 4 3 2 4 5 7 11
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 840*x
x + y 10 8 6 2 4 3 2 4 5 7 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
x + y
1820*x
x + y 10 8 6 2 4 3 2 4 5 6 12
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 7800*x
x + y 10 8 6 2 4 3 2 4 5 6 11
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 720*x
x + y 10 8 6 2 4 3 2 4 5 6 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 560*x
x + y 10 8 6 2 4 3 2 4 5 6 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 30*x
x + y 10 8 6 2 4 3 2 4 5 6 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 120*x
x + y 10 8 6 2 4 3 2 4 5 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
x + y
455*x
x + y 10 8 6 2 4 3 2 4 5 5 12
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 360*x
x + y 10 8 6 2 4 3 2 4 5 5 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 60*x
x + y 10 8 6 2 4 3 2 4 5 5 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 455*x
x + y 10 8 6 2 4 3 2 4 5 4 13
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 1300*x
x + y 10 8 6 2 4 3 2 4 5 4 12
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 315*x
x + y 10 8 6 2 4 3 2 4 5 4 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 160*x
x + y 10 8 6 2 4 3 2 4 5 4 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 50*x
x + y 10 8 6 2 4 3 2 4 5 4 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 120*x
x + y 10 8 6 2 4 3 2 4 5 4 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 70*x
x + y 10 8 6 2 4 3 2 4 5 3 13
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 90*x
x + y 10 8 6 2 4 3 2 4 5 3 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 3 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 70*x
x + y 10 8 6 2 4 3 2 4 5 2 14
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 100*x
x + y 10 8 6 2 4 3 2 4 5 2 13
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 80*x
x + y 10 8 6 2 4 3 2 4 5 2 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 20*x
x + y 10 8 6 2 4 3 2 4 5 2 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 35*x
x + y 10 8 6 2 4 3 2 4 5 2 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 40*x
x + y 10 8 6 2 4 3 2 4 5 2 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 5*x
x + y 10 8 6 2 4 3 2 4 5 14
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y
x + y
- 10*x
x + y 10 8 6 2 4 3 2 4 5 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y +
x + y
10*x
x + y 10 8 6 2 4 3 2 4 5 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y +
x + y
5*x
x + y 10 8 6 2 4 3 2 4 5 15
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y -
x + y
9*x
x + y 10 8 6 2 4 3 2 4 5 10
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y +
x + y
9*x
x + y 10 8 6 2 4 3 2 4 5 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y +
x + y x + y 10 8 6 2 4 3 2 4 5
x *kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 11 x + y
*log(x) *x + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 9 x + y
*log(x) *x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 7 2 x + y
*log(x) *x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 5 3 x + y
*log(x) *x *y + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 3 4 x + y
*log(x) *x *y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 5 x + y
*log(x) *x*y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
20 x + y
*x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
19 x + y
*x + 180*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
18 x + y
*x *y + 90*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
17 x + y
*x *y + 720*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
16 2 x + y
*x *y + 360*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
15 2 x + y
*x *y + 1680*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
14 3 x + y
*x *y + 840*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
13 3 x + y
*x *y + 2520*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
12 4 x + y
*x *y + 1260*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 4 x + y
*x *y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
11 x + y
*x + 2520*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 5 x + y
*x *y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 x + y
*x *y - 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 x + y
*x + 1260*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 5 x + y
*x *y + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 x + y
*x *y - 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 x + y
*x + 1680*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 6 x + y
*x *y + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 2 x + y
*x *y - 160*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 x + y
*x *y + 840*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 6 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 2 x + y
*x *y - 80*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 x + y
*x *y + 720*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 7 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 3 x + y
*x *y - 240*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 2 x + y
*x *y + 360*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 7 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 3 x + y
*x *y - 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 2 x + y
*x *y + 180*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 8 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 4 x + y
*x *y - 160*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 3 x + y
*x *y + 90*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 8 x + y
*x *y + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 4 x + y
*x *y - 80*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 3 x + y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 9 x + y
*x *y + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 5 x + y
*x *y - 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 4 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 x + y
*x*y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 x + y
*x*y - 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 x + y
*x*y + x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 x + y
*y + 10*x
x + y 10 8 6 2 4 3 2 4 5 19
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
10*x
x + y 10 8 6 2 4 3 2 4 5 18
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 100*x
x + y 10 8 6 2 4 3 2 4 5 18
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
90*x
x + y 10 8 6 2 4 3 2 4 5 17
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 90*x
x + y 10 8 6 2 4 3 2 4 5 16 2
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 800*x
x + y 10 8 6 2 4 3 2 4 5 16
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 360*x
x + y 10 8 6 2 4 3 2 4 5 15 2
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 360*x
x + y 10 8 6 2 4 3 2 4 5 14 3
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 2800*x
x + y 10 8 6 2 4 3 2 4 5 14 2
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 840*x
x + y 10 8 6 2 4 3 2 4 5 13 3
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 840*x
x + y 10 8 6 2 4 3 2 4 5 12 4
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 5600*x
x + y 10 8 6 2 4 3 2 4 5 12 3
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 1260*x
x + y 10 8 6 2 4 3 2 4 5 11 4
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 1260*x
x + y 10 8 6 2 4 3 2 4 5 10 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 7000*x
x + y 10 8 6 2 4 3 2 4 5 10 4
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ x
x + y 10 8 6 2 4 3 2 4 5 10
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
1260*x
x + y 10 8 6 2 4 3 2 4 5 9 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 20*x
x + y 10 8 6 2 4 3 2 4 5 9
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
1260*x
x + y 10 8 6 2 4 3 2 4 5 8 6
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 5600*x
x + y 10 8 6 2 4 3 2 4 5 8 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 15*x
x + y 10 8 6 2 4 3 2 4 5 8
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
x + y
40*x
x + y 10 8 6 2 4 3 2 4 5 8
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
840*x
x + y 10 8 6 2 4 3 2 4 5 7 6
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 80*x
x + y 10 8 6 2 4 3 2 4 5 7
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
x + y
840*x
x + y 10 8 6 2 4 3 2 4 5 6 7
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 2800*x
x + y 10 8 6 2 4 3 2 4 5 6 6
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 70*x
x + y 10 8 6 2 4 3 2 4 5 6 2
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 120*x
x + y 10 8 6 2 4 3 2 4 5 6
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
x + y
360*x
x + y 10 8 6 2 4 3 2 4 5 5 7
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 120*x
x + y 10 8 6 2 4 3 2 4 5 5 2
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 360*x
x + y 10 8 6 2 4 3 2 4 5 4 8
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 800*x
x + y 10 8 6 2 4 3 2 4 5 4 7
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 110*x
x + y 10 8 6 2 4 3 2 4 5 4 3
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 120*x
x + y 10 8 6 2 4 3 2 4 5 4 2
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 90*x
x + y 10 8 6 2 4 3 2 4 5 3 8
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 80*x
x + y 10 8 6 2 4 3 2 4 5 3 3
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 90*x
x + y 10 8 6 2 4 3 2 4 5 2 9
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 100*x
x + y 10 8 6 2 4 3 2 4 5 2 8
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
- 75*x
x + y 10 8 6 2 4 3 2 4 5 2 4
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 2 3
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 10*x
x + y 10 8 6 2 4 3 2 4 5 9
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y -
x + y
20*x
x + y 10 8 6 2 4 3 2 4 5 4
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y +
x + y
10*x
x + y 10 8 6 2 4 3 2 4 5 10
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y -
x + y
19*x
x + y 10 8 6 2 4 3 2 4 5 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y +
x + y
20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
10 x + y
*x + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 x + y
*x + 80*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
8 x + y
*x *y + 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 x + y
*x *y + 120*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
6 2 x + y
*x *y + 60*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 2 x + y
*x *y + 80*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 3 x + y
*x *y + 40*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 3 x + y
*x *y + 20*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2 4 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
4 x + y
*x*y + 10*x
x + y 10 8 6 2 4 3 2 4 5 9
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
10*x
x + y 10 8 6 2 4 3 2 4 5 8
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
x + y
40*x
x + y 10 8 6 2 4 3 2 4 5 7
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
x + y
40*x
x + y 10 8 6 2 4 3 2 4 5 6 2
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 60*x
x + y 10 8 6 2 4 3 2 4 5 5 2
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 60*x
x + y 10 8 6 2 4 3 2 4 5 4 3
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 3 3
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 2 4
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 10*x
x + y 10 8 6 2 4 3 2 4 5 4
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y +
x + y
10*x
x + y 10 8 6 2 4 3 2 4 5 5
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y +
x + y 10 8 6 2 4 3 2 4 5 38
50*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
x + y 10 8 6 2 4 3 2 4 5
+ 900*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
36
*x *y + 7650
x + y 10 8 6 2 4 3 2 4 5 34 2
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
+ 40800
x + y 10 8 6 2 4 3 2 4 5 32 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
+ 153000
x + y 10 8 6 2 4 3 2 4 5 30 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
+ 428400
x + y 10 8 6 2 4 3 2 4 5 28 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 60*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
28
*x + 928200
x + y 10 8 6 2 4 3 2 4 5 26 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 780*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
26
*x *y + 1591200
x + y 10 8 6 2 4 3 2 4 5 24 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 4680*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
24 2
*x *y + 2187900
x + y 10 8 6 2 4 3 2 4 5 22 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
- 17160
x + y 10 8 6 2 4 3 2 4 5 22 3
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
+ 2431000
x + y 10 8 6 2 4 3 2 4 5 20 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
- 42900
x + y 10 8 6 2 4 3 2 4 5 20 4
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
+ 2187900
x + y 10 8 6 2 4 3 2 4 5 18
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
10
*y - 77220
x + y 10 8 6 2 4 3 2 4 5 18 5
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 60*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
18
*x + 1591200
x + y 10 8 6 2 4 3 2 4 5 16
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
11
*y - 102960
x + y 10 8 6 2 4 3 2 4 5 16 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 480*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
16
*x *y + 928200
x + y 10 8 6 2 4 3 2 4 5 14
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
12
*y - 102960
x + y 10 8 6 2 4 3 2 4 5 14 7
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 1680*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
14 2
*x *y + 428400
x + y 10 8 6 2 4 3 2 4 5 12
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
13
*y - 77220
x + y 10 8 6 2 4 3 2 4 5 12 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 3360*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
12 3
*x *y + 153000
x + y 10 8 6 2 4 3 2 4 5 10
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
14
*y - 42900
x + y 10 8 6 2 4 3 2 4 5 10 9
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 4200*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
10 4
*x *y + 40800
x + y 10 8 6 2 4 3 2 4 5 8 15
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
- 17160
x + y 10 8 6 2 4 3 2 4 5 8 10
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 3360*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
8 5
*x *y + 7650
x + y 10 8 6 2 4 3 2 4 5 6 16
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 4680*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
6 11
*x *y + 1680
x + y 10 8 6 2 4 3 2 4 5 6 6
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 900*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
4 17
*x *y - 780
x + y 10 8 6 2 4 3 2 4 5 4 12
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 480*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
4 7
*x *y + 50
x + y 10 8 6 2 4 3 2 4 5 2 18
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 60*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 13
*x *y + 60
x + y 10 8 6 2 4 3 2 4 5 2 8
*kummerm(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 100*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
28
*x + 1300
x + y 10 8 6 2 4 3 2 4 5 26
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 7800*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
24 2
*x *y + 28600
x + y 10 8 6 2 4 3 2 4 5 22 3
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
+ 71500
x + y 10 8 6 2 4 3 2 4 5 20 4
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
+ 128700
x + y 10 8 6 2 4 3 2 4 5 18 5
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 160*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
18
*x + 171600
x + y 10 8 6 2 4 3 2 4 5 16 6
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 1280*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
16
*x *y + 171600
x + y 10 8 6 2 4 3 2 4 5 14 7
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 4480*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
14 2
*x *y + 128700
x + y 10 8 6 2 4 3 2 4 5 12 8
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 8960*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
12 3
*x *y + 71500
x + y 10 8 6 2 4 3 2 4 5 10 9
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
- 11200
x + y 10 8 6 2 4 3 2 4 5 10 4
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
+ 28600
x + y 10 8 6 2 4 3 2 4 5 8 10
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 8960*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
8 5
*x *y + 7800
x + y 10 8 6 2 4 3 2 4 5 6 11
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 4480*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
6 6
*x *y + 1300
x + y 10 8 6 2 4 3 2 4 5 4 12
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 1280*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
4 7
*x *y + 100
x + y 10 8 6 2 4 3 2 4 5 2 13
*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 160*kummerm(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
2 8
*x *y + 100
x + y 10 8 6 2 4 3 2 4 5 18
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y 10 8 6 2 4 3 2 4 5 16
800*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
*y + 2800
x + y 10 8 6 2 4 3 2 4 5 14 2
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 5600*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
12 3
*x *y + 7000
x + y 10 8 6 2 4 3 2 4 5 10 4
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 5600*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
8 5
*x *y + 2800
x + y 10 8 6 2 4 3 2 4 5 6 6
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 800*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
4 7
*x *y + 100
x + y 10 8 6 2 4 3 2 4 5 2 8
*kummerm(0,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y )
2*x + 2*y x + y 10 x + y 8 x + y 6 2
/(x*(x + 2*x *x + 10*x *x *y + 20*x *x *y
x + y 4 3 x + y 2 4 x + y 5 20 18
+ 20*x *x *y + 10*x *x *y + 2*x *y + x + 10*x *y
16 2 14 3 12 4 10 5 8 6 6 7
+ 45*x *y + 120*x *y + 210*x *y + 252*x *y + 210*x *y + 120*x *y
4 8 2 9 10
+ 45*x *y + 10*x *y + y ))
df(kummerU(a,b,z),x);
2*x + 2*y
(x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
2*x + 2*y
*x + x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
2*x + 2*y
x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y +
x + y x + y 10 8 6 2 4 3 2 4 5
x *kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
11 x + y
*log(x)*x + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
9 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
7 2 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 3 x + y
*x *y + 5*x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
3 4 x + y
*x *y + x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
5 x + y
*x*y - x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
x + y
*x + x
x + y 10 8 6 2 4 3 2 4 5 11
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y x + y 10 8 6 2 4 3 2 4 5
x *kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
10 x + y
*x *y + 10*x
x + y 10 8 6 2 4 3 2 4 5 10
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x +
x + y
5*x
x + y 10 8 6 2 4 3 2 4 5 9
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
x + y
5*x
x + y 10 8 6 2 4 3 2 4 5 8 2
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 8
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y +
x + y
10*x
x + y 10 8 6 2 4 3 2 4 5 7 2
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 10*x
x + y 10 8 6 2 4 3 2 4 5 6 3
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 60*x
x + y 10 8 6 2 4 3 2 4 5 6 2
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 10*x
x + y 10 8 6 2 4 3 2 4 5 5 3
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 10*x
x + y 10 8 6 2 4 3 2 4 5 4 4
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 40*x
x + y 10 8 6 2 4 3 2 4 5 4 3
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 5*x
x + y 10 8 6 2 4 3 2 4 5 3 4
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 5*x
x + y 10 8 6 2 4 3 2 4 5 2 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ 10*x
x + y 10 8 6 2 4 3 2 4 5 2 4
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y
+ x
x + y 10 8 6 2 4 3 2 4 5 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x*y -
x + y x + y 10 8 6 2 4 3 2 4 5
x *kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
x + y
*x + x
x + y 10 8 6 2 4 3 2 4 5 6
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*y -
x + y x + y 10 8 6 2 4 3 2 4 5
x *kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
x + y
*y - x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*log(x)
x + y
*x - x
x + y 10 8 6 2 4 3 2 4 5
*kummeru(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x -
x + y x + y 10 8 6 2 4 3 2 4 5
x *kummeru(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
x + y 10 8 6 2 4 3 2 4 5
*y + 10*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
20
*x + 90
x + y 10 8 6 2 4 3 2 4 5 18
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 360*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
16 2
*x *y + 840
x + y 10 8 6 2 4 3 2 4 5 14 3
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
+ 1260*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
12 4
*x *y + 1260
x + y 10 8 6 2 4 3 2 4 5 10 5
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 10*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
10
*x + 840
x + y 10 8 6 2 4 3 2 4 5 8 6
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5 8
- 40*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
*y + 360
x + y 10 8 6 2 4 3 2 4 5 6 7
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5 6
- 60*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
2
*y + 90
x + y 10 8 6 2 4 3 2 4 5 4 8
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5 4
- 40*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
3
*y + 10
x + y 10 8 6 2 4 3 2 4 5 2 9
*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5 2
- 10*kummeru(2,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
4
*y - 10
x + y 10 8 6 2 4 3 2 4 5 10
*kummeru(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x -
x + y 10 8 6 2 4 3 2 4 5 8
40*kummeru(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5
- 60*kummeru(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )
6 2
*x *y - 40
x + y 10 8 6 2 4 3 2 4 5 4 3
*kummeru(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x *y
x + y 10 8 6 2 4 3 2 4 5 2
- 10*kummeru(1,3,x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y )*x
4 x + y 10 8 6 2 4 3 2 4 5
*y )/(x*(x + x + 5*x *y + 10*x *y + 10*x *y + 5*x *y + y ))
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
% 9. Bessel functions
% =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
% Lengthy test of the Bessel functions. This isn't even
% remotely exhaustive of the special cases -- though a
% real person with lots of time could no doubt come up
% with a better lot of tests than this automated rubbish.
% Again, compromise by only actually doing one in five or
% nine. If you want a really thorough test, you can
% easily change this to get it; but it'll take hours to
% run.
clear p, q;
hankel1(p,q);
hankel1(p,q)
r := df(ws,q);
hankel1(p - 1,q)*q - hankel1(p,q)*p
r := -------------------------------------
q
on complex;
r;
besselj(p - 1,q)*q - besselj(p,q)*p + i*bessely(p - 1,q)*q - i*bessely(p,q)*p
-------------------------------------------------------------------------------
q
p:=3/8;
3
p := ---
8
r;
- 5 3*pi - 5 3*pi
(8*i*besselj(------,q)*cos(------)*q + 8*besselj(------,q)*sin(------)*q
8 8 8 8
- 3 5
+ 3*i*besselj(------,q) + 8*i*besselj(---,q)*q
8 8
3 3*pi 3 3*pi
- 3*i*besselj(---,q)*cos(------) - 3*besselj(---,q)*sin(------))/(8
8 8 8 8
3*pi
*sin(------)*q)
8
q := pi;
q := pi
r;
- 5 3*pi - 5 3*pi
(8*i*besselj(------,pi)*cos(------)*pi + 8*besselj(------,pi)*sin(------)*pi
8 8 8 8
- 3 5
+ 3*i*besselj(------,pi) + 8*i*besselj(---,pi)*pi
8 8
3 3*pi 3 3*pi
- 3*i*besselj(---,pi)*cos(------) - 3*besselj(---,pi)*sin(------))/(8
8 8 8 8
3*pi
*sin(------)*pi)
8
on rounded;
*** Domain mode complex changed to complex-rounded
r;
- 0.119366207318921501826662822529
*hankel1(0.375,3.14159265358979323846264338328)
+ hankel1( - 0.625,3.14159265358979323846264338328)
off complex, rounded;
*** Domain mode complex-rounded changed to rounded
df(df(besselj(pp,qq)+rr * hankel1(pp*2,qq) *
bessely(pp-qq,qq),qq),qq);
2
(besselj(pp - 2,qq)*qq - 2*besselj(pp - 1,qq)*pp*qq + besselj(pp - 1,qq)*qq
2
+ besselj(pp,qq)*pp + besselj(pp,qq)*pp
2
+ bessely(pp - qq,qq)*hankel1(2*pp - 2,qq)*qq *rr
- 4*bessely(pp - qq,qq)*hankel1(2*pp - 1,qq)*pp*qq*rr
+ bessely(pp - qq,qq)*hankel1(2*pp - 1,qq)*qq*rr
2
+ 4*bessely(pp - qq,qq)*hankel1(2*pp,qq)*pp *rr
+ 2*bessely(pp - qq,qq)*hankel1(2*pp,qq)*pp*rr
2
+ df(bessely(pp - qq,qq),qq,2)*hankel1(2*pp,qq)*qq *rr
2
+ 2*df(bessely(pp - qq,qq),qq)*hankel1(2*pp - 1,qq)*qq *rr
2
- 4*df(bessely(pp - qq,qq),qq)*hankel1(2*pp,qq)*pp*qq*rr)/qq
% Possible values for real args
t!*bes!*vr := { 1, pi, -pi, 26 }$
% Possible values for real and imaginary parts of complex args
t!*bes!*vc := { 0, 3, -41/2 }$
array s!*bes(4)$
s!*bes(1) := "BesselJ"$
s!*bes(2) := "BesselY"$
s!*bes(3) := "BesselI"$
s!*bes(4) := "BesselK"$
pre := 16;
pre := 16
precision pre;
30
preord := 10**pre;
preord := 10000000000000000
t!*b!*c := 3;
t*b*c := 3
algebraic procedure do!*one!*bessel(s,n,z);
(if s = 1 then besselj(n,z)
else if s = 2 then bessely(n,z)
else if s = 3 then besseli(n,z)
else besselk(n,z));
do*one*bessel
algebraic procedure pr!*bessel(s,n,z,k);
<< if t!*b!*c = k
then
<< on rounded;
bes1 := do!*one!*bessel(s,n,z);
precision(pre+5);
bes2 := do!*one!*bessel(s,n,z);
if bes1 neq 0
then disc := floor abs(100*(bes2-bes1)*preord/bes1)
else disc := 0;
precision pre;
write s!*bes(s), "(", n, ",", z, ") = ", bes1;
if not numberp disc then
<< precom := !*complex;
on complex;
disc := disc;
if not precom then off complex >>;
if disc neq 0
then write " (discrepancy ", disc, "% of one s.f.)";
if numberp disc and disc > 200 then
<< write "***** WARNING Significant Inaccuracy.";
write " Lower precision result:";
write " ", bes1;
write " Higher precision result:";
precision(pre+5); write " ", bes2; precision pre >>;
off rounded;
t!*b!*c := 0 >>;
t!*b!*c := t!*b!*c + 1 >>;
pr*bessel
% About to begin Bessel test. We have a list of possible
% values, and we test every Bessel, with every value on the
% list as both order and argument. Every Bessel is computed
% twice, to different precisions (differing by 3), and any
% discrepancy is reported. The value reported is the diff-
% erence between the two computed values, expressed as a
% percentage of the unit of the least significant displayed
% digit. A discrepancy between 100% and 200% means that the
% last digit of the displayed value was found to differ at
% higher precision; values greater than 200% are cause for
% concern. An ideal discrepancy would be between about 1%
% and 20%; if the value is found to be zero, no discrepancy
% is reported.
off msg;
for s := 1:4 do
<< write(" ... Testing ", s!*bes(s), " for real domains ... ");
for each n in t!*bes!*vr do
for each z in t!*bes!*vr do
pr!*bessel(s, n, z, 5) >>;
... Testing BesselJ for real domains ...
BesselJ(1, - 3.141592653589793) = - 0.2846153431797528
BesselJ(3.141592653589793,26) = - 0.006989220174690161
(discrepancy 5% of one s.f.)
BesselJ(26,1) = 3.660826744416803e-35
... Testing BesselY for real domains ...
BesselY(1,3.141592653589793) = 0.358872916776719
BesselY(3.141592653589793, - 3.141592653589793) =
6.283185307179586
0.1545613960392598*( - 1.570796326794897) - 4.829362563540275
--------------------------------------------------------------------------------
3.141592653589793
( - 1.570796326794897)
BesselY( - 3.141592653589793,26) = - 0.1386083623177915
... Testing BesselI for real domains ...
BesselI(1,1) = 0.565159103992485
BesselI(3.141592653589793,3.141592653589793) = 1.011423335928613
BesselI( - 3.141592653589793, - 3.141592653589793) = - 0.8856101155917482
+ 0.4221616153281286*i
BesselI(26,26) = 68397.86776155122
... Testing BesselK for real domains ...
BesselK(3.141592653589793,1) = 9.025908765806763
BesselK( - 3.141592653589793,3.141592653589793) = 0.1107526602738113
(discrepancy 1% of one s.f.)
BesselK(26, - 3.141592653589793) = besselk(26, - 3.141592653589793)
on complex;
for s := 1:3 do
<< write (" ... Testing ", s!*bes(s), " for complex domains ... ");
for each nr in t!*bes!*vc do
for each ni in t!*bes!*vc do
for each zr in t!*bes!*vc do
for each zi in t!*bes!*vc do
pr!*bessel(s, nr+ni*i, zr+zi*i, 9) >>;
... Testing BesselJ for complex domains ...
BesselJ(0, - 20.5 + 3.0*i) = 1.05389016561334 + 1.410918160335249*i
BesselJ(3*i, - 20.5 + 3.0*i) = 0.01225787392170983 + 0.01066256817009466*i
BesselJ( - 20.5*i, - 20.5 + 3.0*i) = - 6.607837931625446e+38
+ 7.203284455482089e+38*i
BesselJ(3, - 20.5 + 3.0*i) = 1.568613483726435 - 0.7011991107137573*i
BesselJ(3 + 3*i, - 20.5 + 3.0*i) = 0.007904103001381543 - 0.006566520928092784*i
BesselJ(3.0 - 20.5*i, - 20.5 + 3.0*i) = - 7.069920310202644e+37
- 1.753271554229047e+37*i
BesselJ( - 20.5, - 20.5 + 3.0*i) = 0.1758742246345278 - 0.332739860634916*i
BesselJ( - 20.5 + 3.0*i, - 20.5 + 3.0*i) = 0.08815299110072903
- 0.1369698556512304*i
BesselJ( - 20.5 - 20.5*i, - 20.5 + 3.0*i) = - 5.364748129151297e+46
+ 2.608178375230083e+47*i
... Testing BesselY for complex domains ...
BesselY(0, - 20.5 + 3.0*i) = - 1.404746667469566 + 1.060048452645186*i
BesselY(3*i, - 20.5 + 3.0*i) = 0.4973091982659732 + 0.7985114801567726*i
BesselY( - 20.5*i, - 20.5 + 3.0*i) = - 7.203284455482089e+38
- 6.607837931625446e+38*i
(discrepancy 10% of one s.f.)
BesselY(3, - 20.5 + 3.0*i) = 0.6963128100601111 + 1.576222640523309*i
BesselY(3 + 3*i, - 20.5 + 3.0*i) = - 1.117333163968302 + 0.9789575771194796*i
BesselY(3.0 - 20.5*i, - 20.5 + 3.0*i) = 1.753271554229047e+37
- 7.069920310202644e+37*i
BesselY( - 20.5, - 20.5 + 3.0*i) = 0.2353954565395826 + 0.144691313932682*i
BesselY( - 20.5 + 3.0*i, - 20.5 + 3.0*i) = - 0.1527215881543493
+ 0.2371137974094512*i
BesselY( - 20.5 - 20.5*i, - 20.5 + 3.0*i) = - 2.608178375230083e+47
- 5.364748129151297e+46*i
... Testing BesselI for complex domains ...
BesselI(0, - 20.5 + 3.0*i) = - 6.891185608459107e+7 - 1.506065792318474e+7*i
BesselI(3*i, - 20.5 + 3.0*i) = - 6879.511500081906 - 1745.250262122384*i
BesselI( - 20.5*i, - 20.5 + 3.0*i) = 4.756052726395977e+40
- 1.836844915663626e+40*i
BesselI(3, - 20.5 + 3.0*i) = 5.56813934269752e+7 + 1.026348768124686e+7*i
BesselI(3 + 3*i, - 20.5 + 3.0*i) = 5917.619873410601 - 1360.623977225443*i
BesselI(3.0 - 20.5*i, - 20.5 + 3.0*i) = 1.125266362182384e+40
- 3.487618645712341e+39*i
BesselI( - 20.5, - 20.5 + 3.0*i) = - 0.0005935206748136158
- 0.001045512507192035*i
BesselI( - 20.5 + 3.0*i, - 20.5 + 3.0*i) = 0.0001367212493519011
- 0.00008899048591346719*i
BesselI( - 20.5 - 20.5*i, - 20.5 + 3.0*i) = - 4.705541285798881e+45
+ 6.968045629188714e+45*i
off complex;
on msg;
write (" ...");
...
write ("Bessel test complete.");
Bessel test complete.
end;
Time for test: 4029 ms, plus GC time: 281 ms