#include <stdlib.h>
#include <stdio.h>
#include <string.h>
#include <math.h>
 
float Series(float x, unsigned int N)
{
    x = fmod(x,2*3.1415925358979323846);
    float term = x*x*x;
    x *= -x;
    float sum = term, tr = 8;
 
    for(int n = 2; n <= N; ++n)
    {
        term /= tr;
        tr = (tr+1)*9-1;
        term *= x*tr/(2*n*(2*n+1));
        sum += term;
    }
    return sum;
}
 
 
int main(int argc, const char * argv[])
{
    for(float x = 0.1; x < 1500; x = x*2)
    {
        printf("%f  %f\n",Series(x,30),sin(x)*sin(x)*sin(x));
    }
}
 

I2luY2x1ZGUgPHN0ZGxpYi5oPgojaW5jbHVkZSA8c3RkaW8uaD4KI2luY2x1ZGUgPHN0cmluZy5oPgojaW5jbHVkZSA8bWF0aC5oPgoKZmxvYXQgU2VyaWVzKGZsb2F0IHgsIHVuc2lnbmVkIGludCBOKQp7CiAgICB4ID0gZm1vZCh4LDIqMy4xNDE1OTI1MzU4OTc5MzIzODQ2KTsKICAgIGZsb2F0IHRlcm0gPSB4KngqeDsKICAgIHggKj0gLXg7CiAgICBmbG9hdCBzdW0gPSB0ZXJtLCB0ciA9IDg7CgogICAgZm9yKGludCBuID0gMjsgbiA8PSBOOyArK24pCiAgICB7CiAgICAgICAgdGVybSAvPSB0cjsKICAgICAgICB0ciA9ICh0cisxKSo5LTE7CiAgICAgICAgdGVybSAqPSB4KnRyLygyKm4qKDIqbisxKSk7CiAgICAgICAgc3VtICs9IHRlcm07CiAgICB9CiAgICByZXR1cm4gc3VtOwp9CgoKaW50IG1haW4oaW50IGFyZ2MsIGNvbnN0IGNoYXIgKiBhcmd2W10pCnsKICAgIGZvcihmbG9hdCB4ID0gMC4xOyB4IDwgMTUwMDsgeCA9IHgqMikKICAgIHsKICAgICAgICBwcmludGYoIiVmICAlZlxuIixTZXJpZXMoeCwzMCksc2luKHgpKnNpbih4KSpzaW4oeCkpOwogICAgfQp9Cg==