using Winston

summand = zeros(11)

# หาค่าแต่ละพจน์ของ summation ก่อน
for i=0:10
    summand[i+1] += 1/factorial(i)
end

# ค่า error ของผลบวกถึงพจน์ที่ n เทียบกับค่า e
err = zeros(10)
for i=1:10
    err[i] = abs(e - sum(summand[1:i]))
end

plot(1:10,err, "b")
xlabel("N")
ylabel("Error")
title("Numerical computation of Euler Number")