static inline
double eval_K_rtheta(double a, double L, double r_0, double x_0, double x_1) {

   double eval_K_rtheta_result;
   eval_K_rtheta_result = (3.0L/512.0L)*sqrt(143)*a*(5*pow(x_0, 3)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) - 16*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 16*pow(r_0 - x_0, 4)) + 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 4) - pow(pow(L, 2) - pow(r_0 - x_0, 2), 5)*(r_0 - x_0))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 7)*pow(x_0, 4));
   return eval_K_rtheta_result;

}

static inline
double eval_dK_rtheta_r(double a, double L, double r_0, double x_0, double x_1) {

   double eval_dK_rtheta_r_result;
   eval_dK_rtheta_r_result = (3.0L/512.0L)*sqrt(143)*a*(5*pow(x_0, 3)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(-16*(pow(L, 2) - pow(r_0 - x_0, 2))*(-2*r_0 + 2*x_0) + (pow(L, 2) - pow(r_0 - x_0, 2))*(4*r_0 - 4*x_0) - 64*pow(r_0 - x_0, 3) - 16*pow(r_0 - x_0, 2)*(2*r_0 - 2*x_0)) + 5*pow(x_0, 3)*(pow(L, 2) - pow(r_0 - x_0, 2))*(4*r_0 - 4*x_0)*(pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) - 16*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 16*pow(r_0 - x_0, 4)) - 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2)) + 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(r_0 - x_0)*(22*r_0 - 22*x_0) + 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0)*(6*r_0 - 6*x_0) + 15*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) - 16*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 16*pow(r_0 - x_0, 4)) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(8*r_0 - 8*x_0) - x_0*pow(pow(L, 2) - pow(r_0 - x_0, 2), 4)*(22*r_0 - 22*x_0) + 10*x_0*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) - (pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 4) + pow(pow(L, 2) - pow(r_0 - x_0, 2), 5) - pow(pow(L, 2) - pow(r_0 - x_0, 2), 4)*(r_0 - x_0)*(10*r_0 - 10*x_0))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 7)*pow(x_0, 4)) - 3.0L/128.0L*sqrt(143)*a*(5*pow(x_0, 3)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) - 16*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 16*pow(r_0 - x_0, 4)) + 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 4) - pow(pow(L, 2) - pow(r_0 - x_0, 2), 5)*(r_0 - x_0))*sin(2*x_1)/(sqrt(M_PI)*pow(L, 7)*pow(x_0, 5));
   return eval_dK_rtheta_r_result;

}

static inline
double eval_dK_rtheta_theta(double a, double L, double r_0, double x_0, double x_1) {

   double eval_dK_rtheta_theta_result;
   eval_dK_rtheta_theta_result = (3.0L/256.0L)*sqrt(143)*a*(5*pow(x_0, 3)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 2)*(pow(pow(L, 2) - pow(r_0 - x_0, 2), 2) - 16*(pow(L, 2) - pow(r_0 - x_0, 2))*pow(r_0 - x_0, 2) + 16*pow(r_0 - x_0, 4)) + 5*pow(x_0, 2)*pow(pow(L, 2) - pow(r_0 - x_0, 2), 3)*(3*pow(L, 2) - 11*pow(r_0 - x_0, 2))*(r_0 - x_0) - x_0*(pow(L, 2) - 11*pow(r_0 - x_0, 2))*pow(pow(L, 2) - pow(r_0 - x_0, 2), 4) - pow(pow(L, 2) - pow(r_0 - x_0, 2), 5)*(r_0 - x_0))*cos(2*x_1)/(sqrt(M_PI)*pow(L, 7)*pow(x_0, 4));
   return eval_dK_rtheta_theta_result;

}