#include      <wfl.h>
using namespace wfl;

void main()
{
    Start();
    WaveField a(256, 256, 2e-6);

    //===== ‰~Œ`ŠJŒû‚ða‚É�¶�¬‚·‚é
    int i, j;
    for (i = 0; i < a.Nx(); i++)
    {
        for (j = 0; j < a.Ny(); j++)
        {
            double x = a.itox(i);
            double y = a.jtoy(j);
            double r = sqrt(x*x + y*y);
            if (r < 0.1e-3)							//”¼Œa0.1mm‚̉~Œ`ŠJŒû
                a.SetPixel(i, j, Complex(1.0, 0));	//ŠJŒû“à‚Í1.0
            else
                a.SetPixel(i, j, Complex(0, 0));	//ŠJŒûŠO‚Í0.0

        }
    }
    a.SaveAsGrayBmp("Aperture.bmp", AMPLITUDE);		//ŠJŒûŒ`�ó‚Ì•Û‘¶

    //===== ‘шæ�§ŒÀŠpƒXƒyƒNƒgƒ‹–@‚É‚æ‚é‰ñ�Ü“`”ÀŒvŽZ
    a.AsmProp(1e-3);								//1mm‚Ì“`”À‰ñ�ÜŒvŽZ‚ðŽÀ�s
    a.Normalize();
    a.SaveAsGrayBmp("Diffraction.bmp", INTENSITY);	//‰ñ�Ü‘œ(�U•�‘œ)‚Ì•Û‘¶
}