#include    "wfl.h"
using namespace wfl;

void main()
{
    Start();
    WaveField::SetDefault(512, 1);  // デフォルト値の設定
    WaveField::SetDefaultPx(0.01);  // デフォルト値の設定

    WaveField a, b;                 // デフォルト値で生成するコンストラクタ

    double f = 0.1;                 // 空間周波数 0.1 [1/m]
    double omega = 2 * Pi * f;      // 角周波数
    int i;
    b.Clear();                      // ゼロクリア
    for (i = 0; i < a.GetNx(); i++) // i < 512 のループ
    {
        Complex c = (0, omega*a.X(i));
        a.SetPixel(i, 0, exp(c));   // 点iのサンプリング値に複素数値exp(i2πfx)を設定
        b.SetReal(i, 0, b.X(i)*b.X(i)); // 点iの実部にxの2乗を設定
    }
    a *= b;                         // a = a * b
    a *= 3.0;                       // a = a * 3
    a.SaveAsCsv("function.csv");    // CSVファイルとして保存
}