scilab source for the finite Fourier Integral

// 方形パルスのフーリエ変換説明
// Kreyszig 2.8節・例２ (p.96-98) 図2.22に対応
// p.99の2行目の式を数値計算、プロット
x = -3:0.05:3;
y = zeros(1,length(x))

a =[8 16 32 64]
for p=1:4
    for i=1:length(x)
        result = intg(0,a(p)*(x(i)+1),sinc)    // sinc(x) = sin(x)/x
        result2 = intg(0,a(p)*(x(i)-1),sinc)
        y(i) = (result(1)-result2(1))/%pi
    end
    subplot(2,2,p)
    plot(x,y,"linewidth",3)
    title('a='+string(a(p)), 'fontsize',4);
end