#include <stdio.h>
#include <stdlib.h>
#include <math.h>

#define MAXSUM 10000

double li4half = 0.51747906167389938633;
double li5half = 0.50840057924226870746;
double li6half = 0.50409539780398855069;
double li7half = 0.50201456332470849457;
double s6 = 0.98744142640329971377;
double z2 = 1.64493406684822643647;
double z3 = 1.20205690315959428540;
double z4 = 1.08232323371113819152;
double z5 = 1.03692775514336992633;
double z6 = 1.01734306198444913971;
double z7 = 1.00834927738192282684;
double s7a = -0.95296007575629860341;
double s7b = 1.02912126296432453422;

double s1[MAXSUM+1];
double s11[MAXSUM+1];
double sm21[MAXSUM+1];
 
double sum1[MAXSUM+1];
double sum2[MAXSUM+1];
double sum3[MAXSUM+1];

int main()
{
	double ln2 = log(2.0);
	double ln22 = ln2*ln2;
	double ln23 = ln22*ln2;
	double ln24 = ln23*ln2;
	double ln25 = ln24*ln2;
	double ln26 = ln25*ln2;
	double ln27 = ln26*ln2;
	double a1 =
       - 1./2.*z2*ln24
       - 129./140.*z2*z2*z2
       + 7./6.*z3*ln23
       - 37./16.*z3*z3
       - 31./8.*z5*ln2
       + 8.*ln2*li5half
       + 4.*ln22*li4half
       + 1./9.*ln26
       + 8.*li6half
       + 2.*s6
			;
	double a2 =
       + 105./4.*z2*z3*ln22
       + 609./8.*z2*z5
       + 30.*z2*ln2*li4half
       + z2*ln25
       + 30.*z2*li5half
       + 717./20.*z2*z2*z3
       - 5.*z2*z2*ln23
       - 21./4.*z2*z2*z2*ln2
       - 35./16.*z3*ln24
       - 105./2.*z3*li4half
       - 75./2.*z3*z3*ln2
       - 30621./128.*z7
       + 30.*ln2*s6
       - 30.*s7a
       - 30.*s7b
		;
	double a3 =
       - 7./4.*z2*z3*ln22
       - 5673./448.*z2*z5
       - 5.*z2*ln2*li4half
       - 17./120.*z2*ln25
       - 5.*z2*li5half
       + 1517./1120.*z2*z2*z3
       + 5./6.*z2*z2*ln23
       - 1./84.*z2*z2*z2*ln2
       - 7./96.*z3*ln24
       - 3./4.*z3*li4half
       - 1563./448.*z3*z3*ln2
       - 93./32.*z5*ln22
       + 74415./1792.*z7
       - 18.*ln2*li6half
       - 43./14.*ln2*s6
       - 6.*ln22*li5half
       - ln23*li4half
       - 1./84.*ln27
       - 24.*li7half
       - 45./7.*s7a
       + 32./7.*s7b
		;
double a4 =
       + 25./8.*z2*z3*ln2
       + 31./4.*z2*z3
       + 6.*z2*ln2
       - 3.*z2*ln22
       - 11./6.*z2*ln23
       + 9./8.*z2*ln24
       + 19.*z2*li4half
       + 7./20.*z2*z2*ln2
       - 197./40.*z2*z2*ln22
       - 2853./280.*z2*z2*z2
       + 15./4.*z3*ln2
       + 39./8.*z3*ln22
       - 35./24.*z3*ln23
       - 27./4.*z3
       - 191./16.*z3*z3
       - 403./16.*z5*ln2
       - 1991./64.*z5
       + 16.*ln2*li4half
       + 24.*ln2*li5half
       + 3.*ln22*li4half
       + 31./60.*ln25
       - 1./72.*ln26
       + 18.*li5half
       + 44.*li6half
       + 36.*s6
		;

	printf("answer1 = %20.10f\n",a1);
	printf("answer2 = %20.10f\n",a2);
	printf("answer3 = %20.10f\n",a3);
	printf("answer4 = %20.10f\n",a4);
	return(0);
}