N a p l o g

見え方のテスト。

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

#define M_PI        3.141592

double term( int seriese, double x )
{
    return sin((2*seriese+1)*x) / (2*seriese+1);
}

int main(int argc, char *argv[])
{
    double  x, y;
    FILE    *fp;
    int     seriese, seriese_max;
    
    if( argc == 2 )
    {
        seriese_max = atoi(argv[1]);
        if( seriese_max == 0 )
            seriese_max = 10;
    }
    else
    {
        seriese_max = 10;
    }
    
    fp = fopen( "fourier.dat", "w" );
    
    for( x = -M_PI; x < M_PI; x += 0.001 )
    {
        y = 0;
        for( seriese = 0; seriese <= seriese_max; seriese++ )
        {
            y += term( seriese, x );
        }
        fprintf( fp, "%lf %lf\n", x, 4/M_PI*y );
    } 

    fclose(fp);
    
    return EXIT_SUCCESS;
}

ここで終わり。
そろそろちゃんとgnuplotの覚書を書いていこう。。。