Число пи ...

[mazaila]

Вот пытался перевестиhttp://ru.wikipedia.org/wiki/Алгоритм_Чудновского На язык паскаль ,но вместо числа пи у меня бог знает что ;
Может поможите мне разобраться с кодом ????

function ymnk(l:integer):integer;
var d,h:integer;
begin
if l=0 then ymnk:=0 else
if l=1 then ymnk:=-1 else
h:=-1;
for d:=1 to l do
h:=h*(-1);
ymnk:=h;

end;


function Factorial(n:word): integer;
var
  fact: integer;
begin
  fact := 1;
  if n=0 then exit;
  for n:=1 to n do
  fact := fact*n;
  Factorial := fact;
end;
var pi,s,n,g:extended; fak1,fak2,fak3,a,b,c,k:integer;

begin
for k:=0 to 10 do begin
a:=6*k;
b:=3*k;
c:=k;
fak1:=Factorial(a);
writeln('fak1:  ',fak1);
fak2:=Factorial(b);
writeln('fak:  ',fak2);
fak3:=Factorial(c);
writeln('fak:  ',fak3);

n:=ymnk(k)*(fak1)*(13591409+545140134*k);
writeln('^:  ',ymnk(k));
writeln('n:=',n);
g:=fak2*(exp(3*ln(fak3)))*exp((3/2+3*k)*ln(640320));
writeln('g:=',g);
s:=n/g;
writeln('s:=',s);
s:=s*12;
writeln('s:=',s);
pi:=1/s;
 writeln('pi',pi)

end;
writeln('pi-full:=',pi);
readln;
end.

___________
Код нужно оформлять по правилам:
тегом [CODE]..[/СODE] (это кнопочка с решёточкой #)
Не забывайте об этом!
Модератор.
Спс ...

[mazaila]

Код небольшой поэтому без коментариев )))

[Sparkman]

Посидел и чего-то наскреб с инета
made in delphi

function Power(X, N : real) : extended;
var t : longint;
    r : real;
    isInteger : boolean;
begin
     if N = 0 then
     begin
          result := 1.0;
          exit;
     end;

     if X = 1.0 then
     begin
          result := 1.0;
          exit;
     end;

     if X = 0.0 then
     begin
          if N > 0.0 then
          begin
               result := 0.0;
               exit;
          end else raise EPowerException.Create('Результат - бесконечность');
     end;

     if (X > 0) then
     try
        result := exp(N * ln(X));
        exit;
     except
        raise EPowerException.Create('Результат - переполнение или потеря значимости');
     end;

{ X - отрицательный, но мы все еще можем вычислить результат, если n целое. }
{ пытаемся получить целую часть n с использованием типа longint, вычисление }
{ четности n не займет много времени }
     try
        t := trunc(n);
        if (n - t) = 0 then
           isInteger := true
        else
            isInteger := False;
     except
        { Лишний бит может вызвать переполнение или потерю значимости }
        r := int(n);
        if (n - r) = 0 then
        begin
             isInteger := true;
             if frac(r/2) = 0.5 then
                t := 1
             else
                t := 2;
        end else
            isInteger := False;
     end;

     if isInteger then
     begin
          {n целое}
          if odd(t) then
             {n нечетное}
             try
                result := -exp(N * ln(-X));
                exit;
             except
                raise EPowerException.Create('Результат - переполнение или потеря значимости');
             end else
                {n четное}
                try
                   result := exp(N * ln(-X));
                   exit;
                except
                   raise EPowerException.Create('Результат - переполнение или потеря значимости');
                end;
     end else
         raise EPowerException.Create('Результат невычисляем');
end;

function Factorial(n:integer): extended;
var fact: extended;
    i:integer;
begin
     fact:=1;
     result := 1;
     if n=0 then exit;
     for i:=1 to n do
     fact := fact*i;
     result := fact;
end;

function GetPi(tochnost : integer) : extended;
var tmp : extended;
    i : integer;
begin
     tmp:=0;
     for i:=0 to tochnost do
          tmp:=tmp+(power((-1),i)*factorial(6*i)*(13591409 + 545140134*i))/(factorial(3*i)*power(factorial(i),3)*power(640320,3*i+3/2));
     tmp := 1/(tmp * 12);
     result:=tmp;
end;

хрень какая то получилась... до tochnost=166 постоянно выдает 3.14159265358979, дальше переполнение.

mazaila, протестил твой код - ругается на "деление на ноль"

[Johnson]

//hint Power(N,M); это стандартная процедура модуля Math;
Позволяет N возвести в степень M

[_-Re@l-_]

Sparkman, ни хрена себе...
А чё, никто не помнит, что в Паскале есть стандартная переменная Pi?
Так-то вот...

Begin
        Writeln(Pi);
        Readln;
End.

[mazaila]

Сообщение от _-Re@l-_

Sparkman, ни хрена себе...
А чё, никто не помнит, что в Паскале есть стандартная переменная Pi?
Так-то вот...

Begin
        Writeln(Pi);
        Readln;
End.

Не, не, не мне надо чтобы теоритически можно было все число пи вычеслить ))))

[mazaila]

Сообщение от Sparkman

Посидел и чего-то наскреб с инета
made in delphi

function Power(X, N : real) : extended;
var t : longint;
    r : real;
    isInteger : boolean;
begin
     if N = 0 then
     begin
          result := 1.0;
          exit;
     end;

     if X = 1.0 then
     begin
          result := 1.0;
          exit;
     end;

     if X = 0.0 then
     begin
          if N > 0.0 then
          begin
               result := 0.0;
               exit;
          end else raise EPowerException.Create('Результат - бесконечность');
     end;

     if (X > 0) then
     try
        result := exp(N * ln(X));
        exit;
     except
        raise EPowerException.Create('Результат - переполнение или потеря значимости');
     end;

{ X - отрицательный, но мы все еще можем вычислить результат, если n целое. }
{ пытаемся получить целую часть n с использованием типа longint, вычисление }
{ четности n не займет много времени }
     try
        t := trunc(n);
        if (n - t) = 0 then
           isInteger := true
        else
            isInteger := False;
     except
        { Лишний бит может вызвать переполнение или потерю значимости }
        r := int(n);
        if (n - r) = 0 then
        begin
             isInteger := true;
             if frac(r/2) = 0.5 then
                t := 1
             else
                t := 2;
        end else
            isInteger := False;
     end;

     if isInteger then
     begin
          {n целое}
          if odd(t) then
             {n нечетное}
             try
                result := -exp(N * ln(-X));
                exit;
             except
                raise EPowerException.Create('Результат - переполнение или потеря значимости');
             end else
                {n четное}
                try
                   result := exp(N * ln(-X));
                   exit;
                except
                   raise EPowerException.Create('Результат - переполнение или потеря значимости');
                end;
     end else
         raise EPowerException.Create('Результат невычисляем');
end;

function Factorial(n:integer): extended;
var fact: extended;
    i:integer;
begin
     fact:=1;
     result := 1;
     if n=0 then exit;
     for i:=1 to n do
     fact := fact*i;
     result := fact;
end;

function GetPi(tochnost : integer) : extended;
var tmp : extended;
    i : integer;
begin
     tmp:=0;
     for i:=0 to tochnost do
          tmp:=tmp+(power((-1),i)*factorial(6*i)*(13591409 + 545140134*i))/(factorial(3*i)*power(factorial(i),3)*power(640320,3*i+3/2));
     tmp := 1/(tmp * 12);
     result:=tmp;
end;

хрень какая то получилась... до tochnost=166 постоянно выдает 3.14159265358979, дальше переполнение.

mazaila, протестил твой код - ругается на "деление на ноль"

А у меня 'Работает'...

[Anatole]

Не, не, не мне надо чтобы теоритически можно было все число пи вычеслить ))))

Теоретически число пи бесконечно. поэтому "всё" не вычислишь даже теоретически

[mazaila]

[QUOTE=Anatole;763815]Теоретически число пи бесконечно. поэтому "всё" вычислишь даже теоретически[/QUOT

Нет нечего бесконечного )))

[Sparkman]

Сообщение от mazaila

Нет нечего бесконечного )))

Тогда назовите мне, пожалуйста, самое максимальное число