Вычисление функции (Pascal)

[lloudb]

Задача:

Код программы:

function f(x: integer): integer;
	var i,s: integer;
	begin
		s:=1;
		for i:=1 to x do s:=s*i;
		f:=s;
	end;
function s(x: real): real;
	var 
		k: integer;
		a,b: real;
	begin
		a:=0; b:=0;
		for k:=1 to 10 do a:=a+((exp(ln(x) * (2*k+1)))/f(2*k+1));
		for k:=1 to 10 do b:=b+((exp(ln(x) * (2*k)))/f(2*k));
		s:=a / b;
	end;
var x: real;
begin
	read(x);
	writeln(s(x));
end.

Вывод:

Runtime error 200 at $000000000040027D
  $000000000040027D
  $0000000000400367
  $0000000000400180

В чем ошибка? Как читать такие Runtime Error?

[type_Oleg]

Сначала посмотрите, что такое Integer:
-2147483648..2147483647
А 21! = 5.10909E+19

В Integer влезет максимум факториал 12

А вообще - лучше последующее слагаемое вычислять через предыдущее, не вычисляя факториалы и степени.

[lloudb]

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

Сначала посмотрите, что такое Integer:
-2147483648..2147483647
А 21! = 5.10909E+19

В Integer влезет максимум факториал 12

А вообще - лучше последующее слагаемое вычислять через предыдущее, не вычисляя факториалы и степени.

Поменял на тип real - ошибки те же.

[type_Oleg]

Вот так наверное:

program Pr_lloudb;

const NN=10;

function fun1(x:Real):Real;
var t,s:Real;
    i:Integer;
begin
 t:=x*x*x/(2*3);
 s:=t;
 for i:=1 to NN-1 do
  begin
   t:=t*x*x/((2*i+2)*(2*i+3));
   s:=s+t;
  end;
 fun1:=s;
end;
function fun2(x:Real):Real;
var t,s:Real;
    i:Integer;
begin
 t:=x*x/2;
 s:=t;
 for i:=1 to NN-1 do
  begin
   t:=t*x*x/((2*i+1)*(2*i+2));
   s:=s+t;
  end;
 fun2:=s;
end;
function fun3(x:Real):Real;
begin
 fun3:=fun1(x)/fun2(x);
end;

var x:Real;
begin
 ReadLn(x);
 WriteLn(fun3(x));
 ReadLn;
end.

[lloudb]

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

Вот так наверное:

program Pr_lloudb;

const NN=10;

function fun1(x:Real):Real;
var t,s:Real;
    i:Integer;
begin
 t:=x*x*x/(2*3);
 s:=t;
 for i:=1 to NN-1 do
  begin
   t:=t*x*x/((2*i+2)*(2*i+3));
   s:=s+t;
  end;
 fun1:=s;
end;
function fun2(x:Real):Real;
var t,s:Real;
    i:Integer;
begin
 t:=x*x/2;
 s:=t;
 for i:=1 to NN-1 do
  begin
   t:=t*x*x/((2*i+1)*(2*i+2));
   s:=s+t;
  end;
 fun2:=s;
end;
function fun3(x:Real):Real;
begin
 fun3:=fun1(x)/fun2(x);
end;

var x:Real;
begin
 ReadLn(x);
 WriteLn(fun3(x));
 ReadLn;
end.

В выражение присутствует степень и факториал. Где вы вычисляете их?

[type_Oleg]

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

В выражение присутствует степень и факториал. Где вы вычисляете их?

Нигде. Я вычисляю последующее слагаемое через предыдущее.

[lloudb]

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

Нигде. Я вычисляю последующее слагаемое через предыдущее.

Не понимаю)
Объясните пожалуйста, как Вы решили

[Аватар]

f(k) = x^(2k+1)/(2k+1)!
f(k-1) = x^(2k-1)/(2k-1)!
f(k)/f(k-1) = x^(2k+1)*(2k-1)!/((2k+1)!*x^(2k-1)) = x^2/2k/(2k+1) -->

f(k) = f(k-1)*x^2/2k/(2k+1)