> restart:
alpha:=0.0;
s:=3-alpha;
assume(ecc in RealRange(0,1));
#ecc:=0.87;
eps:=sqrt(1-ecc^2);
prefaktor:=-2*k/(s-1)*(a/RA)^(s)*(sqrt(1-ecc^2)/ecc) ;
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen1:=convert( (2*s+1)*ecc*eps^(2*s-1)*int((1+ecc*cos(phi))^(-s-1),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen2:=convert( (s+1)*eps^(2*s+1)*int(cos(phi)*(1+ecc*cos(phi))^(-s-2),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
rozdil:=convert( (clen1+clen2) ,EllipticK);
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX;
 

0. 

3. 

(1-ecc^2)^(1/2) 

-1.000000000*k*(a/RA)^3.*(1-ecc^2)^(1/2)/ecc 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

7.*ecc*(1-ecc^2)^2.500000000*int(1/(1+ecc*cos(phi))^4., phi = 0 .. Pi) 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

4.*(1-ecc^2)^3.500000000*int(cos(phi)/(1+ecc*cos(phi))^5., phi = 0 .. Pi) 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

7.*ecc*(1-ecc^2)^2.500000000*int(1/(1+ecc*cos(phi))^4., phi = 0 .. Pi)+4.*(1-ecc^2)^3.500000000*int(cos(phi)/(1+ecc*cos(phi))^5., phi = 0 .. Pi)
7.*ecc*(1-ecc^2)^2.500000000*int(1/(1+ecc*cos(phi))^4., phi = 0 .. Pi)+4.*(1-ecc^2)^3.500000000*int(cos(phi)/(1+ecc*cos(phi))^5., phi = 0 .. Pi) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX		 (1)
 

>
 

> restart:
alpha:=0.5;
s:=3-alpha;
assume(ecc in RealRange(0,1));
#ecc:=0.87;
eps:=sqrt(1-ecc^2);
prefaktor:=-2*k/(s-1)*(a/RA)^(s)*(sqrt(1-ecc^2)/ecc) ;
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen1:=convert( (2*s+1)*ecc*eps^(2*s-1)*int((1+ecc*cos(phi))^(-s-1),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen2:=convert( (s+1)*eps^(2*s+1)*int(cos(phi)*(1+ecc*cos(phi))^(-s-2),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
rozdil:=convert( (clen1+clen2) ,EllipticK);
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX;
 

.5 

2.5 

(1-ecc^2)^(1/2) 

-1.333333333*k*(a/RA)^2.5*(1-ecc^2)^(1/2)/ecc 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562... 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

3.5*(1-ecc^2)^3.000000000*((-0.5714285714e-1*(5.+27.*ecc^2)/(ecc*(1.+ecc)^3*(-1.+ecc)^3)+0.5714285714e-1*(21.*ecc^4+102.*ecc^2+5.)/(ecc*(1.+ecc)^4*(-1.+ecc)^3)-.1142857143*(21.*ecc^4+102.*ecc^2+5.)/((...
3.5*(1-ecc^2)^3.000000000*((-0.5714285714e-1*(5.+27.*ecc^2)/(ecc*(1.+ecc)^3*(-1.+ecc)^3)+0.5714285714e-1*(21.*ecc^4+102.*ecc^2+5.)/(ecc*(1.+ecc)^4*(-1.+ecc)^3)-.1142857143*(21.*ecc^4+102.*ecc^2+5.)/((...
3.5*(1-ecc^2)^3.000000000*((-0.5714285714e-1*(5.+27.*ecc^2)/(ecc*(1.+ecc)^3*(-1.+ecc)^3)+0.5714285714e-1*(21.*ecc^4+102.*ecc^2+5.)/(ecc*(1.+ecc)^4*(-1.+ecc)^3)-.1142857143*(21.*ecc^4+102.*ecc^2+5.)/((...
3.5*(1-ecc^2)^3.000000000*((-0.5714285714e-1*(5.+27.*ecc^2)/(ecc*(1.+ecc)^3*(-1.+ecc)^3)+0.5714285714e-1*(21.*ecc^4+102.*ecc^2+5.)/(ecc*(1.+ecc)^4*(-1.+ecc)^3)-.1142857143*(21.*ecc^4+102.*ecc^2+5.)/((...
3.5*(1-ecc^2)^3.000000000*((-0.5714285714e-1*(5.+27.*ecc^2)/(ecc*(1.+ecc)^3*(-1.+ecc)^3)+0.5714285714e-1*(21.*ecc^4+102.*ecc^2+5.)/(ecc*(1.+ecc)^4*(-1.+ecc)^3)-.1142857143*(21.*ecc^4+102.*ecc^2+5.)/((...
3.5*(1-ecc^2)^3.000000000*((-0.5714285714e-1*(5.+27.*ecc^2)/(ecc*(1.+ecc)^3*(-1.+ecc)^3)+0.5714285714e-1*(21.*ecc^4+102.*ecc^2+5.)/(ecc*(1.+ecc)^4*(-1.+ecc)^3)-.1142857143*(21.*ecc^4+102.*ecc^2+5.)/((...
3.5*(1-ecc^2)^3.000000000*((-0.5714285714e-1*(5.+27.*ecc^2)/(ecc*(1.+ecc)^3*(-1.+ecc)^3)+0.5714285714e-1*(21.*ecc^4+102.*ecc^2+5.)/(ecc*(1.+ecc)^4*(-1.+ecc)^3)-.1142857143*(21.*ecc^4+102.*ecc^2+5.)/((...
3.5*(1-ecc^2)^3.000000000*((-0.5714285714e-1*(5.+27.*ecc^2)/(ecc*(1.+ecc)^3*(-1.+ecc)^3)+0.5714285714e-1*(21.*ecc^4+102.*ecc^2+5.)/(ecc*(1.+ecc)^4*(-1.+ecc)^3)-.1142857143*(21.*ecc^4+102.*ecc^2+5.)/((...
3.5*(1-ecc^2)^3.000000000*((-0.5714285714e-1*(5.+27.*ecc^2)/(ecc*(1.+ecc)^3*(-1.+ecc)^3)+0.5714285714e-1*(21.*ecc^4+102.*ecc^2+5.)/(ecc*(1.+ecc)^4*(-1.+ecc)^3)-.1142857143*(21.*ecc^4+102.*ecc^2+5.)/((... 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562...
6.0*ecc*(1-ecc^2)^2.000000000*((-1.066666667/((1.+ecc)^2*(-1.+ecc)^2)+.1333333333*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^2)-.2666666667*ecc*(23.+9.*ecc^2)/((1.+ecc)^3*(-1.+ecc)^3))*EllipticK(1.414213562... 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX		 (2)
 

>
 

> restart:
alpha:=1.0;
s:=3-alpha;
assume(ecc in RealRange(0,1));
#ecc:=0.87;
eps:=sqrt(1-ecc^2);
prefaktor:=-2*k/(s-1)*(a/RA)^(s)*(sqrt(1-ecc^2)/ecc) ;
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen1:=convert( (2*s+1)*ecc*eps^(2*s-1)*int((1+ecc*cos(phi))^(-s-1),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen2:=convert( (s+1)*eps^(2*s+1)*int(cos(phi)*(1+ecc*cos(phi))^(-s-2),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
rozdil:=convert( (clen1+clen2) ,EllipticK);
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX;
 

1.0 

2.0 

(1-ecc^2)^(1/2) 

-2.000000000*k*(a/RA)^2.0*(1-ecc^2)^(1/2)/ecc 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

5.0*ecc*(1-ecc^2)^1.500000000*int(1/(1+ecc*cos(phi))^3.0, phi = 0 .. Pi) 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

3.0*(1-ecc^2)^2.500000000*int(cos(phi)/(1+ecc*cos(phi))^4.0, phi = 0 .. Pi) 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

5.0*ecc*(1-ecc^2)^1.500000000*int(1/(1+ecc*cos(phi))^3.0, phi = 0 .. Pi)+3.0*(1-ecc^2)^2.500000000*int(cos(phi)/(1+ecc*cos(phi))^4.0, phi = 0 .. Pi)
5.0*ecc*(1-ecc^2)^1.500000000*int(1/(1+ecc*cos(phi))^3.0, phi = 0 .. Pi)+3.0*(1-ecc^2)^2.500000000*int(cos(phi)/(1+ecc*cos(phi))^4.0, phi = 0 .. Pi) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX		 (3)
 

>
 

> restart:
alpha:=1.5;
s:=3-alpha;
assume(ecc in RealRange(0,1));
#ecc:=0.87;
eps:=sqrt(1-ecc^2);
prefaktor:=-2*k/(s-1)*(a/RA)^(s)*(sqrt(1-ecc^2)/ecc) ;
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen1:=convert( (2*s+1)*ecc*eps^(2*s-1)*int((1+ecc*cos(phi))^(-s-1),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen2:=convert( (s+1)*eps^(2*s+1)*int(cos(phi)*(1+ecc*cos(phi))^(-s-2),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
rozdil:=convert( (clen1+clen2) ,EllipticK);
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX;
 

1.5 

1.5 

(1-ecc^2)^(1/2) 

-4.000000000*k*(a/RA)^1.5*(1-ecc^2)^(1/2)/ecc 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2.... 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

2.5*(1-ecc^2)^2.000000000*((.1333333333*(3.+5.*ecc^2)/(ecc*(1.+ecc)^2*(ecc-1.)^2)-.1333333333*(29.*ecc^2+3.)/((1.+ecc)^3*ecc*(ecc-1.)^2)+.2666666667*(29.*ecc^2+3.)/((1.+ecc)^3*(ecc-1.)^3))*EllipticK(1...
2.5*(1-ecc^2)^2.000000000*((.1333333333*(3.+5.*ecc^2)/(ecc*(1.+ecc)^2*(ecc-1.)^2)-.1333333333*(29.*ecc^2+3.)/((1.+ecc)^3*ecc*(ecc-1.)^2)+.2666666667*(29.*ecc^2+3.)/((1.+ecc)^3*(ecc-1.)^3))*EllipticK(1...
2.5*(1-ecc^2)^2.000000000*((.1333333333*(3.+5.*ecc^2)/(ecc*(1.+ecc)^2*(ecc-1.)^2)-.1333333333*(29.*ecc^2+3.)/((1.+ecc)^3*ecc*(ecc-1.)^2)+.2666666667*(29.*ecc^2+3.)/((1.+ecc)^3*(ecc-1.)^3))*EllipticK(1...
2.5*(1-ecc^2)^2.000000000*((.1333333333*(3.+5.*ecc^2)/(ecc*(1.+ecc)^2*(ecc-1.)^2)-.1333333333*(29.*ecc^2+3.)/((1.+ecc)^3*ecc*(ecc-1.)^2)+.2666666667*(29.*ecc^2+3.)/((1.+ecc)^3*(ecc-1.)^3))*EllipticK(1...
2.5*(1-ecc^2)^2.000000000*((.1333333333*(3.+5.*ecc^2)/(ecc*(1.+ecc)^2*(ecc-1.)^2)-.1333333333*(29.*ecc^2+3.)/((1.+ecc)^3*ecc*(ecc-1.)^2)+.2666666667*(29.*ecc^2+3.)/((1.+ecc)^3*(ecc-1.)^3))*EllipticK(1...
2.5*(1-ecc^2)^2.000000000*((.1333333333*(3.+5.*ecc^2)/(ecc*(1.+ecc)^2*(ecc-1.)^2)-.1333333333*(29.*ecc^2+3.)/((1.+ecc)^3*ecc*(ecc-1.)^2)+.2666666667*(29.*ecc^2+3.)/((1.+ecc)^3*(ecc-1.)^3))*EllipticK(1... 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2....
4.0*ecc*(1-ecc^2)^1.000000000*((.6666666667/((1.+ecc)*(ecc-1.))-2.666666667/((1.+ecc)^2*(ecc-1.))+5.333333333*ecc/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)-2.... 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX		 (4)
 

>
 

> restart:
alpha:=2.5;
s:=3-alpha;
assume(ecc in RealRange(0,1));
#ecc:=0.87;
eps:=sqrt(1-ecc^2);
prefaktor:=-2*k/(s-1)*(a/RA)^(s)*(sqrt(1-ecc^2)/ecc) ;
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen1:=convert( (2*s+1)*ecc*eps^(2*s-1)*int((1+ecc*cos(phi))^(-s-1),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen2:=convert( (s+1)*eps^(2*s+1)*int(cos(phi)*(1+ecc*cos(phi))^(-s-2),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
rozdil:=convert( (clen1+clen2) ,EllipticK);
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX;
 

2.5 

.5 

(1-ecc^2)^(1/2) 

4.000000000*k*(a/RA)^.5*(1-ecc^2)^(1/2)/ecc 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414... 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

1.5*(1-ecc^2)^1.000000000*((-.6666666667/(ecc*(1.+ecc)*(ecc-1.))+.6666666667*(3.*ecc^2+1.)/(ecc*(1.+ecc)^2*(ecc-1.))-1.333333333*(3.*ecc^2+1.)/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+e...
1.5*(1-ecc^2)^1.000000000*((-.6666666667/(ecc*(1.+ecc)*(ecc-1.))+.6666666667*(3.*ecc^2+1.)/(ecc*(1.+ecc)^2*(ecc-1.))-1.333333333*(3.*ecc^2+1.)/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+e...
1.5*(1-ecc^2)^1.000000000*((-.6666666667/(ecc*(1.+ecc)*(ecc-1.))+.6666666667*(3.*ecc^2+1.)/(ecc*(1.+ecc)^2*(ecc-1.))-1.333333333*(3.*ecc^2+1.)/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+e...
1.5*(1-ecc^2)^1.000000000*((-.6666666667/(ecc*(1.+ecc)*(ecc-1.))+.6666666667*(3.*ecc^2+1.)/(ecc*(1.+ecc)^2*(ecc-1.))-1.333333333*(3.*ecc^2+1.)/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+e...
1.5*(1-ecc^2)^1.000000000*((-.6666666667/(ecc*(1.+ecc)*(ecc-1.))+.6666666667*(3.*ecc^2+1.)/(ecc*(1.+ecc)^2*(ecc-1.))-1.333333333*(3.*ecc^2+1.)/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+e...
1.5*(1-ecc^2)^1.000000000*((-.6666666667/(ecc*(1.+ecc)*(ecc-1.))+.6666666667*(3.*ecc^2+1.)/(ecc*(1.+ecc)^2*(ecc-1.))-1.333333333*(3.*ecc^2+1.)/((1.+ecc)^2*(ecc-1.)^2))*EllipticK(1.414213562*(ecc/(1.+e... 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414...
2.0*ecc*((2./(1.+ecc)-4.*ecc/((1.+ecc)*(ecc-1.)))*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/(1.+ecc)^(1/2)+2.*EllipticK(1.414213562*(ecc/(1.+ecc))^(1/2))/((1.+ecc)^(1/2)*(ecc-1.))-2.*EllipticE(1.414... 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX		 (5)
 

>
 

> restart:
alpha:=3.0;
s:=3-alpha;
assume(ecc in RealRange(0,1));
#ecc:=0.87;
eps:=sqrt(1-ecc^2);
prefaktor:=-2*k/(s-1)*(a/RA)^(s)*(sqrt(1-ecc^2)/ecc) ;
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen1:=convert( (2*s+1)*ecc*eps^(2*s-1)*int((1+ecc*cos(phi))^(-s-1),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
clen2:=convert( (s+1)*eps^(2*s+1)*int(cos(phi)*(1+ecc*cos(phi))^(-s-2),phi=0..Pi) ,EllipticK);
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB;
rozdil:=convert( (clen1+clen2) ,EllipticK);
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX;
 

3.0 

0. 

(1-ecc^2)^(1/2) 

2.*k*(1-ecc^2)^(1/2)/ecc 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_prefaktoruBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

3.141592654*ecc/((1-ecc^2)^.5000000000*(1.-1.*ecc^2)^(1/2)) 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu1BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

1.*(1-ecc^2)^.5000000000*int(cos(phi)/(1+ecc*cos(phi))^2., phi = 0 .. Pi) 

BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBpo_clenu2BBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 

3.141592654*ecc/((1-ecc^2)^.5000000000*(1.-1.*ecc^2)^(1/2))+1.*(1-ecc^2)^.5000000000*int(cos(phi)/(1+ecc*cos(phi))^2., phi = 0 .. Pi)
3.141592654*ecc/((1-ecc^2)^.5000000000*(1.-1.*ecc^2)^(1/2))+1.*(1-ecc^2)^.5000000000*int(cos(phi)/(1+ecc*cos(phi))^2., phi = 0 .. Pi) 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX		 (6)
 

>