Найти сумму ряда (pascal)

[loctar]

паскаль составить программу Найти сумму ряда с точностью =10-3 , общий член которого an = ((3^n)*n! )/ (2n)!

Код:

uses crt;
const
eps=0.001;

function fact(n : integer) : longint;
begin
if n <= 1
then  fact := 1
else   fact := n * fact(n - 1);
end;

var
t,s:real;
i:integer;
begin
writeln('*** Alphaeus is thinking... ***');
writeln('***          OK             ***');
writeln;
writeln('Программа вычисляет частичную сумму ряда с заданной точностью');
writeln('  Общий член ряда:  an = ((3^n)*n! )/ (2n)!');
writeln;
s:=0;   i:=1;   t:=1;
while t>eps do
begin
t:=(power(3,i)*fact(i))/(2*fact(2*i));
s:=s+t;
writeln('a(',i,')=',t:0:6,'  S=',s:0:6);
i:=i+1;
end;
writeln();
writeln('Сумма ряда равна',s:0:6);
end.

почему то выдает ошибку unknown identifier, и указывает прямо перед "power" . хотя вроде все объявлено ((

[Serge_Bliznykov]

в TurboPascal нет функции Power()

либо заменяйте возведение в степень через экспоненту, либо пишите функцию возведения в степень путём циклического умножения...

... (2*fact(2*i))

а почему удваиваете факториал?!
в формуле не так...

[Аватар]

Нет необходимости каждый раз пересчитывать общий член ряда. Там рекуррентная зависимость n-го члена от предыдущего:
a(n)=a(n-1)*3/2/(2*n-1)

[loctar]

сделал, заработало, спасибо ...
работает и лааадно...)