polytail
return (x*(t-c1)-(ch*t-polytail(c1/x)));
static long double polytail(long double);
return x*(t-c1)-(ch*t-polytail(c1/x));