static double const g6 = +1.589413637225924000396747428344197673e-10;
static double const h0 = -26.268288778188055898859139820039082208557;
static double const k0 = -1.08992067927738266857142057187739904444e9;
static double const k1 = +7.86385819193028103421370536822805774891e7;
static double const k2 = +4.41099062172012274030491314715096929537e6;
static double const k3 = -151605.59508390068630347221193128225766282;
static double const k4 = +6987.6926038225639201637412277215418962752;

double sin9(double x)
{
    double xp = g6 * x;
    double y = x * x + h0;
    double y2 = y * y;
    //__asm__("" : "+x" (y), "+x" (y2), "+x" (xp));
    return (k0 + k1 * y + y2 * (k2 + k3 * y + y2 * (k4 + y2))) * xp;
}