In[]:=
$Assumptions=0≤≤a&&0≤≤a&&>&&0≤≤b&&0≤≤b&&>&&a>0&&b>0&&a>b;
x
1
x
2
x
2
x
1
y
1
y
2
y
2
y
1
In[]:=
K1=Log[+]
1
2
a
2
b
b
∫
0
b
∫
0
a
∫
0
a
∫
0
1
2
2
(-)
x
2
x
1
2
(-)
y
2
y
1
x
1
x
2
y
1
y
2
1
12
2
a
2
b
2
b
3
b
a
b
2
a
2
b
b
a
4
a
4
b
4
a
2
a
2
b
4
b
2
a
2
b
In[]:=
a=5;b=2;
In[]:=
K1//N
0.448187
In[]:=
a=.;b=.;
In[]:=
NU=Numerator[K1]//Expand
Out[]=
-25+8aπ-8aArcTan+8bArcTan-16aArcTan+2Log[a]+2Log[b]-Log[+]+6Log[+]-Log[+]
2
a
2
b
3
b
3
b
a
b
3
a
b
a
3
b
b
a
4
a
4
b
4
a
2
a
2
b
2
a
2
b
2
a
2
b
4
b
2
a
2
b
In[]:=
NU=NU/.2Log[a]-Log[+]->-Log1+;
4
a
4
a
2
a
2
b
4
a
2
b
2
a
In[]:=
NU=NU/.2Log[b]-Log[+]->-Log1+;
4
b
4
b
2
a
2
b
4
b
2
a
2
b
In[]:=
K1=
NU
Denominator[K1]
Out[]=
1
12
2
a
2
b
2
a
2
b
3
b
3
b
a
b
3
a
b
a
3
b
b
a
4
b
2
a
2
b
2
a
2
b
2
a
2
b
4
a
2
b
2
a
In[]:=
a=5;b=2;
In[]:=
K1//N
Out[]=
0.448187
In[]:=
a=.;b=.;
In[]:=
Limit[K1,b0]//Expand
Out[]=
-+Log[a]
3
2
In[]:=
K=K1/.ba//FullSimplify
Out[]=
1
12
In[]:=
K//N
Out[]=
-0.805087+Log[a]