//=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*
//   Symplectic法の展開係数リスト
//=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*=*

//---- 使用する展開係数のもののみのコメントを外す。
#define USE_10_2
//#define USE_10_4
//#define USE_10_6
//#define USE_10_8
//#define USE_11_4
//#define USE_11_6
//#define USE_11_8



#ifdef USE_10_2
// (10)式による2次近似展開の展開係数

const int symp_m = 2;
const int symp_r = 2;

const double symp_k[2] = {
	+0.50000000000000000000,
	+0.50000000000000000000
};

const double symp_u[2] = {
	+0.00000000000000000000,
	+1.00000000000000000000
};

#endif

#ifdef USE_10_4
// (10)式による4次近似展開の展開係数

const int symp_m = 4;
const int symp_r = 4;

const double symp_k[4] = {
	+0.67560359597982888591,
	-0.17560359597982888591,
	-0.17560359597982888591,
	+0.67560359597982888591
};

const double symp_u[4] = {
	+0.00000000000000000000,
	+1.35120719195965777182,
	-1.70241438391931554364,
	+1.35120719195965777182
};

#endif

#ifdef USE_10_6
// (10)式による6次近似展開の展開係数

const int symp_m = 6;
const int symp_r = 10;

const double symp_k[10] = {
	+0.79361246386112149764,
	-0.20627658481643984523,
	-0.20627658481643984523,
	-0.11800886788129261173,
	+0.23694957365305080454,
	+0.23694957365305080454,
	-0.11800886788129261173,
	-0.20627658481643984523,
	-0.20627658481643984523,
	+0.79361246386112149764
};

const double symp_u[10] = {
	+0.00000000000000000000,
	+1.58722492772224299529,
	-1.99977809735512268574,
	+1.58722492772224299529,
	-1.82324266348482821876,
	+2.29714181079092982785,
	-1.82324266348482821876,
	+1.58722492772224299529,
	-1.99977809735512268574,
	+1.58722492772224299529
};

#endif

#ifdef USE_10_8
// (10)式による8次近似展開の展開係数

const int symp_m = 8;
const int symp_r = 28;

const double symp_k[28] = {
	+0.88581669259776807390,
	-0.23024240475441615228,
	-0.23024240475441615228,
	-0.13171948501820229804,
	+0.26447907159195915838,
	+0.26447907159195915838,
	-0.13171948501820229804,
	-0.23024240475441615228,
	-0.23024240475441615228,
	-0.09220422873664646524,
	+0.25420822469239234831,
	+0.25420822469239234831,
	+0.14543010215511198435,
	-0.29200856953086762324,
	-0.29200856953086762324,
	+0.14543010215511198435,
	+0.25420822469239234831,
	+0.25420822469239234831,
	-0.09220422873664646524,
	-0.23024240475441615228,
	-0.23024240475441615228,
	-0.13171948501820229804,
	+0.26447907159195915838,
	+0.26447907159195915838,
	-0.13171948501820229804,
	-0.23024240475441615228,
	-0.23024240475441615228,
	+0.88581669259776807390
};

const double symp_u[28] = {
	+0.00000000000000000000,
	+1.77163338519553614780,
	-2.23211819470436845236,
	+1.77163338519553614780,
	-2.03507235523194074389,
	+2.56403049841585906066,
	-2.03507235523194074389,
	+1.77163338519553614780,
	-2.23211819470436845236,
	+1.77163338519553614780,
	-1.95604184266882907828,
	+2.46445829205361377490,
	-1.95604184266882907828,
	+2.24690204697905304698,
	-2.83091918604078829347,
	+2.24690204697905304698,
	-1.95604184266882907828,
	+2.46445829205361377490,
	-1.95604184266882907828,
	+1.77163338519553614780,
	-2.23211819470436845236,
	+1.77163338519553614780,
	-2.03507235523194074389,
	+2.56403049841585906066,
	-2.03507235523194074389,
	+1.77163338519553614780,
	-2.23211819470436845236,
	+1.77163338519553614780
};


#endif

#ifdef USE_11_4
// (11)式による4次近似展開の展開係数

const int symp_m = 4;
const int symp_r = 6;

const double symp_k[6] = {
	+0.20724538589718785597,
	+0.41449077179437571194,
	-0.12173615769156356792,
	-0.12173615769156356792,
	+0.41449077179437571194,
	+0.20724538589718785597
};

const double symp_u[6] = {
	+0.00000000000000000000,
	+0.41449077179437571194,
	+0.41449077179437571194,
	-0.65796308717750284778,
	+0.41449077179437571194,
	+0.41449077179437571194
};


#endif

#ifdef USE_11_6
// (11)式による6次近似展開の展開係数

const int symp_m = 6;
const int symp_r = 26;

const double symp_k[26] = {
	+0.07731617143363592182,
	+0.15463234286727184363,
	-0.04541560043427139703,
	-0.04541560043427139703,
	+0.15463234286727184363,
	+0.15463234286727184363,
	+0.15463234286727184363,
	-0.04541560043427139703,
	-0.04541560043427139703,
	+0.15463234286727184363,
	-0.02470312840371993723,
	-0.20403859967471171810,
	+0.05992624404552196471,
	+0.05992624404552196471,
	-0.20403859967471171810,
	-0.02470312840371993723,
	+0.15463234286727184363,
	-0.04541560043427139703,
	-0.04541560043427139703,
	+0.15463234286727184363,
	+0.15463234286727184363,
	+0.15463234286727184363,
	-0.04541560043427139703,
	-0.04541560043427139703,
	+0.15463234286727184363,
	+0.07731617143363592182
};

const double symp_u[26] = {
	+0.00000000000000000000,
	+0.15463234286727184363,
	+0.15463234286727184363,
	-0.24546354373581463770,
	+0.15463234286727184363,
	+0.15463234286727184363,
	+0.15463234286727184363,
	+0.15463234286727184363,
	-0.24546354373581463770,
	+0.15463234286727184363,
	+0.15463234286727184363,
	-0.20403859967471171810,
	-0.20403859967471171810,
	+0.32389108776575564752,
	-0.20403859967471171810,
	-0.20403859967471171810,
	+0.15463234286727184363,
	+0.15463234286727184363,
	-0.24546354373581463770,
	+0.15463234286727184363,
	+0.15463234286727184363,
	+0.15463234286727184363,
	+0.15463234286727184363,
	-0.24546354373581463770,
	+0.15463234286727184363,
	+0.15463234286727184363
};


#endif

#ifdef USE_11_8
// (11)式による8次近似展開の展開係数

const int symp_m = 8;
const int symp_r = 126;

const double symp_k[126] = {
	+0.02780170839404786168,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	-0.00888286576489946342,
	-0.07336914831789465019,
	+0.02154855745197108152,
	+0.02154855745197108152,
	-0.07336914831789465019,
	-0.00888286576489946342,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	-0.00888286576489946342,
	-0.07336914831789465019,
	+0.02154855745197108152,
	+0.02154855745197108152,
	-0.07336914831789465019,
	-0.00888286576489946342,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	-0.00608895374850765628,
	-0.06778132428511103591,
	+0.01990741059443593836,
	+0.01990741059443593836,
	-0.06778132428511103591,
	-0.06778132428511103591,
	-0.06778132428511103591,
	+0.01990741059443593836,
	+0.01990741059443593836,
	-0.06778132428511103591,
	+0.01082833465587792338,
	+0.08943799359686688266,
	-0.02626798576236234750,
	-0.02626798576236234750,
	+0.08943799359686688266,
	+0.01082833465587792338,
	-0.06778132428511103591,
	+0.01990741059443593836,
	+0.01990741059443593836,
	-0.06778132428511103591,
	-0.06778132428511103591,
	-0.06778132428511103591,
	+0.01990741059443593836,
	+0.01990741059443593836,
	-0.06778132428511103591,
	-0.00608895374850765628,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	-0.00888286576489946342,
	-0.07336914831789465019,
	+0.02154855745197108152,
	+0.02154855745197108152,
	-0.07336914831789465019,
	-0.00888286576489946342,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	-0.00888286576489946342,
	-0.07336914831789465019,
	+0.02154855745197108152,
	+0.02154855745197108152,
	-0.07336914831789465019,
	-0.00888286576489946342,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.01633075275717682864,
	-0.01633075275717682864,
	+0.05560341678809572336,
	+0.02780170839404786168
};

const double symp_u[126] = {
	+0.00000000000000000000,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.07336914831789465019,
	-0.07336914831789465019,
	+0.11646626322183681324,
	-0.07336914831789465019,
	-0.07336914831789465019,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.07336914831789465019,
	-0.07336914831789465019,
	+0.11646626322183681324,
	-0.07336914831789465019,
	-0.07336914831789465019,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.06778132428511103591,
	-0.06778132428511103591,
	+0.10759614547398291262,
	-0.06778132428511103591,
	-0.06778132428511103591,
	-0.06778132428511103591,
	-0.06778132428511103591,
	+0.10759614547398291262,
	-0.06778132428511103591,
	-0.06778132428511103591,
	+0.08943799359686688266,
	+0.08943799359686688266,
	-0.14197396512159157766,
	+0.08943799359686688266,
	+0.08943799359686688266,
	-0.06778132428511103591,
	-0.06778132428511103591,
	+0.10759614547398291262,
	-0.06778132428511103591,
	-0.06778132428511103591,
	-0.06778132428511103591,
	-0.06778132428511103591,
	+0.10759614547398291262,
	-0.06778132428511103591,
	-0.06778132428511103591,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.07336914831789465019,
	-0.07336914831789465019,
	+0.11646626322183681324,
	-0.07336914831789465019,
	-0.07336914831789465019,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.07336914831789465019,
	-0.07336914831789465019,
	+0.11646626322183681324,
	-0.07336914831789465019,
	-0.07336914831789465019,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	+0.05560341678809572336,
	-0.08826492230244938064,
	+0.05560341678809572336,
	+0.05560341678809572336
};
#endif