Tự tính lấy số pi

Thuật toán Gauxơ-La-grăng (Gauss-Legendre)

1. Khởi tạo giá trị
   CODE
   a = 1
   b = 1 / sqrt(2)
   t = 1/4
   x = 1
2. Lặp, cho đến khi chênh lệnh giữa a và b nằm trong sai số cần thiết:
   CODE
   y = a
   a = (a+b) / 2
   b = sqrt(b*y)
   t = t - x * (y-a)^2
   x = 2 * x
3. Pi được tính xấp xỉ như sau:
   CODE
   pi = ((a+b)^2) / (4*t)

The algorithm has second order convergent nature. Then if you want to calculate up to n digits, iteration count of the order log2 n is sufficient. E.g. 19 times for 1 million decimal digits, 31 times for 3.2 billon decimal digits.

Chương trình tính số pi

PiW

Download

SuperPI 1.1e (33,55 triệu chữ số)

Tháng Tám năm 1995, một siêu máy tính tại Đại học Tokyo đã tính pi với 4,294,960,000 chữ số. Chương trình viết bởi D.Takahashi và người cộng sự  Dr. Y.Kanada tại trung tâm máy tính ĐH Tokyo. Kỉ lục này có lẽ được xem là kỉ lục thế giới hiện tại (bài viết tháng 7/97). Chương trình này đã được dịch sang cho các máy tính cá nhân như môi trường Windows 95/NT. Để tính 33,55 triệu chữ số, cần 3 ngày trên máy Pentium 90MHz, 40MB bộ nhớ và 340MB đĩa cứng trống.

Chương trình này miễn phí và bạn cũng được quyền phân phát nó miễn phí.

Download