#include    <wfl.h>
using namespace wfl;

void main(void)
{
    Start();
    WaveField::SetDefault(1024, 1, 0.01);
    WaveField a, b;
    int i;

    // 矩形関数をaに設定
    a.Clear();
    for (i = a.I(-0.5); i < a.I(+0.5); i++)
        a.SetReal(i, 0, 1.0);
    a.SaveAsCsv("Rect.csv");
    a.Fft(-1);                      // 高速フーリエ変換
    a.SaveAsCsv("FFT-Rect.csv");

    // Λ(ラムダ)関数をbに設定
    b.Clear();
    for (i = b.I(-1.0); i < b.I(+1.0); i++)
        b.SetReal(i, 0, 1.0 - fabs(b.X(i)));
    b.SaveAsCsv("lambda.csv");
    b.Fft(-1);                      // 高速フーリエ変換
    b.SaveAsCsv("FFT-lambda.csv");
}