float_t
float_t p, q;
float_t p, q;
float_t w,s1,s2,z;
float_t y;
float_t ft;
float_t s,P,Q;
float_t s,R,S;
float_t y,hi,lo,c,t,e,hxs,hfx,r1,twopk;
float_t z;
float_t z,r,s;
float_t s,r,z;
float_t z,r,s;
float_t s,r,z;
float_t hfsq,f,s,z,R,w,t1,t2,dk,hi,lo;
float_t hfsq,f,c,s,z,R,w,t1,t2,dk;
float_t x2 = x * x;
static const float_t toint = 1/EPS;
float_t y;
static const float_t toint = 1/EPS;
float_t y;
float_t y = x;
float_t s, c;