@echo off&title 计算圆周率 by plp626@cn-dos.net

setlocal enabledelayedexpansion

if not %1.==. (set c=%1) else set c=100

set/a c=(c*100/3)+70,cc=c/10,count=0

for /l %%a in (1 1 %cc%)do set/a f_%%a=2000

for /l %%a in (%c% -132 100)do (set/a n=%%a/10,m=2*n-1

set/a "d=f_!n!*10000,f_!b!=d%%m,d=d/m,n-=1"

for /l %%b in (!n! -1 1)do (set/a n=%%b,m=2*n-1

   set/a "d=d*n+f_!n!*10000,f_!n!=d%%m,d=d/m,n-=1"

)

set/a "an=e+d/10000,e=d%%10000"

if !an! lss 1000 set an=000!an!

set/p=!an:~-4!<nul

)

echo\&echo/&set/p=数学爱好者加我QQ:************，

Originally posted by Climbing at 2008-4-30 02:11 PM:
其实，没有所谓的纯批处理，所谓的纯只是限定用系统内置的命令（包括外部命令）。所以，只要将相关的外部命令（例如汇编程序）用批处理写入磁 ...

Atn 函数返回数值的反正切值。

Atn(number)

number 参数可以是任意有效的数值表达式。

说明
Atn 函数计算直角三角形两个边的比值 (number) 并返回对应角的弧度值。此比值是该角对边的长度与邻边长度之比。 结果的范围是从 -pi/2 到 pi/2 弧度。

弧度变换为角度的方法是将弧度乘以 pi/180。反之，角度变换为弧度的方法是将角度乘以180/pi 。

下面的示例利用 Atn 来计算 pi 的值:

Dim pipi = 4 * Atn(1) ' 计算 pi 的值。注意 Atn 是 Tan（将角作为参数返回直角三角形两边的比值）的反三角函数。不要混淆 Atn 与余切（正切的倒数 (1/tangent)）函数。

k = 1: s = 1: t = 1: m = 1

time_ = timer

While Abs(t) > 0.0000001

k = k + 2  

m = -m  

t = m / k  

s = s + t  

Wend

tt = timer - time_

MsgBox tt & vbCrLf & "圆周率π=" & 4 * s

#include < stdio.h>

long a=10000, b, c=2800, d, e, f, g; 

main(){

for(;b-c;)f=a/5; 

for(;d=0,g=c*2;c-=14,printf("%.4d",e+d/a),e=d%a) 

for(b=c;d+=f*a,f=d%--g,d/=g--,--b;d*=b);

} 

#include < stdio.h>

long b, c=2800, d, e, f; 

main(){

    while(b-c!=0){

        f=2000;

        b++;

        }



    while(c!=0){

        b=c;

        d=f*10000;

        f=d%(b*2-1);

        d=d/(b*2-1);

        --b;

        while(b!=0){

            d=d*b+f*10000;

            f=d%(b*2-1);

            d=d/(b*2-1);

            --b;

            }

        c-=14;

        printf("%.4d",e+d/10000);

        e=d%10000;

        }

    } 

Originally posted by slore at 2008-4-30 04:35 PM:
C 语言的：
这个强悍，一瞬间（不到1秒）输出2800位：


代码不知道你调试了没？反正看着是有错误。代码输出的是余数？貌似全是0
for的那部分有 ...

n=100

for i = 1 to n

  str=2+(n+1-i)/((n+1-i)*2+1)*str

next

msgbox str

@echo off

setlocal EnableDelayedExpansion

set str=1

for /l %%i in (100,1,1) do set/a str=2+%%i/(%%i*2+1)*!str!

echo %str%

pause

Originally posted by zh159 at 2008-4-30 05:11 PM:
这个VBS的很快，不过受VBS限制位数

n=100

for i = 1 to n

  str=2+(n+1-i)/((n+1-i)*2+1)*str

next

msgbox str

如果转为批处理则为：

@echo off

setlocal E ... 

'将那个C的代码用VBS处理了下,速度很快,因为是
'每位用一个数处理的,所以没有位数限制。
'恩,记得用CScript调用……

'这个是连续运算最后输出,比一行一行的输出速度快。

Dim b, c, d, e,n, f()
n = 1000 '在这里定义位数
c = n / 4 * 14
ReDim f(c)
For b = 0 To c
f(b) = 2000
Next

Do While(c)
b = c
d = f(b) * 10000
f(b) = d Mod (b * 2 - 1)
d = d \ (b * 2 - 1)
b = b - 1
Do While(b)
d = d * b + f(b) * 10000
f(b) = d Mod (b * 2 - 1)
d = d \ (b * 2 - 1)
b = b - 1
Loop
c = c - 14
i = e + d \ 10000
i = Right("0000" & i,4)
If k < 10 Then
k = k + 1
ShowLine = ShowLine & i
Else
k = 0
'WSH.Echo ShowLine
ShowLine = ShowLine & vbCrLf
End If
e = d Mod 10000
Loop
WSH.Echo ShowLine

Originally posted by xtanbmy at 2008-5-1 03:52:
大虾们，能不能把结果输出到一个文本文件中呀？

355/113=3.1415926535897...

Originally posted by PPdos at 2008-5-1 08:54:
古算术法

Originally posted by PPdos at 2008-5-1 08:54:
古算术法