#include    "wfl.h"
using namespace wfl;

void main()
{
    Start();
    WaveField a(256, 1, 0.01), b(1024, 1, 0.02);

    double f = 0.1;                     // 空間周波数 0.1 [1/m]
    double omega = 2 * Pi * f;          // 角周波数

    int i;
    a.Clear();                          // aをゼロクリア
    for (i = 0; i < a.GetNx(); i++)     // i < 256 のループ
    {
        double x = a.X(i);              // サンプリング点iの物理座標x
        a.SetReal(i, 0, cos(omega*x));  // 点iのサンプリング値の実部としてcos(x)を設定
    }

    b.Clear();
    for (i = 0; i < b.GetNx(); i++)     // i < 1024 のループ
    {
        double x = b.X(i);              // サンプリング点iの物理座標x
        b.SetImag(i, 0, sin(omega*x));  // 点iのサンプリング値の虚部としてsin(x)を設定
    }

    a.SaveAsCsv("cos関数.csv");         // CSVファイルとして保存
    b.SaveAsCsv("sin関数.csv");
}