[数值计算] 批处理计算圆周率任意位

plp626

直接双击:默认计算100位,
命令行下:第一个参数(正整数)为计算圆周率的位数

 
@echo off&title 计算圆周率 by plp626
setlocal enabledelayedexpansion
if %1.==. (set c=100) else set/a c=%1
set/a c=(c*100/3)+70,cc=c/10,n=cc
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_!n!=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:275928264，验证号:plp626
复制代码

所用公式:
pi=2*1+2*1!/3!!+2*2!/5!!+2*3!/7!!+...+2*k!/(2*k+1)!!+...

yslyxqysl

最后少了个<nul&pause.

[ 本帖最后由 yslyxqysl 于 2009-2-22 21:58 编辑 ]

plp626

no,no,no
不需要....

lockmove

有 set/p=  就停了吧.
pi=2*1+2*1!/3!!+2*2!/5!!+2*3!/7!!+...+2*k!/(2*k+1)!!+...

这公式在后面越来越慢,计算也有限,不过挺好了.

batman

plp626，兄弟，我们有好久没见了，我是cn-dos的ZW750516，你能来到我们论坛，我真是太高兴了。。。

plp626

zw兄，还记得俺，感觉很亲切，其实我早来这里了，只是没时间上网罢了，

BBCC

昨晚正想问你公式是啥的...
之前见est用python写过一个,好像速度更快...

Batcher

python的优势之一是它的数据类型可以存放N大的数据^_^

zqz0012005

plp626这个“老古董”居然重新冒出来了。

523066680

核潜水艇浮起来拉 ～
话说这里成为迎宾大道了  欢迎～

plp626

一楼代码有一处笔误，导致后面精度有所丢失，两年多了，今儿才发现，现在已改正过来。

CrLf

这个牛，顶一下古董兄

zz100001

为什么叫“古董兄”？

▄︻┻═┳一

13# zz100001


你看他的注册日期就明白了。

applba

关键是效率……