Расчёт числа Пи по формуле Уоллиса

[egoregor]

/**********************************************************
Расчёт числа  ПИ = 3.141 по формуле Уоллиса ( 1616 - 1703 )

PI     2^2     4^2     6^2        (2 * n) ^2
--- = ----- * ----- * ----- * ------------------
 2     1*3     3*5     5*7    (2n - 1) * ( 2n + 1)

**********************************************************/
using System;

class __exe {
public static void Main () { Console.Clear();

double Pi = 3.14159265358979323846; // 26433832795;
double PI = PI_WALLES () ;
double pi = pi_walles () ;

Console.WriteLine("{0} {1}" ,Pi ,"= windows calc.exe" );
Console.WriteLine("{0} {1}" ,PI ,"= PI_WALLES       " );
Console.WriteLine("{0} {1}" ,pi ,"= pi_walles       " );

} // Main (end)
/*********************************************************/
public static double PI_WALLES () {
double x = 0;
return (double) (x*4) ; }
/*********************************************************/
public static double pi_walles ( ) {
double n    = 1 ;         // Число  n  из формулы Уоллиса
double pi   = 1 ; 
double lim  = 123456789 ; // Предел для счётчика for()
double chis = 0 ;         // Числитель дроби.   Делимое
double znam = 0 ;         // Знаменатель дроби. Делитель
double r    = 0 ;         // Результат деления дроби

for ( n = 1 ; n <= lim ; n ++ ) {
chis = znam = 0 ;
chis = ( 2*n) ;
znam = ((2*n) - 1) * 
       ((2*n) + 1) ;
chis = chis * chis;
  r  = chis / znam;
  pi = pi * r;                  } // for (end)

return (double) (pi*2); }

} // __exe (end)
/************************** END **************************/

[Shad0wF1rst]

А зачем ты 3 раза умножаешь 2 * n когда у тебя уже готовй результат есть в переменной chis. И вообще к чему это программа? Она не работает совсвем, или работает но не так или не компилируется?

И еще просто можно формулу чутка упростить, тем самым сократив количество действий:

public static double pi_walles ( ) {
int n    = 1 ;         // Число  n  из формулы Уоллиса
double pi   = 1 ; 
double lim  = 123456789 ; // Предел для счётчика for()

for ( ; n <= lim ; n ++ ) pi *= (1 + 1.0/((n*n << 2) - 1));
return (double) (pi*2); }

}