if defined PiDebug ( set /a PrintMinQuadIndex=0 ) else ( set /a PrintMinQuadIndex=2 ) set /a PrintMaxQuadIndex = MaxQuadIndex - 1 for /L %%i in (!PrintMinQuadIndex!, 1, !PrintMaxQuadIndex!) do ( set PrintDigit=!%1_%%i! if !PrintDigit! lss 1000 ( if !PrintDigit! lss 100 ( if !PrintDigit! lss 10 ( set PrintDigit=000!PrintDigit! ) else ( set PrintDigit=00!PrintDigit! ) ) else ( set PrintDigit=0!PrintDigit! ) ) set PrintBuffer=!PrintDigit!!PrintBuffer! ) set PrintBuffer=!%1_%MaxQuadIndex%!.!PrintBuffer:x=! echo.%1 = !PrintBuffer! goto :EOF :DivideByInteger if defined PiDebug echo.DivideByInteger %1 %2 set /a DBI_Carry = 0 for /L %%i in (!MaxQuadIndex!, -1, 0) do ( set /a DBI_Digit = DBI_Carry*10000 + %1_%%i set /a DBI_Carry = DBI_Digit %% %2 set /a %1_%%i = DBI_Digit / %2 ) goto :EOF :MultiplyByInteger if defined PiDebug echo.MultiplyByInteger %1 %2 set /a MBI_Carry = 0 for /L %%i in (0, 1, !MaxQuadIndex!) do ( set /a MBI_Digit = %1_%%i * %2 + MBI_Carry set /a %1_%%i = MBI_Digit %% 10000 set /a MBI_Carry = MBI_Digit / 10000 ) goto :EOF :ArctanRecip if defined PiDebug echo.ArctanRecip %1 %2 call :SetToInteger %1 1 call :DivideByInteger %1 %2 call :CopyValue AR_Recip %1 set /a AR_Toggle = -1 set /a AR_K = 3 :ArctanLoop if defined PiDebug ( echo.

goto :EOF ) call :ArctanRecip PiTerm %2 set /a PiEngineFactor=%1 if !PiEngineFactor! lss 0 ( set /a PiEngineFactor *= -1 set Formula=!Formula! call :MultiplyByInteger PiTerm !PiEngineFactor! call :Subtract Pi PiTerm set Operator=- ) else ( call :MultiplyByInteger PiTerm %1 call :Add Pi PiTerm set Operator=+ ) if defined Formula ( set Formula=!Formula! !Operator! !PiEngineFactor!*arctan^(1/%2^) ) else ( set Formula=pi = %1*arctan^(1/%2^) ) shift shift goto PiTermLoop :SetToInteger for /L %%i in (0, 1, !MaxQuadIndex!) do ( set /a %1_%%i = 0 ) set /a %1_!MaxQuadIndex! = %2 goto :EOF :Print set PrintBuffer=x REM Omit a couple of least significant quads, because they will have roundoff errors. call :PiEngine 48 18 32 57 -20 239 call :PiEngine 16 5 -4 239 goto :EOF :PiEngine call :SetToInteger Pi 0 set Formula= :PiTermLoop if '%1' = '' ( call :Print pi echo. Defined talk (echo on) else (echo off) setlocal EnableDelayedExpansion echo.pi.bat - By Don Cross - set /a NumQuads = 30 set /a MaxQuadIndex = NumQuads - 1 echo.

Calculate pi of protein