Вычисление последовательностей, помогите плиз Опции

[Beauty_]

Даны числовой ряд и некоторое число E. найти сумму тех членов ряда, модуль которых больше или равен заданному E. общий член ряда имеет вид

М

Beauty, пожалуйста выбирай более экономные форматы. Я отредактировал (сравни: твой .bmp был 92К, этот .gif - всего лишь 1.4К).

[kess]

как решить такую задачу, со степенью и факториалом в одной формуле?

[Lapp]

как решить такую задачу, со степенью и факториалом в одной формуле?

Сделай функции для вычисления всего этого.
Например:

function Factorial(n:integer):LongInt;
begin
  if n>1 then Factorial:=n*Factorial(n-1) else Factorial:=1
end;

function IntPower(a:real; n:integer):real;
begin
  if n>0 then Power:=a*Power(a,n-1) else Power:=1
end;

Дальше понятно?

[kess]

Дальше понятно?

не оч (((

[Lapp]

Хорошо, продвигаемся вместе..
Вот примерно так это можно сделать (с учетом функций, приведенных выше)

S:=0;
n:=1;
repeat
  a:=Power(2/n,n)*Factorial(n);
  if a>=e then S:=S+a else break;
  Inc(n)
until false;

Теперь стало понятнее?


Добавил увеличение n - забыл в первый раз, извините

Сообщение отредактировано: Lapp - 2.05.2007 8:23

[kess]

Хорошо, продвигаемся вместе..
Вот примерно так это можно сделать (с учетом функций, приведенных выше)

S:=0;
n:=1;
repeat
  a:=Power(2/n,n)*Factorial(n);
  if a>=e then S:=S+a else break
until false;

Теперь стало понятнее?

ну так

[Lapp]

ну так

?
что это значит?

[kess]

?
что это значит?

я уже мало что соображаю

[Lapp]

я уже мало что соображаю

Скорее еще, а не уже..
Написать всю прогу?

[kess]

Скорее еще, а не уже..
Написать всю прогу?

если не влом, то да))) буду очень благодарна)

[Lapp]

Вот так - вроде должно работать (не проверял)

function Factorial(n:integer):LongInt;
begin
  if n>1 then Factorial:=n*Factorial(n-1) else Factorial:=1
end;

function IntPower(a:real; n:integer):real;
begin
  if n>0 then IntPower:=a*IntPower(a,n-1) else IntPower:=1
end;

var
  a,S,e:real;
  n:integer;

begin
  Write('Введите точность E: ');
  ReadLn(e);
  S:=0;
  n:=1;
  repeat
    a:=IntPower(2/n,n)*Factorial(n);
    if a>=e then S:=S+a else break;
    Inc(n)
  until false;
  WriteLn('Сумма ',n,' членов ряда (меньших ',e,') равна ',S)
end.

Я в предыдущем фрагменте забыл увеличивать n . Сейчас испрвлю там тоже..

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

Сообщение отредактировано: Lapp - 2.05.2007 8:34

[kess]

Вот так - вроде должно работать (не проверял)

function Factorial(n:integer):LongInt;
begin
  if n>1 then Factorial:=n*Factorial(n-1) else Factorial:=1
end;

function IntPower(a:real; n:integer):real;
begin
  if n>0 then Power:=a*Power(a,n-1) else Power:=1
end;

var
  a,S,e:real;
  n:integer;

begin
  Write('Введите точность E: ');
  ReadLn(e);
  S:=0;
  n:=1;
  repeat
    a:=Power(2/n,n)*Factorial(n);
    if a>=e then S:=S+a else break;
    Inc(n)
  until false;
  WriteLn('Сумма ',n,' членов ряда (меньших ',e,') равна ',S)
end.

Я в предыдущем фрагменте забыл увеличивать n . Сейчас испрвлю там тоже..

спасибо

Сообщение отредактировано: kess - 2.05.2007 8:43

[Lapp]

Но вообще-то это не лучший способ считать, поскольку числитель и знаменатель очень большие и переполняются, хотя их частное вполне нормальное.. Поэтому количество членов ряда сильно ограничено.
Вот так будет лучше:

var
  a,S,e:real;
  i,n:integer;

begin
  Write('Введите точность E: ');
  ReadLn(e);
  S:=0;
  n:=1;
  repeat
    a:=1;
    for i:=1 to n do a:=a*2*i/n;
    if a>=e then S:=S+a else break;
    Inc(n)
  until false;
  WriteLn('Сумма ',n,' членов ряда (меньших ',e,') равна ',S)
end.

[-Volvo-]

Вот так будет лучше:

Возможно, только это совсем не то же самое, что требуется по условию... Условие требует найти:

(2¹ * 1!) / 1¹ + (2² * 2!) / 2² + (2³ * 3!) / 3³ + ...

А что ищется в последнем примере?

(2¹ * 1!) / 1 + (2² * 2!) / (1*2) + (2³ * 3!) / (1*2*3) + ...

[Lapp]

А что ищется в последнем примере?

(2¹ * 1!) / 1 + (2² * 2!) / (1*2) + (2³ * 3!) / (1*2*3) + ...

volvo, не догоняю.. Почему?..
Внутренний цикл по i, а делю я на n всякий раз.. У тебя получается внизу тоже факториал. Значит, факториал сокращается, остается только 2ⁿ. Эта величина всегда больше единицы! Я гонял свою прогу, она выдает тот же самый ответ, что и первый вариант (конечно, при малых n).

Я что-то упускаю?..

[volvo]

Lapp, извини, немного не так прочитал цикл, потом уже в машине понял что ошибся, но времени ответить не было... Все верно в твоей программе.

[kess]

а если так?? там возведение в степень +факториал.......

Program posledovatelnost;
Uses CRT;
var
  a,S,E:real;
  k,n:integer;

begin
ClrScr;
  WriteLn('VVedite to4nost E: ');
  ReadLn(E);
  WriteLn('Vvedite n ');
  ReadLn(n);

    for k:=1 to n-1 do
    a:=(exp(n*ln(2)))*(k*(k+1))/exp(n*Ln(n));
    If abs(a)>=E Then S:=S+a else S:=0;
    WriteLn('Symma ',n,'-x 4lenov rjada pavna= ', s:4:2);
    Repeat until Keypressed;
end.

[volvo]

А откуда ты знаешь, сколько элементов тебе надо будет просуммировать? Тебе по условию N не дано...

[kess]

А откуда ты знаешь, сколько элементов тебе надо будет просуммировать? Тебе по условию N не дано...

обычно любой можно задать

[Lapp]

обычно любой можно задать

Это как это??
Обычно много чего можно... Но речь идет не о том, что обычно можно, а о том, что сейчас нужно!

Но это не единственная твоя ошибка..

а если так?? там возведение в степень +факториал.......

Попытаюсь сделать первый проход по ошибкам (типа нулевое приближение ).

    { 1. Перед циклом весьма желательно обнулять переменную для суммы }
    for k:=1 to n-1 do  { 2. Тут нужен begin, видимо, иначе оператор If не в цикле.  И почему цикл до n-1? }
    a:=(exp(n*ln(2)))*(k*(k+1))/exp(n*Ln(n));  { 3. Расчет факториала неверный }
    If abs(a)>=E Then S:=S+a else S:=0;   { 4. Зачем обнуление суммы здесь?..  }
    { 5. А тут нужен end к begin'у (см. п.2 }

И вообще, нужен, конечно, не for, а repeat (или while).
Извини, просто любопытство: а чем тебе не понравился мой последний вариант? Что-то непонятно?