Nemtsovapq.txt

Nemtsovapq.txt is C++ code, related to the asymptotic expansion of function SuNem for positive parameters \(p\) and \(q\) and large negative values of the argument.

After calculation, "Power" had been replaced to "pow" in order to provide compatibility with the C++ compilers.

C++ code

int NEMTSOVM=18;
DB APQ[18][9];
DB NEMTSOVp;
DB NEMTSOVq;
int nemtsovapq(DB p, DB q)
{
NEMTSOVp=p;
NEMTSOVq=q;

APQ[0][0]=1.;

APQ[1][0]=(-1.*q)/p;

APQ[2][0]=0.;
APQ[2][1]=(0.25*(3.*pow(p,3) + 2.*pow(q,2)))/pow(p,2);

APQ[3][0]=q + (3.*pow(q,3))/pow(p,3);
APQ[3][1]=-1.5*q - (1.*pow(q,3))/pow(p,3);

APQ[4][0]=(0.125*(5.*pow(p,6) - 4.*pow(p,3)*pow(q,2) - 44.*pow(q,4)))/pow(p,4);
APQ[4][1]=(-0.125*pow(3.*pow(p,3) + 2.*pow(q,2),2))/pow(p,4);
APQ[4][2]=(0.09375*pow(3.*pow(p,3) + 2.*pow(q,2),2))/pow(p,4);

APQ[5][0]=(-0.08333333333333333*(39.*pow(p,6)*q + 104.*pow(p,3)*pow(q,3) + 4.*pow(q,5)))/pow(p,5);
APQ[5][1]=5.25*p*q + (14.*pow(q,3))/pow(p,2) + (7.*pow(q,5))/pow(p,5);
APQ[5][2]=(-0.25*q*pow(3.*pow(p,3) + 2.*pow(q,2),2))/pow(p,5);

APQ[6][0]=-1.59375*pow(p,3) + 2.0208333333333335*pow(q,2) + (25.333333333333332*pow(q,4))/pow(p,3) + (21.666666666666668*pow(q,6))/pow(p,6);
APQ[6][1]=(0.03125*(129.*pow(p,9) + 98.*pow(p,6)*pow(q,2) - 628.*pow(p,3)*pow(q,4) - 424.*pow(q,6)))/pow(p,6);
APQ[6][2]=(-0.125*pow(3.*pow(p,3) + 2.*pow(q,2),3))/pow(p,6);
APQ[6][3]=(0.0390625*pow(3.*pow(p,3) + 2.*pow(q,2),3))/pow(p,6);

APQ[7][0]=(0.004166666666666667*(2613.*pow(p,9)*q + 8606.*pow(p,6)*pow(q,3) - 32.*pow(p,3)*pow(q,5) - 9904.*pow(q,7)))/pow(p,7);
APQ[7][1]=(-0.5*(45.*pow(p,9)*q + 150.*pow(p,6)*pow(q,3) + 104.*pow(p,3)*pow(q,5) + 16.*pow(q,7)))/pow(p,7);
APQ[7][2]=(0.0625*q*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(27.*pow(p,3) + 46.*pow(q,2)))/pow(p,7);
APQ[7][3]=(-0.125*q*pow(3.*pow(p,3) + 2.*pow(q,2),3))/pow(p,7);

APQ[8][0]=(-0.00017361111111111112*(-27945.*pow(p,12) + 63024.*pow(p,9)*pow(q,2) + 741328.*pow(p,6)*pow(q,4) + 952224.*pow(p,3)*pow(q,6) + 37008.*pow(q,8)))/pow(p,8);
APQ[8][1]=(0.0026041666666666665*(-5535.*pow(p,12) - 1380.*pow(p,9)*pow(q,2) + 63916.*pow(p,6)*pow(q,4) + 92896.*pow(p,3)*pow(q,6) + 34208.*pow(q,8)))/pow(p,8);
APQ[8][2]=(0.01953125*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(89.*pow(p,6) + 44.*pow(p,3)*pow(q,2) - 284.*pow(q,4)))/pow(p,8);
APQ[8][3]=(-0.09244791666666667*pow(3.*pow(p,3) + 2.*pow(q,2),4))/pow(p,8);
APQ[8][4]=(0.01708984375*pow(3.*pow(p,3) + 2.*pow(q,2),4))/pow(p,8);

APQ[9][0]=(0.0000992063492063492*(-395109.*pow(p,12)*q - 1.479384e6*pow(p,9)*pow(q,3) + 197084.*pow(p,6)*pow(q,5) + 3.924576e6*pow(p,3)*pow(q,7) + 1.863648e6*pow(q,9)))/pow(p,9);
APQ[9][1]=(0.008333333333333333*(11889.*pow(p,12)*q + 47244.*pow(p,9)*pow(q,3) + 35476.*pow(p,6)*pow(q,5) - 22096.*pow(p,3)*pow(q,7) - 18848.*pow(q,9)))/pow(p,9);
APQ[9][2]=(-0.03125*q*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(321.*pow(p,6) + 832.*pow(p,3)*pow(q,2) + 220.*pow(q,4)))/pow(p,9);
APQ[9][3]=(0.14583333333333334*q*pow(3.*pow(p,3) + 2.*pow(q,2),3)*(9.*pow(p,3) + 14.*pow(q,2)))/pow(p,9);
APQ[9][4]=(-0.0625*q*pow(3.*pow(p,3) + 2.*pow(q,2),4))/pow(p,9);

APQ[10][0]=(0.0000496031746031746*(-317079.*pow(p,15) + 1.141368e6*pow(p,12)*pow(q,2) + 1.3035486e7*pow(p,9)*pow(q,4) + 2.154074e7*pow(p,6)*pow(q,6) + 1.728336e6*pow(p,3)*pow(q,8) - 6.964704e6*pow(q,10)))/pow(p,10);
APQ[10][1]=(-0.00013020833333333333*(-417555.*pow(p,15) + 247446.*pow(p,12)*pow(q,2) + 8.939976e6*pow(p,9)*pow(q,4) + 1.7083184e7*pow(p,6)*pow(q,6) + 8.930576e6*pow(p,3)*pow(q,8) + 906208.*pow(q,10)))/pow(p,10);
APQ[10][2]=(-0.0009765625*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(8205.*pow(p,9) + 790.*pow(p,6)*pow(q,2) - 70016.*pow(p,3)*pow(q,4) - 56992.*pow(q,6)))/pow(p,10);
APQ[10][3]=(0.0003255208333333333*pow(3.*pow(p,3) + 2.*pow(q,2),3)*(5283.*pow(p,6) + 3684.*pow(p,3)*pow(q,2) - 12212.*pow(q,4)))/pow(p,10);
APQ[10][4]=(-0.060546875*pow(3.*pow(p,3) + 2.*pow(q,2),5))/pow(p,10);
APQ[10][5]=(0.0076904296875*pow(3.*pow(p,3) + 2.*pow(q,2),5))/pow(p,10);

APQ[11][0]=(-8.267195767195768e-6*(-1.7834835e7*pow(p,15)*q - 7.3008222e7*pow(p,12)*pow(q,3) + 3.0118484e7*pow(p,9)*pow(q,5) + 3.60190216e8*pow(p,6)*pow(q,7) + 2.98650816e8*pow(p,3)*pow(q,9) + 1.1423616e7*pow(q,11)))/pow(p,11);
APQ[11][1]=(0.0000496031746031746*(-8.922663e6*pow(p,15)*q - 4.004469e7*pow(p,12)*pow(q,3) - 2.8714524e7*pow(p,9)*pow(q,5) + 6.0447704e7*pow(p,6)*pow(q,7) + 7.5662304e7*pow(p,3)*pow(q,9) + 2.1802944e7*pow(q,11)))/pow(p,11);
APQ[11][2]=(0.005208333333333333*q*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(10815.*pow(p,9) + 35626.*pow(p,6)*pow(q,2) + 13148.*pow(p,3)*pow(q,4) - 17528.*pow(q,6)))/pow(p,11);
APQ[11][3]=(-0.005208333333333333*q*pow(3.*pow(p,3) + 2.*pow(q,2),3)*(1983.*pow(p,6) + 5000.*pow(p,3)*pow(q,2) + 1492.*pow(q,4)))/pow(p,11);
APQ[11][4]=(0.0026041666666666665*q*pow(3.*pow(p,3) + 2.*pow(q,2),4)*(351.*pow(p,3) + 514.*pow(q,2)))/pow(p,11);
APQ[11][5]=(-0.03125*q*pow(3.*pow(p,3) + 2.*pow(q,2),5))/pow(p,11);

APQ[12][0]=(7.233796296296296e-7*(7.3691343e7*pow(p,18) - 3.94881996e8*pow(p,15)*pow(q,2) - 4.453983252e9*pow(p,12)*pow(q,4) - 8.6891488e9*pow(p,9)*pow(q,6) - 7.032344e8*pow(p,6)*pow(q,8) + 7.13340864e9*pow(p,3)*pow(q,10) + 2.34211776e9*pow(q,12)))/pow(p,12);
APQ[12][1]=(3.1001984126984127e-6*(-6.8160393e7*pow(p,18) + 1.20809412e8*pow(p,15)*pow(q,2) + 2.394735756e9*pow(p,12)*pow(q,4) + 5.442220032e9*pow(p,9)*pow(q,6) + 3.403845488e9*pow(p,6)*pow(q,8) - 3.35072064e8*pow(p,3)*pow(q,10) - 5.74833216e8*pow(q,12)))/pow(p,12);
APQ[12][2]=(0.000016276041666666666*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(2.269485e6*pow(p,12) - 1.364592e6*pow(p,9)*pow(q,2) - 3.7748624e7*pow(p,6)*pow(q,4) - 5.0972192e7*pow(p,3)*pow(q,6) - 8.403664e6*pow(q,8)))/pow(p,12);
APQ[12][3]=(-0.00008138020833333333*pow(3.*pow(p,3) + 2.*pow(q,2),3)*(121953.*pow(p,9) + 51926.*pow(p,6)*pow(q,2) - 828712.*pow(p,3)*pow(q,4) - 675760.*pow(q,6)))/pow(p,12);
APQ[12][4]=(0.0007120768229166666*pow(3.*pow(p,3) + 2.*pow(q,2),4)*(2043.*pow(p,6) + 1668.*pow(p,3)*pow(q,2) - 3668.*pow(q,4)))/pow(p,12);
APQ[12][5]=(-0.03714599609375*pow(3.*pow(p,3) + 2.*pow(q,2),6))/pow(p,12);
APQ[12][6]=(0.0035247802734375*pow(3.*pow(p,3) + 2.*pow(q,2),6))/pow(p,12);

APQ[13][0]=(3.0062530062530064e-7*(-1.896359886e9*pow(p,18)*q - 8.211727518e9*pow(p,15)*pow(q,3) + 7.330312149e9*pow(p,12)*pow(q,5) + 6.859224718e10*pow(p,9)*pow(q,7) + 7.693914458e10*pow(p,6)*pow(q,9) + 7.0978104e9*pow(p,3)*pow(q,11) - 1.01427504e10*pow(q,13)))/pow(p,13);
APQ[13][1]=(0.0000496031746031746*(4.0136247e7*pow(p,18)*q + 1.96336866e8*pow(p,15)*pow(q,3) + 1.10946972e8*pow(p,12)*pow(q,5) - 6.3236084e8*pow(p,9)*pow(q,7) - 9.821076e8*pow(p,6)*pow(q,9) - 4.2415296e8*pow(p,3)*pow(q,11) - 3.322656e7*pow(q,13)))/pow(p,13);
APQ[13][2]=(-0.00029761904761904765*q*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(1.027449e6*pow(p,12) + 3.965799e6*pow(p,9)*pow(q,2) + 1.467231e6*pow(p,6)*pow(q,4) - 5.599756e6*pow(p,3)*pow(q,6) - 3.032108e6*pow(q,8)))/pow(p,13);
APQ[13][3]=(0.0026041666666666665*q*pow(3.*pow(p,3) + 2.*pow(q,2),3)*(27579.*pow(p,9) + 90218.*pow(p,6)*pow(q,2) + 40772.*pow(p,3)*pow(q,4) - 32072.*pow(q,6)))/pow(p,13);
APQ[13][4]=(-0.0026041666666666665*q*pow(3.*pow(p,3) + 2.*pow(q,2),4)*(3501.*pow(p,6) + 8622.*pow(p,3)*pow(q,2) + 2752.*pow(q,4)))/pow(p,13);
APQ[13][5]=(0.0046875*q*pow(3.*pow(p,3) + 2.*pow(q,2),5)*(127.*pow(p,3) + 178.*pow(q,2)))/pow(p,13);
APQ[13][6]=(-0.015625*q*pow(3.*pow(p,3) + 2.*pow(q,2),6))/pow(p,13);

APQ[14][0]=(4.6972703222703225e-9*(-3.9679860033e10*pow(p,21) + 3.0029217681e11*pow(p,18)*pow(q,2) + 3.38342616324e12*pow(p,15)*pow(q,4) + 7.408413837904e12*pow(p,12)*pow(q,6) + 1.46146770288e11*pow(p,9)*pow(q,8) - 1.1559854405216e13*pow(p,6)*pow(q,10) - 7.117415870976e12*pow(p,3)*pow(q,12) - 2.64421433856e11*pow(q,14)))/pow(p,14);
APQ[14][1]=(2.5834986772486774e-8*(3.2386607379e10*pow(p,21) - 1.07957407482e11*pow(p,18)*pow(q,2) - 1.733355516588e12*pow(p,15)*pow(q,4) - 4.449730470744e12*pow(p,12)*pow(q,6) - 3.03916666448e12*pow(p,9)*pow(q,8) + 1.47128341088e12*pow(p,6)*pow(q,10) + 2.08135714752e12*pow(p,3)*pow(q,12) + 4.9524541824e11*pow(q,14)))/pow(p,14);
APQ[14][2]=(-7.750496031746032e-7*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(2.19169629e8*pow(p,15) - 3.5553309e8*pow(p,12)*pow(q,2) - 6.166193808e9*pow(p,9)*pow(q,4) - 1.0890851408e10*pow(p,6)*pow(q,6) - 3.048459168e9*pow(p,3)*pow(q,8) + 1.691731008e9*pow(q,10)))/pow(p,14);
APQ[14][3]=(1.3563368055555556e-6*pow(3.*pow(p,3) + 2.*pow(q,2),3)*(4.0479075e7*pow(p,12) - 5.749176e6*pow(p,9)*pow(q,2) - 5.62200632e8*pow(p,6)*pow(q,4) - 7.73179616e8*pow(p,3)*pow(q,6) - 1.49793232e8*pow(q,8)))/pow(p,14);
APQ[14][4]=(-5.086263020833333e-6*pow(3.*pow(p,3) + 2.*pow(q,2),4)*(2.014419e6*pow(p,9) + 1.311338e6*pow(p,6)*pow(q,2) - 1.1310016e7*pow(p,3)*pow(q,4) - 9.2984e6*pow(q,6)))/pow(p,14);
APQ[14][5]=(1.0172526041666667e-6*pow(3.*pow(p,3) + 2.*pow(q,2),5)*(1.093887e6*pow(p,6) + 978036.*pow(p,3)*pow(q,2) - 1.595908e6*pow(q,4)))/pow(p,14);
APQ[14][6]=(-0.021883138020833335*pow(3.*pow(p,3) + 2.*pow(q,2),7))/pow(p,14);
APQ[14][7]=(0.001636505126953125*pow(3.*pow(p,3) + 2.*pow(q,2),7))/pow(p,14);

APQ[15][0]=(1.4453139453139454e-9*(1.552277202213e12*pow(p,21)*q + 6.896620479678e12*pow(p,18)*pow(q,3) - 1.1419431315332e13*pow(p,15)*pow(q,5) - 9.2933135666136e13*pow(p,12)*pow(q,7) - 1.26765283992576e14*pow(p,9)*pow(q,9) - 1.6668516962048e13*pow(p,6)*pow(q,11) + 4.4046690143232e13*pow(p,3)*pow(q,13) + 1.1169030341632e13*pow(q,15)))/pow(p,15);
APQ[15][1]=(1.5031265031265032e-7*(-5.9690999871e10*pow(p,21)*q - 3.09427026462e11*pow(p,18)*pow(q,3) - 8.0737547043e10*pow(p,15)*pow(q,5) + 1.819467675906e12*pow(p,12)*pow(q,7) + 3.2708358359e12*pow(p,9)*pow(q,9) + 1.76025326572e12*pow(p,6)*pow(q,11) + 4.27892112e10*pow(p,3)*pow(q,13) - 1.310337408e11*pow(q,15)))/pow(p,15);
APQ[15][2]=(0.00001240079365079365*q*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(1.30639407e8*pow(p,15) + 5.63112732e8*pow(p,12)*pow(q,2) + 1.43085632e8*pow(p,9)*pow(q,4) - 1.638610392e9*pow(p,6)*pow(q,6) - 1.553645952e9*pow(p,3)*pow(q,8) - 1.89063552e8*pow(q,10)))/pow(p,15);
APQ[15][3]=(-0.00008680555555555556*q*pow(3.*pow(p,3) + 2.*pow(q,2),3)*(5.350851e6*pow(p,12) + 2.0750376e7*pow(p,9)*pow(q,2) + 1.0410204e7*pow(p,6)*pow(q,4) - 2.2619104e7*pow(p,3)*pow(q,6) - 1.3090592e7*pow(q,8)))/pow(p,15);
APQ[15][4]=(0.0003255208333333333*q*pow(3.*pow(p,3) + 2.*pow(q,2),4)*(235065.*pow(p,9) + 762998.*pow(p,6)*pow(q,2) + 387404.*pow(p,3)*pow(q,4) - 202488.*pow(q,6)))/pow(p,15);
APQ[15][5]=(-0.0005208333333333333*q*pow(3.*pow(p,3) + 2.*pow(q,2),5)*(13968.*pow(p,6) + 33723.*pow(p,3)*pow(q,2) + 11234.*pow(q,4)))/pow(p,15);
APQ[15][6]=(0.000390625*q*pow(3.*pow(p,3) + 2.*pow(q,2),6)*(949.*pow(p,3) + 1286.*pow(q,2)))/pow(p,15);
APQ[15][7]=(-0.0078125*q*pow(3.*pow(p,3) + 2.*pow(q,2),7))/pow(p,15);

APQ[16][0]=(6.45229439872297e-12*(1.03363573937481e14*pow(p,24) - 1.059403995089136e15*pow(p,21)*pow(q,2) - 1.200470977464264e16*pow(p,18)*pow(q,4) - 2.845633966726925e16*pow(p,15)*pow(q,6) + 3.633958789784512e15*pow(p,12)*pow(q,8) + 7.547842672660685e16*pow(p,9)*pow(q,10) + 6.526436741685504e16*pow(p,6)*pow(q,12) + 5.97487843266048e15*pow(p,3)*pow(q,14) - 4.297771790862592e15*pow(q,16)))/pow(p,16);
APQ[16][1]=(-4.6972703222703225e-9*(7.13587936248e11*pow(p,24) - 3.793210310007e12*pow(p,21)*pow(q,2) - 5.5209684415527e13*pow(p,18)*pow(q,4) - 1.54964083625592e14*pow(p,15)*pow(q,6) - 1.06753897521072e14*pow(p,12)*pow(q,8) + 1.2437493276664e14*pow(p,9)*pow(q,10) + 1.9203309181448e14*pow(p,6)*pow(q,12) + 7.188859917504e13*pow(p,3)*pow(q,14) + 4.70701055424e12*pow(q,16)))/pow(p,16);
APQ[16][2]=(3.390842013888889e-7*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(2.293658127e9*pow(p,18) - 6.886986636e9*pow(p,15)*pow(q,2) - 1.00381971732e11*pow(p,12)*pow(q,4) - 2.1104075392e11*pow(p,9)*pow(q,6) - 7.366310032e10*pow(p,6)*pow(q,8) + 9.4438672512e10*pow(p,3)*pow(q,10) + 3.9241486464e10*pow(q,12)))/pow(p,16);
APQ[16][3]=(-9.68812003968254e-8*pow(3.*pow(p,3) + 2.*pow(q,2),3)*(3.041756865e9*pow(p,15) - 3.109356546e9*pow(p,12)*pow(q,2) - 7.3548639816e10*pow(p,9)*pow(q,4) - 1.33016094864e11*pow(p,6)*pow(q,6) - 4.4454764176e10*pow(p,3)*pow(q,8) + 1.4820204832e10*pow(q,10)))/pow(p,16);
APQ[16][4]=(3.814697265625e-6*pow(3.*pow(p,3) + 2.*pow(q,2),4)*(1.7521863e7*pow(p,12) + 3.392176e6*pow(p,9)*pow(q,2) - 2.08271792e8*pow(p,6)*pow(q,4) - 2.9140336e8*pow(p,3)*pow(q,6) - 6.2328944e7*pow(q,8)))/pow(p,16);
APQ[16][5]=(-2.5431315104166666e-7*pow(3.*pow(p,3) + 2.*pow(q,2),5)*(3.6779607e7*pow(p,9) + 2.9913834e7*pow(p,6)*pow(q,2) - 1.75313544e8*pow(p,3)*pow(q,4) - 1.45859632e8*pow(q,6)))/pow(p,16);
APQ[16][6]=(8.0108642578125e-6*pow(3.*pow(p,3) + 2.*pow(q,2),6)*(99109.*pow(p,6) + 94012.*pow(p,3)*pow(q,2) - 121196.*pow(q,4)))/pow(p,16);
APQ[16][7]=(-0.01254136221749442*pow(3.*pow(p,3) + 2.*pow(q,2),8))/pow(p,16);
APQ[16][8]=(0.0007671117782592773*pow(3.*pow(p,3) + 2.*pow(q,2),8))/pow(p,16);

APQ[17][0]=(-3.4412236793189175e-11*(2.59576629518043e14*pow(p,24)*q + 1.147050922259076e15*pow(p,21)*pow(q,3) - 3.246345313354224e15*pow(p,18)*pow(q,5) - 2.4292905239971984e16*pow(p,15)*pow(q,7) - 3.797143695129176e16*pow(p,12)*pow(q,9) - 5.387799652978176e15*pow(p,9)*pow(q,11) + 2.513458232213197e16*pow(p,6)*pow(q,13) + 1.2472178662207488e16*pow(p,3)*pow(q,15) + 4.48343243565056e14*pow(q,17)))/pow(p,17);
APQ[17][1]=(5.7812557812557815e-9*(6.984780601608e12*pow(p,24)*q + 3.7414035872124e13*pow(p,21)*pow(q,3) - 9.271185963951e12*pow(p,18)*pow(q,5) - 3.70498332766288e14*pow(p,15)*pow(q,7) - 7.41030880483656e14*pow(p,12)*pow(q,9) - 4.57227727967776e14*pow(p,9)*pow(q,11) + 5.136767236008e13*pow(p,6)*pow(q,13) + 1.2559826771136e14*pow(p,3)*pow(q,15) + 2.5744937944064e13*pow(q,17)))/pow(p,17);
APQ[17][2]=(-3.757816257816258e-8*q*pow(3.*pow(p,3) + 2.*pow(q,2),2)*(2.23842506121e11*pow(p,18) + 1.040873210778e12*pow(p,15)*pow(q,2) - 7.392887904e9*pow(p,12)*pow(q,4) - 5.43977603636e12*pow(p,9)*pow(q,6) - 6.91477030408e12*pow(p,6)*pow(q,8) - 1.563648912e12*pow(p,3)*pow(q,10) + 4.6174399104e11*pow(q,12)))/pow(p,17);
APQ[17][3]=(2.066798941798942e-6*q*pow(3.*pow(p,3) + 2.*pow(q,2),3)*(1.382218869e9*pow(p,15) + 6.036911286e9*pow(p,12)*pow(q,2) + 2.6720437e9*pow(p,9)*pow(q,4) - 1.409665812e10*pow(p,6)*pow(q,6) - 1.420387728e10*pow(p,3)*pow(q,8) - 2.06231328e9*pow(q,10)))/pow(p,17);
APQ[17][4]=(-0.000010850694444444445*q*pow(3.*pow(p,3) + 2.*pow(q,2),4)*(5.3384733e7*pow(p,12) + 2.07389868e8*pow(p,9)*pow(q,2) + 1.23604752e8*pow(p,6)*pow(q,4) - 1.78442672e8*pow(p,3)*pow(q,6) - 1.10799376e8*pow(q,8)))/pow(p,17);
APQ[17][5]=(0.00026041666666666666*q*pow(3.*pow(p,3) + 2.*pow(q,2),5)*(276969.*pow(p,9) + 892103.*pow(p,6)*pow(q,2) + 488552.*pow(p,3)*pow(q,4) - 180020.*pow(q,6)))/pow(p,17);
APQ[17][6]=(-0.00002170138888888889*q*pow(3.*pow(p,3) + 2.*pow(q,2),6)*(249111.*pow(p,6) + 591372.*pow(p,3)*pow(q,2) + 202892.*pow(q,4)))/pow(p,17);
APQ[17][7]=(0.00011160714285714285*q*pow(3.*pow(p,3) + 2.*pow(q,2),7)*(2003.*pow(p,3) + 2642.*pow(q,2)))/pow(p,17);
APQ[17][8]=(-0.00390625*q*pow(3.*pow(p,3) + 2.*pow(q,2),8))/pow(p,17);

return 18;}