Программа для вычисления выражения по формуле (вычисление числового ряда)

[Vitalyas]

Привет всем! Задача такая написать программу для решения этого уравнения (n! - факториал)
Написал программу, но при вводе в программу n = 2 или более ответ расходится от просчитанного в ручную, но при вводе n = 1 ответ правильный...Нужна помощь

Вот текст программы:

uses crt;
var a, b, c, x, s, f:real;
    n, i, j:integer;
   begin
        write('Введите х. x=');
        read(x);
        write('Введите n. n=');
        read(n);
        if n=-0.5
           then WriteLn('Функция не определена')
                else
                    For i:=1 to n do
                       begin
                          If odd(n) then a:=-1
                                    else a:=1;
                          f:=1;
                          for j:=1 to n do
                             begin
                             f:=f*j
                             end;
                          b:=1/f;
                          c:=exp((2*n+1)*ln(x))/(2*n+1);
                          s:=x+a*b*c;
                       end;
                    WriteLn('Сумма=', s:8:2);
   end.

[Stilet]

Сделай ты функциями.

function fac(n:integer):integer;var x:integer;
 begin x:=n; dec(n); repeat x:=x*n; until n=1; fac:=x; end;

И выражение:

double y(x:double;n:integer); var a:integer;
begin
 If odd(n) then a:=-1   else a:=1;end;
 y:=a*exp((2*n+1)*ln(x))/(fac(n)*(2*n+1));

И в цикле уже юзай функцию y();

[Vitalyas]

Извините, но я паскаль начал изучать совсем недавно. Сможете написать код программы, если вам не сложно?

[Serge_Bliznykov]

а я бы воспользовался реккурентой формулой (очевидно, что для нахождения K-го элемента нужно числитель R-1 члена ряда умножать на -X^2, а его знаменатель на K ):
ну и делить на (2*K + 1).

Примерно так

        write('Введите х. x=');
        read(x);
        write('Введите n. n=');
        read(n);
        sum := x;
        f := x;
        for i:=1 to n do begin
             f := -f *x*x/i;
             sum := sum + f / (2*i+1);
        end;
       WriteLn('Сумма=', s:8:2);
end.

не проверял...

[Vitalyas]

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

а я бы воспользовался реккурентой формулой (очевидно, что для нахождения K-го элемента нужно числитель R-1 члена ряда умножать на -X^2, а его знаменатель на K ):
ну и делить на (2*K + 1).

Примерно так

        write('Введите х. x=');
        read(x);
        write('Введите n. n=');
        read(n);
        sum := x;
        f := x;
        for i:=1 to n do begin
             f := -f *x*x/i;
             sum := sum + f / (2*i+1);
        end;
       WriteLn('Сумма=', s:8:2);
end.

не проверял...

Что-то я как-то не понял...Подставил :

        sum := x;
        f := x;
        for i:=1 to n do begin
             f := -f *x*x/i;
             sum := sum + f / (2*i+1);
        end;

заместо С в моем коде, но ответ совсем отличается...Может я что-то не так сделал? Вот код:

uses crt;
var a, b, c, x, s, f, sum:real;
    n, i, j:integer;
   begin
        write('Введите х. x=');
        read(x);
        write('Введите n. n=');
        read(n);
        if n=-0.5
           then WriteLn('Функция не определена')
                else
                    For i:=1 to n do
                       begin
                          If odd(n) then a:=-1
                                    else a:=1;
                          sum:=x;
                          f := x;
                          for j:=1 to n do 
                             begin
                                f := -f *x*x/j;
                                sum := sum + f / (2*j+1);
                             end;
                          b:=1/f;
                          s:=x+a*b*sum;
                       end;
                    WriteLn('Сумма=', s:8:2);
   end.

Попробовал по другому, убрав цикл с j, вот что получилось, но ответ все равно не сходится:

uses crt;
var a, b, c, x, s, f, sum:real;
    n, i, j:integer;
   begin
        write('Введите х. x=');
        read(x);
        write('Введите n. n=');
        read(n);
        if n=-0.5
           then WriteLn('Функция не определена')
                else
                    For i:=1 to n do
                       begin
                          If odd(n) then a:=-1
                                    else a:=1;
                          sum:=x;
                          f := x;
                          begin
                                f := -f *x*x/i;
                                sum := sum + f / (2*i+1);
                          end;
                          b:=1/f;
                          s:=x+a*b*sum;
                       end;
                    WriteLn('Сумма=', s:8:2);
   end.

[Serge_Bliznykov]

я дал ПОЛНОСТЬЮ готовый код!


проверять нужно такой код:

uses crt;
var
  x, f, sum: real;
  n, i: integer;
begin
  write('Введите х. x=');
  readLn(x);
  write('Введите n. n=');
  readLn(n);

  sum := x;
  f := x;
  for i := 1 to n do begin
    f := -f * x * x / i;
    sum := sum + f / (2 * i + 1);
  end;
  WriteLn('Сумма=', sum: 8: 2);
  Readln
end.

[Vitalyas]

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

я дал ПОЛНОСТЬЮ готовый код!


проверять нужно такой код:

uses crt;
var
  x, f, sum: real;
  n, i: integer;
begin
  write('Введите х. x=');
  readLn(x);
  write('Введите n. n=');
  readLn(n);

  sum := x;
  f := x;
  for i := 1 to n do begin
    f := -f * x * x / i;
    sum := sum + f / (2 * i + 1);
  end;
  WriteLn('Сумма=', sum: 8: 2);
  Readln
end.

Спасибо большое! Код правильный! Я просто ваш код пробовал до этого и что-то ответ не сошелся, а вот сейчас вы полностью расписали и сошелся. Спасибо большое еще раз!

[Serge_Bliznykov]

Да всегда пожалуйста! Обращайтесь!

[Vitalyas]

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

Да всегда пожалуйста! Обращайтесь!

А можно вопросик один по этой же функции?
Возможно ли решить эту функцию так как записано в код в самом начале, то есть n! расписав через цикл?

[Serge_Bliznykov]

Ну да, разумеется можно!!
Просто получится сложнее и менее эффективно (для каждого члена ряда заново будут просчитываться степени числа X и факториал N), но, тем не менее, можно и так решить задачу.
Если всё правильно расписать, то результат должен получиться один и тот же!