Changeset 20668


Ignore:
Timestamp:
05/30/16 12:20:23 (9 years ago)
Author:
Mathieu Morlighem
Message:

NEW: more tests for Earl

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/c/shared/FSanalyticals/fsanalyticals.cpp

    r20663 r20668  
    99 * case 6: 3d test with sinusoidal functions, non homogeneous Dirichlet conditions
    1010 *
    11  * case 201: 3d test quadratic functions, EARL flow law
     11 * case 201: 3d test quadratic functions FS, EARL flow law
     12 * case 202: 3d test quadratic functions HO, EARL flow law
    1213 */
    1314
     
    7576                case 201:
    7677                        return pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 1.0L/3.0L) + y - 1.0L/18.0L*(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(2*pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 2*x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 8*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1) + x*pow(y, 2)*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*y*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 8*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1) + pow(x, 2)*y*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + x*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + (2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1) + 2*pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2) - 2*x*pow(y, 2)*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) - pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 4.0L/3.0L) - 2.0L/9.0L*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*((x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2) + x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2) + 2*x*pow(y, 2)*(x - 1)*pow(y - 1, 2)*(z - 1) + x*pow(y, 2)*(2*x - 1)*pow(y - 1, 2)*(z - 1) + pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 2*(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1) + pow(x, 2)*y*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + x*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 8*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*x*pow(y, 2)*z*(x - 1)*(y - 1)*(z - 1) + x*pow(y, 2)*z*(2*x - 1)*(y - 1)*(z - 1) + 2*x*y*z*(x - 1)*pow(y - 1, 2)*(z - 1) + x*y*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) + pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 8*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(pow(x, 2)*y*z*(y - 1)*(2*y - 1)*(z - 1) + 4*x*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 2*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*pow(x, 2)*y*z*(x - 1)*(y - 1)*(z - 1) + pow(x, 2)*y*z*(x - 1)*(2*y - 1)*(z - 1) + pow(x, 2)*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*pow(y - 1, 2)*(z - 1) + 2*x*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + x*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1)))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 4.0L/3.0L) + 2*(pow(x, 2)*y*z*(y - 1)*(2*y - 1)*(z - 1) + 4*x*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 1.0L/3.0L) - 1.0L/18.0L*(2*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1) + x*pow(y, 2)*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*y*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + (pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1) + 2*pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 8*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(z - 1) + 4*x*y*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 8*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(2*pow(x, 2)*y*z*(x - 1)*(y - 1)*(z - 1) + pow(x, 2)*y*z*(x - 1)*(2*y - 1)*(z - 1) + pow(x, 2)*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + 2*x*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + x*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 2*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(z - 1) + 2*pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*(y - 1)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*(y - 1)*(z - 1) - 2*x*y*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*y*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) - y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)))*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 4.0L/3.0L) - 1.0L/2.0L + (2*pow(x, 2)*y*z*pow(x - 1, 2)*(z - 1) + 2*pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*(y - 1)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*(y - 1)*(z - 1) - 2*x*y*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*y*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) - y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 1.0L/3.0L);
    77 
     78                case 202:
     79                        return -x*y*(x - 1)*(y - 1)/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 1.0L/3.0L) - 4*y*z*(y - 1)*(z - 1)/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 1.0L/3.0L) + (1.0L/18.0L)*(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1))*(10*x*y*(x - 1)*(y - 1)*(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1)) + 8*(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(x - 1) + x*y*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1) + x*(x - 1)*(y - 1)*(z - 1)) + 2*(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(y - 1) + x*y*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1) + y*(x - 1)*(y - 1)*(z - 1)) + (x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(x - 1) + 2*x*y*z*(y - 1) + x*y*(x - 1)*(z - 1) + 2*x*y*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1) + x*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1) + 2*y*(x - 1)*(y - 1)*(z - 1)) + 2*(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(2*x*y*z*(x - 1) + x*y*z*(y - 1) + 2*x*y*(x - 1)*(z - 1) + x*y*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1) + 2*x*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1) + y*(x - 1)*(y - 1)*(z - 1)))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 4.0L/3.0L) + (2.0L/9.0L)*(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1))*(4*y*z*(y - 1)*(z - 1)*(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1)) + 5*(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(y - 1) + x*y*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1) + y*(x - 1)*(y - 1)*(z - 1)) + 8*(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + (x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + 4*y*z*(y - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + 4*(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + y*z*(y - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 4.0L/3.0L) + (2.0L/9.0L)*(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(4*y*z*(y - 1)*(z - 1)*(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1)) + 5*(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(y - 1) + x*y*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1) + y*(x - 1)*(y - 1)*(z - 1)) + 8*(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + (x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + 4*y*z*(y - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + 4*(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + y*z*(y - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 4.0L/3.0L) - 2*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 1.0L/3.0L) + (1.0L/18.0L)*(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1))*(16*x*z*(x - 1)*(z - 1)*(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1)) + 5*(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(x - 1) + x*y*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1) + x*(x - 1)*(y - 1)*(z - 1)) + 2*(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + 2*(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(x - 1)*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + 2*(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + 4*x*z*(x - 1)*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 4.0L/3.0L) - (x*y*z*(z - 1) + x*z*(x - 1)*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 1.0L/3.0L);
     80                case 211:
     81                        return (200.0L/707238535107.0L)*pow(17, 2.0L/3.0L)*pow(6, 1.0L/3.0L)*(-9*(169*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))*(2197*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-2197*z/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 26*z*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 2) + sqrt(13)*z*(-78*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 169*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 13*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 169*sqrt(13)*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L)) + 13*(9*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 78*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 13*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))*(9*sqrt(13)*z*(-2197*z/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 26*z*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 2) + sqrt(13)*z*(-78*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 169*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 13*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 117*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 1014*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 2197*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 13*sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L)))/pow((8.0L/17.0L)*pow(fabs((15.0L/26.0L)*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (15.0L/4394.0L)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + (8.0L/17.0L)*pow(fabs(-45.0L/2197.0L*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (30.0L/169.0L)*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 5.0L/169.0L*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + (8.0L/17.0L)*pow(fabs((135.0L/4394.0L)*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 45.0L/169.0L*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (15.0L/338.0L)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + 1, 4.0L/3.0L) - 150.0L/235746178369.0L*pow(17, 2.0L/3.0L)*pow(6, 1.0L/3.0L)*(3*(169*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))*(507*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 3*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-6*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 2) + 3*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 6591*z/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 3*sqrt(13)*z*(6*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - (507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 169*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 169*sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 169*sqrt(13)*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 507*sqrt(13)*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 28561*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))) + 13*(9*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 78*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 13*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))*(9*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 78*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 3*sqrt(13)*z*(-6*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 2) + 3*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 6591*z/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 3*sqrt(13)*z*(6*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - (507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 169*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 169*sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 169*sqrt(13)*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 13*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 507*sqrt(13)*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 2197*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))))/pow((8.0L/17.0L)*pow(fabs((15.0L/26.0L)*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (15.0L/4394.0L)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + (8.0L/17.0L)*pow(fabs(-45.0L/2197.0L*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (30.0L/169.0L)*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 5.0L/169.0L*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + (8.0L/17.0L)*pow(fabs((135.0L/4394.0L)*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 45.0L/169.0L*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (15.0L/338.0L)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + 1, 4.0L/3.0L);
     82
     83                case 212:
     84                        return (16.0L/9.0L)*x*(4*x + 1)/(pow(2*(4*pow(fabs((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z)), 2) + 9)/(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11) + 1, 1.0L/3.0L)*pow((4.0L/3.0L)*pow(x, 2) + (1.0L/3.0L)*pow(2*x + 1, 2) + 11.0L/6.0L, 4.0L/3.0L)) + (64.0L/3.0L)*x*(-(4*x + 1)*(4*pow(fabs((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z)), 2) + 9)/pow(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11, 2) + 2*(pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (y + 3*z)/fabs(y + 3*z))*((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z))/(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11))/(pow(2*(4*pow(fabs((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z)), 2) + 9)/(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11) + 1, 4.0L/3.0L)*pow((4.0L/3.0L)*pow(x, 2) + (1.0L/3.0L)*pow(2*x + 1, 2) + 11.0L/6.0L, 1.0L/3.0L)) + (4.0L/9.0L)*(4*x + 1)/(pow(2*(4*pow(fabs((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z)), 2) + 9)/(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11) + 1, 1.0L/3.0L)*pow((4.0L/3.0L)*pow(x, 2) + (1.0L/3.0L)*pow(2*x + 1, 2) + 11.0L/6.0L, 4.0L/3.0L)) - 4/(pow(2*(4*pow(fabs((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z)), 2) + 9)/(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11) + 1, 1.0L/3.0L)*pow((4.0L/3.0L)*pow(x, 2) + (1.0L/3.0L)*pow(2*x + 1, 2) + 11.0L/6.0L, 1.0L/3.0L)) + (16.0L/3.0L)*(-(4*x + 1)*(4*pow(fabs((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z)), 2) + 9)/pow(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11, 2) + 2*(pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (y + 3*z)/fabs(y + 3*z))*((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z))/(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11))/(pow(2*(4*pow(fabs((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z)), 2) + 9)/(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11) + 1, 4.0L/3.0L)*pow((4.0L/3.0L)*pow(x, 2) + (1.0L/3.0L)*pow(2*x + 1, 2) + 11.0L/6.0L, 1.0L/3.0L));
    7885                default:
    7986                        _error_("FS analytical solution"<<testid<<" not implemented yet");
     
    138145                case 201:
    139146                        return -x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 1.0L/3.0L) + x + (1.0L/18.0L)*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 2*x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 8*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1) + x*pow(y, 2)*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*y*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 8*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1) + pow(x, 2)*y*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + x*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + (2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1) + 2*pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2) - 2*x*pow(y, 2)*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) - pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 4.0L/3.0L) + (2.0L/9.0L)*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1) + x*pow(y, 2)*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*y*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + (pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1) + 2*pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 8*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(z - 1) + 4*x*y*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 8*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(2*pow(x, 2)*y*z*(x - 1)*(y - 1)*(z - 1) + pow(x, 2)*y*z*(x - 1)*(2*y - 1)*(z - 1) + pow(x, 2)*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + 2*x*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + x*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 2*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(z - 1) + 2*pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*(y - 1)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*(y - 1)*(z - 1) - 2*x*y*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*y*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) - y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 4.0L/3.0L) - 2*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(z - 1) + 4*x*y*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 1.0L/3.0L) - 1.0L/18.0L*((x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2) + x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2) + 2*x*pow(y, 2)*(x - 1)*pow(y - 1, 2)*(z - 1) + x*pow(y, 2)*(2*x - 1)*pow(y - 1, 2)*(z - 1) + pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 2*(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1) + pow(x, 2)*y*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + x*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 8*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*x*pow(y, 2)*z*(x - 1)*(y - 1)*(z - 1) + x*pow(y, 2)*z*(2*x - 1)*(y - 1)*(z - 1) + 2*x*y*z*(x - 1)*pow(y - 1, 2)*(z - 1) + x*y*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) + pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 8*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(pow(x, 2)*y*z*(y - 1)*(2*y - 1)*(z - 1) + 4*x*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 2*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*pow(x, 2)*y*z*(x - 1)*(y - 1)*(z - 1) + pow(x, 2)*y*z*(x - 1)*(2*y - 1)*(z - 1) + pow(x, 2)*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*pow(y - 1, 2)*(z - 1) + 2*x*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + x*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1)))*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 4.0L/3.0L) - 1.0L/2.0L + (2*pow(x, 2)*y*z*(x - 1)*(y - 1)*(z - 1) + pow(x, 2)*y*z*(x - 1)*(2*y - 1)*(z - 1) + pow(x, 2)*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*pow(y - 1, 2)*(z - 1) + 2*x*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + x*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 1.0L/3.0L);
     147                case 202:
     148                        return -2*x*y*(x - 1)*(y - 1)/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 1.0L/3.0L) - 8*x*z*(x - 1)*(z - 1)/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 1.0L/3.0L) + (1.0L/9.0L)*(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1))*(10*x*y*(x - 1)*(y - 1)*(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1)) + 8*(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(x - 1) + x*y*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1) + x*(x - 1)*(y - 1)*(z - 1)) + 2*(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(y - 1) + x*y*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1) + y*(x - 1)*(y - 1)*(z - 1)) + (x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(x - 1) + 2*x*y*z*(y - 1) + x*y*(x - 1)*(z - 1) + 2*x*y*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1) + x*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1) + 2*y*(x - 1)*(y - 1)*(z - 1)) + 2*(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(2*x*y*z*(x - 1) + x*y*z*(y - 1) + 2*x*y*(x - 1)*(z - 1) + x*y*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1) + 2*x*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1) + y*(x - 1)*(y - 1)*(z - 1)))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 4.0L/3.0L) + (4.0L/9.0L)*(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1))*(16*x*z*(x - 1)*(z - 1)*(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1)) + 5*(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(x - 1) + x*y*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1) + x*(x - 1)*(y - 1)*(z - 1)) + 2*(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + 2*(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(x - 1)*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + 2*(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + 4*x*z*(x - 1)*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 4.0L/3.0L) + (1.0L/9.0L)*(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(16*x*z*(x - 1)*(z - 1)*(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1)) + 5*(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(x - 1) + x*y*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1) + x*(x - 1)*(y - 1)*(z - 1)) + 2*(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + 2*(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(x - 1)*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + 2*(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + 4*x*z*(x - 1)*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 4.0L/3.0L) - (x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 1.0L/3.0L) + (1.0L/18.0L)*(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1))*(4*y*z*(y - 1)*(z - 1)*(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1)) + 5*(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(y - 1) + x*y*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1) + y*(x - 1)*(y - 1)*(z - 1)) + 8*(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + (x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + 4*y*z*(y - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)) + 4*(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1))*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + y*z*(y - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1)))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 4.0L/3.0L) - 1.0L/2.0L*(x*y*z*(z - 1) + x*z*(y - 1)*(z - 1) + y*z*(x - 1)*(z - 1) + 4*y*z*(y - 1)*(z - 1) + z*(x - 1)*(y - 1)*(z - 1))/pow((5.0L/6.0L)*pow(x*y*z*(x - 1)*(y - 1) + x*y*(x - 1)*(y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*y*z*(x - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(x*y*z*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(x*y*z*(x - 1)*(z - 1) + 2*x*y*z*(y - 1)*(z - 1) + x*z*(x - 1)*(y - 1)*(z - 1) + 2*y*z*(x - 1)*(y - 1)*(z - 1), 2) + (1.0L/3.0L)*pow(2*x*y*z*(x - 1)*(z - 1) + x*y*z*(y - 1)*(z - 1) + 2*x*z*(x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(y - 1)*(z - 1), 2), 1.0L/3.0L);
     149                case 211:
     150                        return (200.0L/235746178369.0L)*pow(17, 2.0L/3.0L)*pow(6, 1.0L/3.0L)*(-9*(169*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))*(2197*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-2197*z/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 26*z*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 2) + sqrt(13)*z*(-78*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 169*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 13*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 169*sqrt(13)*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L)) + 13*(9*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 78*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 13*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))*(9*sqrt(13)*z*(-2197*z/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 26*z*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 2) + sqrt(13)*z*(-78*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 169*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 13*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 117*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 1014*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 2197*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 13*sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L)))/pow((8.0L/17.0L)*pow(fabs((15.0L/26.0L)*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (15.0L/4394.0L)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + (8.0L/17.0L)*pow(fabs(-45.0L/2197.0L*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (30.0L/169.0L)*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 5.0L/169.0L*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + (8.0L/17.0L)*pow(fabs((135.0L/4394.0L)*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 45.0L/169.0L*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (15.0L/338.0L)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + 1, 4.0L/3.0L);
     151
     152                case 212:
     153                                return -8.0L/3.0L*x*(9*(y + 3*z)*(pow(y + 3*z, 2)/pow(fabs(y + 3*z), 3) - 1/fabs(y + 3*z))/fabs(y + 3*z) + 4*((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z))*(3*(2*x + 1)*pow(y + 3*z, 4)/pow(fabs(y + 3*z), 5) - 4*(2*x + 1)*pow(y + 3*z, 2)/pow(fabs(y + 3*z), 3) + (2*x + 1)/fabs(y + 3*z)))/(pow(2*(4*pow(fabs((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z)), 2) + 9)/(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11) + 1, 4.0L/3.0L)*pow((4.0L/3.0L)*pow(x, 2) + (1.0L/3.0L)*pow(2*x + 1, 2) + 11.0L/6.0L, 1.0L/3.0L)*(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11)) - 26.0L/3.0L*(9*(y + 3*z)*(pow(y + 3*z, 2)/pow(fabs(y + 3*z), 3) - 1/fabs(y + 3*z))/fabs(y + 3*z) + 4*((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z))*(3*(2*x + 1)*pow(y + 3*z, 4)/pow(fabs(y + 3*z), 5) - 4*(2*x + 1)*pow(y + 3*z, 2)/pow(fabs(y + 3*z), 3) + (2*x + 1)/fabs(y + 3*z)))/(pow(2*(4*pow(fabs((2*x + 1)*pow(y + 3*z, 3)/pow(fabs(y + 3*z), 3) - (2*x + 1)*(y + 3*z)/fabs(y + 3*z)), 2) + 9)/(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11) + 1, 4.0L/3.0L)*pow((4.0L/3.0L)*pow(x, 2) + (1.0L/3.0L)*pow(2*x + 1, 2) + 11.0L/6.0L, 1.0L/3.0L)*(8*pow(x, 2) + 2*pow(2*x + 1, 2) + 11));
     154
    140155                default:
    141156                        _error_("FS analytical solution"<<testid<<" not implemented yet");
     
    167182                case 201:
    168183                        return (1.0L/18.0L)*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1) + x*pow(y, 2)*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*y*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + (pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1) + 2*pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 8*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(z - 1) + 4*x*y*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 8*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(2*pow(x, 2)*y*z*(x - 1)*(y - 1)*(z - 1) + pow(x, 2)*y*z*(x - 1)*(2*y - 1)*(z - 1) + pow(x, 2)*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + 2*x*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + x*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 2*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(z - 1) + 2*pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*(y - 1)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*(y - 1)*(z - 1) - 2*x*y*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*y*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) - y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 4.0L/3.0L) - 1.0L/18.0L*(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*((x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2) + x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2) + 2*x*pow(y, 2)*(x - 1)*pow(y - 1, 2)*(z - 1) + x*pow(y, 2)*(2*x - 1)*pow(y - 1, 2)*(z - 1) + pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 2*(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1) + pow(x, 2)*y*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + x*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 8*(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*x*pow(y, 2)*z*(x - 1)*(y - 1)*(z - 1) + x*pow(y, 2)*z*(2*x - 1)*(y - 1)*(z - 1) + 2*x*y*z*(x - 1)*pow(y - 1, 2)*(z - 1) + x*y*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) + pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1)) + 8*(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))*(pow(x, 2)*y*z*(y - 1)*(2*y - 1)*(z - 1) + 4*x*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1)) + 2*(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))*(2*pow(x, 2)*y*z*(x - 1)*(y - 1)*(z - 1) + pow(x, 2)*y*z*(x - 1)*(2*y - 1)*(z - 1) + pow(x, 2)*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*pow(y - 1, 2)*(z - 1) + 2*x*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + x*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1)))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 4.0L/3.0L) - (x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1) + x*pow(y, 2)*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*y*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 1.0L/3.0L) + (pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1) + pow(x, 2)*y*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + x*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1))/pow((1.0L/6.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2) + x*pow(y, 2)*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (1.0L/6.0L)*pow(pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1) + pow(x, 2)*y*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (4.0L/3.0L)*pow(x*pow(y, 2)*z*(x - 1)*(2*x - 1)*(y - 1)*(z - 1) + x*y*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2) + (4.0L/3.0L)*pow(pow(x, 2)*y*z*(x - 1)*(y - 1)*(2*y - 1)*(z - 1) + x*y*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1), 2) + (1.0L/6.0L)*pow(2*pow(x, 2)*y*z*pow(x - 1, 2)*(y - 1)*(z - 1) + pow(x, 2)*y*z*pow(x - 1, 2)*(2*y - 1)*(z - 1) + pow(x, 2)*z*pow(x - 1, 2)*(y - 1)*(2*y - 1)*(z - 1) - 2*x*pow(y, 2)*z*(x - 1)*pow(y - 1, 2)*(z - 1) - x*pow(y, 2)*z*(2*x - 1)*pow(y - 1, 2)*(z - 1) - pow(y, 2)*z*(x - 1)*(2*x - 1)*pow(y - 1, 2)*(z - 1), 2), 1.0L/3.0L);
     184                case 211:
     185                        return (200.0L/235746178369.0L)*pow(17, 2.0L/3.0L)*pow(6, 1.0L/3.0L)*(-9*(169*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))*(2197*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-2197*z/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 26*z*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 2) + sqrt(13)*z*(-78*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 169*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 13*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 169*sqrt(13)*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L)) + 13*(9*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 78*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 13*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))*(9*sqrt(13)*z*(-2197*z/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 26*z*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 2) + sqrt(13)*z*(-78*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 169*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 13*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 117*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 1014*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 2197*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 13*sqrt(13)*pow(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y), 2)/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L)))/pow((8.0L/17.0L)*pow(fabs((15.0L/26.0L)*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (15.0L/4394.0L)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + (8.0L/17.0L)*pow(fabs(-45.0L/2197.0L*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (30.0L/169.0L)*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 5.0L/169.0L*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + (8.0L/17.0L)*pow(fabs((135.0L/4394.0L)*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 45.0L/169.0L*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (15.0L/338.0L)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + 1, 4.0L/3.0L) + (100.0L/235746178369.0L)*pow(17, 2.0L/3.0L)*pow(6, 1.0L/3.0L)*(3*(169*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))*(507*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 3*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-6*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 2) + 3*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 6591*z/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 3*sqrt(13)*z*(6*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - (507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 169*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 169*sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 169*sqrt(13)*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 507*sqrt(13)*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 28561*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))) + 13*(9*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 78*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 13*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))*(9*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 78*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 3*sqrt(13)*z*(-6*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 2) + 3*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) + 6591*z/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 3*sqrt(13)*z*(6*sqrt(13)*z*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - (507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 169*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 169*sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 169*sqrt(13)*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 13*(507*z + sqrt(13)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y)))*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/pow(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2), 3.0L/2.0L) - 507*sqrt(13)*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + 2197*sqrt(13)/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))))/pow((8.0L/17.0L)*pow(fabs((15.0L/26.0L)*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (15.0L/4394.0L)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + (8.0L/17.0L)*pow(fabs(-45.0L/2197.0L*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (30.0L/169.0L)*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 5.0L/169.0L*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + (8.0L/17.0L)*pow(fabs((135.0L/4394.0L)*sqrt(13)*z*(-sqrt(13)*z*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + sqrt(13)*z*(6*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - (3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) - 45.0L/169.0L*sqrt(13)*z/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2)) + (15.0L/338.0L)*(3*sqrt(13)*z + 2*sqrt(13)*(2*x + y))/sqrt(9*pow(z, 2) + pow(fabs((3.0L/13.0L)*sqrt(13)*z + (2.0L/13.0L)*sqrt(13)*(2*x + y)), 2))), 2) + 1, 4.0L/3.0L);
    169186
    170187                default:
Note: See TracChangeset for help on using the changeset viewer.