Использование подпрограммы функции и факториала.

[ShaoKhan]

Добрый вечер.
Наткнулся на такую вот задачу!
Дано: y= 1-(1/1!)+(1/2!)-(1/3!)+...+ - (1/n!)
И вся эта сумма ряда должна свестись к значению (1/e), что приблизительно равно 0.37373737...

Вот мой код программы:

program Project2;

{$APPTYPE CONSOLE}

uses
  SysUtils,Math;
  var y,c:real;n,j,k,s,p:integer;
  function sum(var n1:integer):real;
  var i:integer;
  c1:real;
begin
c1:=1;
for i:=1 to n1 do begin
s:=1;                               //
for p:=1 to n1 do begin       //это факториал
s:=s*p;                            //
c1:=(c1+(Power(-1,i))*(1/s)); //это сама функция
end;
end;
sum:=c1;
end;
begin
writeln('Enter n');
readln(n);
for j:=1 to n do begin
y:=sum(j);
writeln('y=',y:2:10);
readln;
end
end.

К сожалению оно не сходится к этому числу.
Подскажите, чем можете!

[Juffin]

код хотя бы отформатировать надо, чтобы можно было разобраться.

program Project2;

{$APPTYPE CONSOLE}

uses SysUtils,Math;

var y,c:real;
     n,j,k,s,p:integer;

function sum(var n1:integer):real;
var i:integer;
     c1:real;
begin

c1:=1;

for i:=1 to n1
do begin
   s:=1;                               
   for p:=1 to n1
   do begin
      s:=s*p;                            
      c1:=(c1+(Power(-1,i))*(1/s));
      end;
   end;

sum:=c1;
end;

begin
writeln('Enter n');
readln(n);
for j:=1 to n 
do begin
   y:=sum(j);
   writeln('y=',y:2:10);
   readln;
   end
end.

Решение, кстати, весьма красивое. Вот только непонятно, зачем внешний цикл. Возводи -1 в степень p и все.

[VIK_aka_TOR]

держи... подправил... немного по другому реализовал... в процедуру запихнул само вычисление обратного факториала...

program Project2;
  var y:real;
  n,j:integer;
  
function factor(n1:integer):real;
var i:integer;
  c1:real;
begin
    c1 := 1;
    for i:=1 to n1 do // cчитаем факториал N!
      c1:= c1 * i;
    if odd(n1) then // если нечетный то
      factor := -1/c1
    else factor := 1/c1;  
end;
begin
  writeln('Enter n');
  readln(n);
  y := 1;
  for j:=1 to n do
    begin
      y := y + factor(j);
    end;
writeln('y=',y:2:10);
readln;
end.

[ShaoKhan]

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

Решение, кстати, весьма красивое. Вот только непонятно, зачем внешний цикл. Возводи -1 в степень p и все.

Если это делать, то ряд только идёт на уменьшение, а не сходится к числу.

[ShaoKhan]

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

держи... подправил... немного по другому реализовал... в процедуру запихнул само вычисление обратного факториала...

Спасибо!
Но вот хотелось бы, что бы он пошагово высчитывал всё, а не сразу выдавал ответ.

[VIK_aka_TOR]

пошагово... это же лишь в цикл запихнуть... и все... )))
к примеру как было у вас...по нажатию Enter ))

begin
  writeln('Enter n');
  readln(n);
  y := 1;
  for j:=1 to n do
  begin
    begin
      y := y + factor(j);
    end;
writeln('y=',y:2:10);
readln;
end;
readln;
end.

[ShaoKhan]

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

пошагово... это же лишь в цикл запихнуть... и все... )))
к примеру как было у вас...по нажатию Enter ))

Как говориться, смотришь в книгу - видишь фигу, так и у меня сейчас.
Не могу понять куда засунуть цикл, и от чего до чего брать..

[VIK_aka_TOR]

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

Как говориться, смотришь в книгу - видишь фигу, так и у меня сейчас.
Не могу понять куда засунуть цикл, и от чего до чего брать..

собственно в предыдущем посте подредактировал... для Enter... можно даже видоизменить и сделать так... будет попроще и красивше...

begin
  writeln('Enter n');
  readln(n);
  y := 1;
  for j:=1 to n do
  begin
   y := y + factor(j); // подсчет на данном шаге
   writeln('y=',y:2:10); // вывод на данном шаге
   readln; // ожидание нажатия Enter
  end;
readln;
end.

[ShaoKhan]

Всё!
Понял, разобрался!
Огромное спасибо!