Помогите с Pascal'ем

Суббота, 28 Октября 2006 г. 23:38 + в цитатник

Хренотень все записи автора Хелп, не могу найти ошибку.
Задание такое:
Найти рекуррентную формулу, используя которую для заданного с клавиатуры значения m и x (2>x>1), найдите сумму S:
S={*тут следует значок суммы от n=1 до m*}(-1)ⁿ*(x-1)²ⁿ/(2n)!!
Я делала

program summa;
var f, m, n, k, i : integer;
x, s : real;
begin
repeat
writeln('vvedi m,x');
readln(m);
readln(x);
until ((x>1) and (x<2));
s:=-(x-1)/2;
writeln('s1=',s);

for n:=1 to m do
begin
f:=1;
k:=2*n;
for i:=1 to k do begin
f:=f*i;
end;
writeln(f);

if (n mod 2)=0 then s:=s+(exp(2*n*ln(x-1)))/f
else s:=s-(exp(2*n*ln(x-1)))/f;

end;
writeln('s=', s:10:4);
readln;
end.

Но во-1, даже при маленьких m, факториал начинает зашкаливать и превращаться в ноль. А во-2, СУММА ВСЁ ВРЕМЯ ПОЛУЧАЕТСЯ ОДИНАКОВАЯ!
Что делать? :(
Понимаю, что наверна это жутко просто, но я тупость невиданная, ибо после гумманитарного класса без информатики ВААПЩЕ поступила в техн. фуз на программирование :)
Хелп!Мне в понедельник сдавать это и мноооогое другое...

Настроение сейчас - мозг вышел погулять

Рубрики:

вопросы

Комментировать

« Пред. запись — К дневнику — След. запись »

Страницы: [1] [Новые]

Никнахуйненужен обратиться по имени Воскресенье, 29 Октября 2006 г. 00:48 (ссылка)

Для переменной, которая используется значение факториала лучше определять тип вещественный (real), т.к. диаппазон integer где-то от -32000 до 32000 (точно не помню, но это и не важно), даже 8!=40320

TLС v.0.7.1

Ответить С цитатой В цитатник

Wethril_Haldar обратиться по имени Воскресенье, 29 Октября 2006 г. 01:19 (ссылка)

1-е - двойной факториал - это не факториал факториала а факториал касающийся перемножения только нечетных членов последовательности
1*3*5*...*(2n-3)*(2n-1)

2-e факториал можно оформить как рекуррентную(т.е. вызывающую саму себя) функцию. так будет работать побыстрей, и красиво это , в конце-концов)))

Ответить С цитатой В цитатник

Dim_K обратиться по имени Помогите с Pascal'ем Воскресенье, 29 Октября 2006 г. 08:38 (ссылка)

Что я могу сказать... Несколько неточностей:
1). Вы написали совсем не рекуррентную формулу. Что такое рекуррентная формула? Вот, например, числа Фибоначчи x[1]=1, x[2]=1, x[i]=x[i-2]+x[i-1], i больше 2, то есть следующее слагаемое вычисляется при помощи предыдущих.
2). Выше уже было сказано про двойной факториал, может там все-таки обычный?

program ex;
var m, n, i : integer;
x, s_pred,s_next,sum : real;

begin
repeat
writeln('vvedi m,x');
readln(m);
readln(x);
until ((x>1) and (x<2) and (m>0));
sum:=-(x-1)*(x-1)/2;
s_pred:=sum;
if m>1 then
for n:=2 to m do
begin
s_next:=-s_pred*(x-1)*(x-1)/(2*n); {если нужен обычный факториал а не двойной то делить не на 2n а на n}
sum:=sum+s_next;
s_pred:=s_next;
end; { for n:=2 to m}
writeln('s=', sum:10:4);
readln;
end.
При рекурретных вычислениях не возникает необходимости заводить пременный для "больших" вычислений (факториал), в сумме они обычно погашаются другими множетелями или слагаемыми

LI 5.09.15

Ответить С цитатой В цитатник

Dim_K обратиться по имени Воскресенье, 29 Октября 2006 г. 08:47 (ссылка)

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

Ответить С цитатой В цитатник

code:

program summa;
type my_type=1..999999999;
var j,m: integer;
x,sum: real;
function st(a:real;n:byte):real;
var s:real;
i:byte;
begin
s:=1;
for i:= 1 to n do
s:=s*a;
st:=s;
end;
function fak(n:byte):my_type;
begin
if n=1 then fak:=1
else fak:=n*fak(n-2);
end;
begin
repeat
writeln('vvedi m,x');
readln(m,x);
until ((x>1) and (x<2));
sum:=0;
for j:=1 to m do
begin
sum:=sum+st(-1,j)*st(x-1,2*j)/fak(2*j-1);
end;
writeln('summa=', sum:10:4);
readln;
end.

LiveInternetLiveInternet

-Рубрики

-Подписка по e-mail

-Поиск по дневнику

-Интересы

-Постоянные читатели

-Статистика

Помогите с Pascal'ем