long double SqrtA01( long double x, long double c )
{
   return 1. / x * ( x * x + c ) * .5;
}
long double SqrtB01( long double x, long double c )
{
   return 2.0 * x * c / ( x * x + c );
}
long double SqrtA02( long double x, long double c )
{
   long double t1 = x * x;
   return x * ( t1 + 3.0 * c ) / ( 3.0 * t1 + c );
}
long double SqrtB02( long double x, long double c )
{
   long double t3 = x * x;
   return 1. / x * c * ( 3.0 * t3 + c ) / ( t3 + 3.0 * c );
}
long double SqrtA03( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t5 = c * c;
   return 1. / x * ( t3 + 6.0 * t2 * c + t5 ) / ( t2 + c ) * .25;
}
long double SqrtB03( long double x, long double c )
{
   long double t2 = x * x;
   long double t4 = t2 * t2;
   long double t6 = c * c;
   return 4.0 * x * c * ( t2 + c ) / ( t4 + 6.0 * t2 * c + t6 );

}
long double SqrtA04( long double x, long double c )
{
   long double t1 = x * x;
   long double t2 = t1 * t1;
   long double t3 = t1 * c;
   long double t4 = c * c;
   return x * ( t2 + 10.0 * t3 + 5.0 * t4 ) / ( 5.0 * t2 + 10.0 * t3 + t4 );
}
long double SqrtB04( long double x, long double c )
{
   long double t3 = x * x;
   long double t4 = t3 * t3;
   long double t5 = t3 * c;
   long double t6 = c * c;
   return 1. / x * c * ( 5.0 * t4 + 10.0 * t5 + t6 ) / ( t4 + 10.0 * t5 + 5.0 * t6 );

}
long double SqrtA05( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t6 = c * c;
   return 1. / x * ( t3 * t2 + 15.0 * t3 * c + 15.0 * t2 * t6 + t6 * c ) /
   ( 3.0 * t3 + 10.0 * t2 * c + 3.0 * t6 ) * .5;
}
long double SqrtB05( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t5 = c * c;
   return 2.0 * x * c * ( 3.0 * t3 + 10.0 * t2 * c + 3.0 * t5 )
   /
   ( t3 * t2 + 15.0 * t3 * c + 15.0 * t2 * t5 + t5 * c );
}
long double SqrtA06( long double x, long double c )
{
   long double t1 = x * x;
   long double t2 = t1 * t1;
   long double t3 = t2 * t1;
   long double t4 = t2 * c;
   long double t5 = c * c;
   long double t6 = t1 * t5;
   long double t7 = t5 * c;
   return x * ( t3 + 21.0 * t4 + 35.0 * t6 + 7.0 * t7 )
   /
   ( 7.0 * t3 + 35.0 * t4 + 21.0 * t6 + t7 );
}
long double SqrtB06( long double x, long double c )
{
   long double t3 = x * x;
   long double t4 = t3 * t3;
   long double t5 = t4 * t3;
   long double t6 = t4 * c;
   long double t7 = c * c;
   long double t8 = t3 * t7;
   long double t9 = t7 * c;
   return 1.
   /
   x * c * ( 7.0 * t5 + 35.0 * t6 + 21.0 * t8 + t9 )
   /
   ( t5 + 21.0 * t6 + 35.0 * t8 + 7.0 * t9 );
}
long double SqrtA07( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t5 = t3 * t2;
   long double t7 = c * c;
   long double t9 = t7 * c;
   long double t11 = t7 * t7;
   return 1.
   /
   x * ( t4 + 28.0 * t5 * c + 70.0 * t3 * t7 + 28.0 * t2 * t9 + t11 )
   /
   ( t5 + 7.0 * t3 * c + 7.0 * t2 * t7 + t9 ) * .125;
}
long double SqrtB07( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t2;
   long double t6 = c * c;
   long double t8 = t6 * c;
   long double t10 = t3 * t3;
   long double t14 = t6 * t6;
   return 8.0 * x * c * ( t4 + 7.0 * t3 * c + 7.0 * t2 * t6 + t8 )
   /
   ( t10 + 28.0 * t4 * c + 70.0 * t3 * t6 + 28.0 * t2 * t8 + t14 );
}
long double SqrtA08( long double x, long double c )
{
   long double t1 = x * x;
   long double t2 = t1 * t1;
   long double t3 = t2 * t2;
   long double t5 = t2 * t1 * c;
   long double t6 = c * c;
   long double t7 = t2 * t6;
   long double t9 = t1 * t6 * c;
   long double t10 = t6 * t6;
   return x * ( t3 + 36.0 * t5 + 126.0 * t7 + 84.0 * t9 + 9.0 * t10 )
   /
   ( 9.0 * t3 + 84.0 * t5 + 126.0 * t7 + 36.0 * t9 + t10 );
}
long double SqrtB08( long double x, long double c )
{
   long double t3 = x * x;
   long double t4 = t3 * t3;
   long double t5 = t4 * t4;
   long double t7 = t4 * t3 * c;
   long double t8 = c * c;
   long double t9 = t4 * t8;
   long double t11 = t3 * t8 * c;
   long double t12 = t8 * t8;
   return 1.
   /
   x * c * ( 9.0 * t5 + 84.0 * t7 + 126.0 * t9 + 36.0 * t11 + t12 )
   /
   ( t5 + 36.0 * t7 + 126.0 * t9 + 84.0 * t11 + 9.0 * t12 );
}
long double SqrtA09( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t7 = t3 * t2;
   long double t8 = c * c;
   long double t10 = t8 * c;
   long double t12 = t8 * t8;
   return 1.
   /
   x * ( t4 * t2 + 45.0 * t4 * c + 210.0 * t7 * t8 + 210.0 * t3 * t10 + 45.0 * t2 * t12 + t12 * c )
   /
   ( 5.0 * t4 + 60.0 * t7 * c + 126.0 * t3 * t8 + 60.0 * t2 * t10 + 5.0 * t12 ) * .5;
}
long double SqrtB09( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t5 = t3 * t2;
   long double t7 = c * c;
   long double t9 = t7 * c;
   long double t11 = t7 * t7;
   return 2.0 * x * c * ( 5.0 * t4 + 60.0 * t5 * c + 126.0 * t3 * t7 + 60.0 * t2 * t9 + 5.0 * t11 )
   /
   ( t4 * t2 + 45.0 * t4 * c + 210.0 * t5 * t7 + 210.0 * t3 * t9 + 45.0 * t2 * t11 + t11 * c );
}
long double SqrtA010( long double x, long double c )
{
   long double t1 = x * x;
   long double t2 = t1 * t1;
   long double t3 = t2 * t2;
   long double t4 = t3 * t1;
   long double t5 = t3 * c;
   long double t7 = c * c;
   long double t8 = t2 * t1 * t7;
   long double t10 = t2 * t7 * c;
   long double t11 = t7 * t7;
   long double t12 = t1 * t11;
   long double t13 = t11 * c;
   return x * ( t4 + 55.0 * t5 + 330.0 * t8 + 462.0 * t10 + 165.0 * t12 + 11.0 * t13 )
   /
   ( 11.0 * t4 + 165.0 * t5 + 462.0 * t8 + 330.0 * t10 + 55.0 * t12 + t13 );
}
long double SqrtB010( long double x, long double c )
{
   long double t3 = x * x;
   long double t4 = t3 * t3;
   long double t5 = t4 * t4;
   long double t6 = t5 * t3;
   long double t7 = t5 * c;
   long double t9 = c * c;
   long double t10 = t4 * t3 * t9;
   long double t12 = t4 * t9 * c;
   long double t13 = t9 * t9;
   long double t14 = t3 * t13;
   long double t15 = t13 * c;
   return 1.
   /
   x * c * ( 11.0 * t6 + 165.0 * t7 + 462.0 * t10 + 330.0 * t12 + 55.0 * t14 + t15 )
   /
   ( t6 + 55.0 * t7 + 330.0 * t10 + 462.0 * t12 + 165.0 * t14 + 11.0 * t15 );
}
long double SqrtA011( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t6 = t4 * t2;
   long double t8 = c * c;
   long double t10 = t3 * t2;
   long double t11 = t8 * c;
   long double t13 = t8 * t8;
   long double t15 = t13 * c;
   return 1.
   /
   x * ( t4 * t3 + 66.0 * t6 * c + 495.0 * t4 * t8 + 924.0 * t10 * t11 + 495.0 * t3 * t13 + 66.0 * t2 * t15 + t13 * t8 )
   /
   ( 3.0 * t6 + 55.0 * t4 * c + 198.0 * t10 * t8 + 198.0 * t3 * t11 + 55.0 * t2 * t13 + 3.0 * t15 ) * .25;
}
long double SqrtB011( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t5 = t4 * t2;
   long double t7 = t3 * t2;
   long double t8 = c * c;
   long double t10 = t8 * c;
   long double t12 = t8 * t8;
   long double t14 = t12 * c;
   return 4.0 * x * c * ( 3.0 * t5 + 55.0 * t4 * c + 198.0 * t7 * t8 + 198.0 * t3 * t10 + 55.0 * t2 * t12 + 3.0 * t14 )
   /
   ( t4 * t3 + 66.0 * t5 * c + 495.0 * t4 * t8 + 924.0 * t7 * t10 + 495.0 * t3 * t12 + 66.0 * t2 * t14 + t12 * t8 );
}
long double SqrtA012( long double x, long double c )
{
   long double t1 = x * x;
   long double t2 = t1 * t1;
   long double t3 = t2 * t2;
   long double t5 = t3 * t1 * c;
   long double t6 = t3 * t2;
   long double t8 = c * c;
   long double t10 = t2 * t1 * t8 * c;
   long double t11 = t3 * t8;
   long double t12 = t8 * t8;
   long double t13 = t12 * t8;
   long double t14 = t2 * t12;
   long double t16 = t1 * t12 * c;
   return x * ( 78.0 * t5 + t6 + 1716.0 * t10 + 715.0 * t11 + 13.0 * t13 + 1287.0 * t14 + 286.0 * t16 )
   /
   ( 286.0 * t5 + 13.0 * t6 + 1716.0 * t10 + 1287.0 * t11 + t13 + 715.0 * t14 + 78.0 * t16 );
}
long double SqrtB012( long double x, long double c )
{
   long double t3 = x * x;
   long double t4 = t3 * t3;
   long double t5 = t4 * t4;
   long double t7 = t5 * t3 * c;
   long double t8 = t5 * t4;
   long double t10 = c * c;
   long double t12 = t4 * t3 * t10 * c;
   long double t13 = t5 * t10;
   long double t14 = t10 * t10;
   long double t15 = t14 * t10;
   long double t16 = t4 * t14;
   long double t18 = t3 * t14 * c;
   return 1.
   /
   x * c * ( 286.0 * t7 + 13.0 * t8 + 1716.0 * t12 + 1287.0 * t13 + t15 + 715.0 * t16 + 78.0 * t18 )
   /
   ( 78.0 * t7 + t8 + 1716.0 * t12 + 715.0 * t13 + 13.0 * t15 + 1287.0 * t16 + 286.0 * t18 );
}
long double SqrtA013( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t5 = t4 * t3;
   long double t7 = t4 * t2;
   long double t8 = c * c;
   long double t10 = t8 * c;
   long double t12 = t3 * t2;
   long double t13 = t8 * t8;
   long double t15 = t13 * c;
   long double t18 = t13 * t8;
   return 1.
   /
   x * ( 91.0 * t5 * c + 1001.0 * t7 * t8 + 3003.0 * t4 * t10 + 3003.0 * t12 * t13 + 1001.0 * t3 * t15 + t4 * t12 + 91.0 * t2 * t18 + t13 * t10 )
   /
   ( 182.0 * t7 * c + 1001.0 * t4 * t8 + 1716.0 * t12 * t10 + 1001.0 * t3 * t13 + 182.0 * t2 * t15 + 7.0 * t5 + 7.0 * t18 ) * .5;
}
long double SqrtB013( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t5 = t4 * t3;
   long double t6 = t4 * t2;
   long double t8 = c * c;
   long double t9 = t8 * t8;
   long double t10 = t9 * c;
   long double t14 = t3 * t2;
   long double t15 = t8 * c;
   long double t17 = t9 * t8;
   return 2.0 * x * c * ( 7.0 * t5 + 182.0 * t6 * c + 182.0 * t2 * t10 + 1001.0 * t4 * t8 + 1001.0 * t3 * t9 + 1716.0 * t14 * t15 + 7.0 * t17 )
   /
   ( 91.0 * t5 * c + 1001.0 * t6 * t8 + 91.0 * t2 * t17 + 3003.0 * t4 * t15 + 1001.0 * t3 * t10 + 3003.0 * t14 * t9 + t9 * t15 + t4 * t14 );
}
long double SqrtA014( long double x, long double c )
{
   long double t1 = x * x;
   long double t2 = t1 * t1;
   long double t3 = t2 * t2;
   long double t5 = t3 * t2 * c;
   long double t7 = c * c;
   long double t8 = t3 * t1 * t7;
   long double t9 = t7 * t7;
   long double t11 = t1 * t9 * t7;
   long double t12 = t7 * c;
   long double t13 = t3 * t12;
   long double t15 = t2 * t9 * c;
   long double t16 = t2 * t1;
   long double t17 = t16 * t9;
   long double t18 = t9 * t12;
   long double t19 = t3 * t16;
   return x * ( 105.0 * t5 + 1365.0 * t8 + 455.0 * t11 + 5005.0 * t13 + 3003.0 * t15 + 6435.0 * t17 + 15.0 * t18 + t19 )
   /
   ( 455.0 * t5 + 3003.0 * t8 + 105.0 * t11 + 6435.0 * t13 + 1365.0 * t15 + 5005.0 * t17 + t18 + 15.0 * t19 );
}
long double SqrtB014( long double x, long double c )
{
   long double t3 = x * x;
   long double t4 = t3 * t3;
   long double t5 = t4 * t4;
   long double t7 = t5 * t4 * c;
   long double t9 = c * c;
   long double t10 = t5 * t3 * t9;
   long double t11 = t9 * t9;
   long double t13 = t3 * t11 * t9;
   long double t14 = t9 * c;
   long double t15 = t5 * t14;
   long double t17 = t4 * t11 * c;
   long double t18 = t4 * t3;
   long double t19 = t18 * t11;
   long double t20 = t11 * t14;
   long double t21 = t5 * t18;
   return 1.
   /
   x * c * ( 455.0 * t7 + 3003.0 * t10 + 105.0 * t13 + 6435.0 * t15 + 1365.0 * t17 + 5005.0 * t19 + t20 + 15.0 * t21 )
   /
   ( 105.0 * t7 + 1365.0 * t10 + 455.0 * t13 + 5005.0 * t15 + 3003.0 * t17 + 6435.0 * t19 + 15.0 * t20 + t21 );
}
long double SqrtA015( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t2;
   long double t5 = t3 * t3;
   long double t6 = t5 * t4;
   long double t8 = t5 * t3;
   long double t9 = c * c;
   long double t11 = t9 * t9;
   long double t12 = t11 * t9;
   long double t14 = t5 * t2;
   long double t15 = t9 * c;
   long double t17 = t11 * c;
   long double t20 = t11 * t15;
   long double t22 = t5 * t5;
   long double t23 = t11 * t11;
   return 1.
   /
   x * ( 120.0 * t6 * c + 1820.0 * t8 * t9 + 1820.0 * t3 * t12 + 8008.0 * t14 * t15 + 8008.0 * t4 * t17 + 12870.0 * t5 * t11 + 120.0 * t2 * t20 + t22 + t23 )
   /
   ( 35.0 * t8 * c + 273.0 * t14 * t9 + 35.0 * t2 * t12 + 715.0 * t5 * t15 + 273.0 * t3 * t17 + 715.0 * t4 * t11 + t20 + t6 ) * .0625;
}
long double SqrtB015( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t5 = t4 * t3;
   long double t7 = t4 * t2;
   long double t8 = c * c;
   long double t10 = t8 * t8;
   long double t11 = t10 * t8;
   long double t13 = t8 * c;
   long double t15 = t10 * c;
   long double t17 = t3 * t2;
   long double t19 = t10 * t13;
   long double t20 = t4 * t17;
   long double t29 = t4 * t4;
   long double t30 = t10 * t10;
   return 16.0 * x * c * ( 35.0 * t5 * c + 273.0 * t7 * t8 + 35.0 * t2 * t11 + 715.0 * t4 * t13 + 273.0 * t3 * t15 + 715.0 * t17 * t10 + t19 + t20 )
   /
   ( 120.0 * t20 * c + 1820.0 * t5 * t8 + 1820.0 * t3 * t11 + 8008.0 * t7 * t13 + 8008.0 * t17 * t15 + 12870.0 * t4 * t10 + 120.0 * t2 * t19 + t29 + t30 );
}
long double SqrtA016( long double x, long double c )
{
   long double t1 = x * x;
   long double t2 = t1 * t1;
   long double t3 = t2 * t2;
   long double t5 = c * c;
   long double t6 = t3 * t2 * t5;
   long double t8 = t5 * c;
   long double t9 = t3 * t1 * t8;
   long double t10 = t5 * t5;
   long double t11 = t3 * t10;
   long double t12 = t2 * t1;
   long double t14 = t12 * t10 * c;
   long double t16 = t2 * t10 * t5;
   long double t18 = t1 * t10 * t8;
   long double t20 = t3 * t12 * c;
   long double t21 = t10 * t10;
   long double t22 = t3 * t3;
   return x * ( 2380.0 * t6 + 12376.0 * t9 + 24310.0 * t11 + 19448.0 * t14 + 6188.0 * t16 + 680.0 * t18 + 136.0 * t20 + 17.0 * t21 + t22 )
   /
   ( 6188.0 * t6 + 19448.0 * t9 + 24310.0 * t11 + 12376.0 * t14 + 2380.0 * t16 + 136.0 * t18 + 680.0 * t20 + t21 + 17.0 * t22 );
}
long double SqrtB016( long double x, long double c )
{
   long double t3 = x * x;
   long double t4 = t3 * t3;
   long double t5 = t4 * t4;
   long double t7 = c * c;
   long double t8 = t4 * t5 * t7;
   long double t10 = t7 * c;
   long double t11 = t5 * t3 * t10;
   long double t12 = t7 * t7;
   long double t13 = t5 * t12;
   long double t14 = t4 * t3;
   long double t16 = t14 * t12 * c;
   long double t18 = t4 * t12 * t7;
   long double t20 = t3 * t12 * t10;
   long double t22 = t5 * t14 * c;
   long double t23 = t12 * t12;
   long double t24 = t5 * t5;
   return 1.
   /
   x * c * ( 6188.0 * t8 + 19448.0 * t11 + 24310.0 * t13 + 12376.0 * t16 + 2380.0 * t18 + 136.0 * t20 + 680.0 * t22 + t23 + 17.0 * t24 )
   /
   ( 2380.0 * t8 + 12376.0 * t11 + 24310.0 * t13 + 19448.0 * t16 + 6188.0 * t18 + 680.0 * t20 + 136.0 * t22 + 17.0 * t23 + t24 );
}
long double SqrtA017( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t2;
   long double t5 = t3 * t3;
   long double t6 = t5 * t4;
   long double t7 = c * c;
   long double t9 = t5 * t3;
   long double t10 = t7 * c;
   long double t12 = t5 * t2;
   long double t13 = t7 * t7;
   long double t15 = t13 * c;
   long double t17 = t13 * t7;
   long double t19 = t13 * t10;
   long double t21 = t5 * t5;
   long double t23 = t13 * t13;
   return 1.
   /
   x * ( 3060.0 * t6 * t7 + 18564.0 * t9 * t10 + 43758.0 * t12 * t13 + 43758.0 * t5 * t15 + 18564.0 * t4 * t17 + 3060.0 * t3 * t19 + 153.0 * t21 * c + 153.0 * t2 * t23 + t21 * t2 + t23 * c )
   /
   ( 4284.0 * t9 * t7 + 15912.0 * t12 * t10 + 24310.0 * t5 * t13 + 15912.0 * t4 * t15 + 4284.0 * t3 * t17 + 408.0 * t2 * t19 + 408.0 * t6 * c + 9.0 * t23 + 9.0 * t21 ) * .5;
}
long double SqrtB017( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t5 = t4 * t2;
   long double t6 = c * c;
   long double t7 = t6 * c;
   long double t9 = t6 * t6;
   long double t11 = t3 * t2;
   long double t12 = t4 * t11;
   long double t14 = t4 * t3;
   long double t16 = t4 * t4;
   long double t17 = t9 * t6;
   long double t19 = t9 * c;
   long double t21 = t9 * t7;
   long double t23 = t9 * t9;
   return 2.0 * x * c * ( 15912.0 * t5 * t7 + 24310.0 * t4 * t9 + 408.0 * t12 * c + 4284.0 * t14 * t6 + 9.0 * t16 + 4284.0 * t3 * t17 + 15912.0 * t11 * t19 + 408.0 * t2 * t21 + 9.0 * t23 )
   /
   ( 43758.0 * t5 * t9 + 43758.0 * t4 * t19 + 3060.0 * t12 * t6 + 18564.0 * t14 * t7 + 153.0 * t16 * c + 3060.0 * t3 * t21 + 18564.0 * t11 * t17 + 153.0 * t2 * t23 + t16 * t2 + t23 * c );
}
long double SqrtA018( long double x, long double c )
{
   long double t1 = x * x;
   long double t2 = t1 * t1;
   long double t3 = c * c;
   long double t4 = t3 * c;
   long double t5 = t3 * t3;
   long double t7 = t2 * t5 * t4;
   long double t8 = t2 * t1;
   long double t10 = t8 * t5 * t3;
   long double t11 = t2 * t2;
   long double t13 = t11 * t5 * c;
   long double t15 = t11 * t2 * t4;
   long double t16 = t11 * t11;
   long double t17 = t16 * c;
   long double t19 = t11 * t8 * t3;
   long double t21 = t11 * t1 * t5;
   long double t22 = t5 * t5;
   long double t23 = t1 * t22;
   long double t24 = t16 * t1;
   long double t25 = t22 * c;
   return x * ( 11628.0 * t7 + 50388.0 * t10 + 92378.0 * t13 + 27132.0 * t15 + 171.0 * t17 + 3876.0 * t19 + 75582.0 * t21 + 969.0 * t23 + t24 + 19.0 * t25 )
   /
   ( 3876.0 * t7 + 27132.0 * t10 + 75582.0 * t13 + 50388.0 * t15 + 969.0 * t17 + 11628.0 * t19 + 92378.0 * t21 + 171.0 * t23 + 19.0 * t24 + t25 );
}
long double SqrtB018( long double x, long double c )
{
   long double t3 = x * x;
   long double t4 = t3 * t3;
   long double t5 = t4 * t4;
   long double t6 = t5 * t5;
   long double t7 = t6 * t3;
   long double t8 = t4 * t3;
   long double t10 = c * c;
   long double t11 = t5 * t8 * t10;
   long double t12 = t6 * c;
   long double t13 = t10 * c;
   long double t14 = t10 * t10;
   long double t16 = t4 * t14 * t13;
   long double t17 = t14 * t14;
   long double t18 = t3 * t17;
   long double t20 = t5 * t4 * t13;
   long double t22 = t5 * t3 * t14;
   long double t24 = t8 * t14 * t10;
   long double t26 = t5 * t14 * c;
   long double t27 = t17 * c;
   return 1.
   /
   x * c * ( 19.0 * t7 + 11628.0 * t11 + 969.0 * t12 + 3876.0 * t16 + 171.0 * t18 + 50388.0 * t20 + 92378.0 * t22 + 27132.0 * t24 + 75582.0 * t26 + t27 )
   /
   ( 171.0 * t12 + 27132.0 * t20 + 3876.0 * t11 + 969.0 * t18 + 19.0 * t27 + 75582.0 * t22 + 92378.0 * t26 + t7 + 11628.0 * t16 + 50388.0 * t24 );
}
long double SqrtA019( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = c * c;
   long double t5 = t4 * t4;
   long double t6 = t5 * t5;
   long double t8 = t3 * t3;
   long double t9 = t8 * t8;
   long double t11 = t8 * t2;
   long double t12 = t5 * c;
   long double t14 = t8 * t3;
   long double t16 = t3 * t2;
   long double t17 = t8 * t16;
   long double t18 = t4 * c;
   long double t21 = t9 * t2;
   long double t23 = t6 * c;
   long double t25 = t5 * t18;
   long double t27 = t5 * t4;
   long double t30 = 4845.0 * t3 * t6 + t9 * t3 + 184756.0 * t11 * t12 + 125970.0 * t14 * t5 + 38760.0 * t17 * t18 + 4845.0 * t9 * t4 + 190.0 * t21 * c + 190.0 * t2 * t23 + 38760.0 * t16 * t25 + 125970.0 * t8 * t27 + t6 * t4;
   return 1.
   /
   x * t30
   /
   ( 285.0 * t2 * t6 + 5.0 * t21 + 41990.0 * t8 * t12 + 41990.0 * t11 * t5 + 19380.0 * t14 * t18 + 3876.0 * t17 * t4 + 285.0 * t9 * c + 5.0 * t23 + 3876.0 * t3 * t25 + 19380.0 * t16 * t27 ) * .25;
}
long double SqrtB019( long double x, long double c )
{
   long double t2 = x * x;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t5 = t4 * t4;
   long double t7 = t5 * t2;
   long double t8 = c * c;
   long double t9 = t8 * t8;
   long double t10 = t9 * t9;
   long double t11 = t10 * c;
   long double t12 = t4 * t2;
   long double t14 = t4 * t3;
   long double t15 = t8 * c;
   long double t17 = t3 * t2;
   long double t18 = t4 * t17;
   long double t21 = t9 * t15;
   long double t23 = t9 * t8;
   long double t25 = t9 * c;
   long double t39 = t5 * t3 + t10 * t8 + 4845.0 * t5 * t8 + 190.0 * t7 * c + 38760.0 * t18 * t15 + 184756.0 * t12 * t25 + 125970.0 * t14 * t9 + 125970.0 * t4 * t23 + 38760.0 * t17 * t21 + 4845.0 * t3 * t10 + 190.0 * t2 * t11;
   return 4.0 * x * c * ( 285.0 * t5 * c + 5.0 * t7 + 5.0 * t11 + 41990.0 * t12 * t9 + 19380.0 * t14 * t15 + 3876.0 * t18 * t8 + 285.0 * t2 * t10 + 3876.0 * t3 * t21 + 19380.0 * t17 * t23 + 41990.0 * t4 * t25 )
   /
   t39;
}
long double SqrtA020( long double x, long double c )
{
   long double t1 = x * x;
   long double t2 = t1 * t1;
   long double t3 = t2 * t2;
   long double t4 = t3 * t3;
   long double t5 = t4 * t2;
   long double t6 = c * c;
   long double t7 = t6 * t6;
   long double t8 = t7 * t7;
   long double t9 = t8 * t6;
   long double t10 = t4 * t6;
   long double t12 = t4 * t1 * c;
   long double t13 = t2 * t1;
   long double t15 = t6 * c;
   long double t16 = t3 * t13 * t15;
   long double t19 = t3 * t1 * t7 * c;
   long double t21 = t3 * t2 * t7;
   long double t23 = t3 * t7 * t6;
   long double t25 = t13 * t7 * t15;
   long double t26 = t2 * t8;
   long double t28 = t1 * t8 * c;
   long double t29 = t5 + 21.0 * t9 + 5985.0 * t10 + 210.0 * t12 + 54264.0 * t16 + 352716.0 * t19 + 203490.0 * t21 + 293930.0 * t23 + 116280.0 * t25 + 20349.0 * t26 + 1330.0 * t28;
   long double t31 = 21.0 * t5 + t9 + 20349.0 * t10 + 1330.0 * t12 + 116280.0 * t16 + 352716.0 * t19 + 293930.0 * t21 + 203490.0 * t23 + 54264.0 * t25 + 5985.0 * t26 + 210.0 * t28;
   return x * t29
   /
   t31;
}
long double SqrtB020( long double x, long double c )
{
   long double t3 = x * x;
   long double t4 = t3 * t3;
   long double t5 = t4 * t4;
   long double t6 = t5 * t5;
   long double t7 = c * c;
   long double t8 = t6 * t7;
   long double t9 = t4 * t3;
   long double t11 = t7 * c;
   long double t12 = t5 * t9 * t11;
   long double t13 = t7 * t7;
   long double t14 = t13 * t13;
   long double t15 = t4 * t14;
   long double t17 = t5 * t13 * t7;
   long double t19 = t9 * t13 * t11;
   long double t21 = t5 * t4 * t13;
   long double t24 = t5 * t3 * t13 * c;
   long double t25 = t6 * t4;
   long double t27 = t6 * t3 * c;
   long double t29 = t3 * t14 * c;
   long double t30 = t14 * t7;
   long double t31 = 20349.0 * t8 + 116280.0 * t12 + 5985.0 * t15 + 203490.0 * t17 + 54264.0 * t19 + 293930.0 * t21 + 352716.0 * t24 + 21.0 * t25 + 1330.0 * t27 + 210.0 * t29 + t30;
   long double t32 = 54264.0 * t12 + 203490.0 * t21 + 1330.0 * t29 + 116280.0 * t19 + 20349.0 * t15 + 352716.0 * t24 + 293930.0 * t17 + 210.0 * t27 + 5985.0 * t8 + t25 + 21.0 * t30;
   return 1.
   /
   x * c * t31
   /
   t32;
}