Задачка, Если не угадал с разделом, подкорректируйте, плз Опции

[Олири]

Девушка попросила помочь, сказала, что решать (как я понял) надо с использованием Паскаля.
Как объяснила: "надо решить сложную задачу на циклические вычислительные процессы..."

Текст задачи: Для заданных значений m и x вычислить бином Ньютона (1+x)^m непосредственно и по формуле разложения в ряд 1+x)^mСⁱ_mXⁱ_m= Для вычисления Сⁱ_m использования рекуррентное соотношение Сⁱ⁺¹_m=Сⁱ_m(m-i)/(i+1),C⁰_m=1

Было бы неплохо, если бы кто помог) я в этом ничего не соображаю)

Сообщение отредактировал Олири - Dec 27 2010, 20:18

[iMagus]

С разделом может и не ошибся.

ЗЫ. Методичка есть?)

Это задача по программированию? На каком языке?

А, узрел - Паскаль

ЗЗЫ. На кого девушка учится? Возможно, ей стоит сменить будущую профессию =\

Сообщение отредактировал iMagus - Dec 27 2010, 20:29

[Олири]

С разделом может и не ошибся.

ЗЫ. Методичка есть?)

Это задача по программированию? На каком языке?

А, узрел - Паскаль

ЗЗЫ. На кого девушка учится? Возможно, ей стоит сменить будущую профессию =\

так я эта девушка и учусь я в на инженера -строителя..
это операторы цикла с параметром

[iMagus]

Зачем ты учишься на инженера-строителя?)

Код

function calc_binom(x: integer, m: integer): integer;
begin
calc_binom := pow(1+x,m);
end;

function calc_binom_series(x: integer, m: integer): integer;
var
  sum: real;
  ci, ci1: real;
  i: integer;
begin
i := 1;
ci := 1;
sum := 1;
while i <= m do
  begin
    ci1 := ci*(m-i)/(i+1);
    sum := sum + ci1*pow(x,i);
    ci := ci1;
    i := i + 1;
  end;
calc_binom_series := sum;
end;

Проверить негде сейчас.

[Олири]

ого..спасибо))))сейчас проверю)) а учусь..ну мне нравиться)))всё....кроме этой задачи))...

Сообщение отредактировал Олири - Dec 28 2010, 14:59

[Олири]

неправильно(( при подстановке x=2 m=3 ответ должен быть 27 а не 5((

[Шантина]

http://www.cyberforum.ru/

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

Иногда там тоже от нечего делать помогаю, правда по C#, и заместо высшей математики мне в инсте в голову вложили только опилки))

[iMagus]

неправильно(( при подстановке x=2 m=3 ответ должен быть 27 а не 5((

Значит ты формулу для вычисления биномиальных коэффициентов неправильно написал
Или я чего-то пропустил.

Я вот про эту:

Для вычисления Сⁱ_m использования рекуррентное соотношение Сⁱ⁺¹_m=Сⁱ_m(m-i)/(i+1),C⁰_m=1

[Олири]

балин в условии и впраду рамс ...вот правильное http://www.cyberforum.ru/pascal/thread209046.html

[iMagus]

Что такое рамс? Регистрироваться на том форуме надо. Чуть позже гляну правильное условие.

[Олири]

да я не знаю как тут напечатать кое что и поэтому кидаю ссылку..

[iMagus]

Да нет, все правильно ты написал в 1м посте. Надо подумать, в чем подвох

ЗЫ. Фото экрана монитора - это круто

[Олири]

Да нет, все правильно ты написал в 1м посте. Надо подумать, в чем подвох

ЗЫ. Фото экрана монитора - это круто

неправильно..там сумма от 0 до m...

[iMagus]

неправильно..там сумма от 0 до m...

Это и так было известно.

Я сейчас эту же рекуррентную "формулу" вывел, но она не работает. Тут наверно какая-то методическая ошибка имеет место быть. Настоящий математик объяснил бы наверно ее.

Сообщение отредактировал iMagus - Dec 28 2010, 22:13

[iMagus]

В общем, из-за хитрых свойств факториала предложенная формула работать не будет. Либо смысл задания в том и есть, чтобы показать, что формула неправильная. Либо следует уточнить задание.

ЗЫ. К слову сказать, не очень понятно, почему у тебя получилось в ответе - 5

неправильно(( при подстановке x=2 m=3 ответ должен быть 27 а не 5((

Должно было получиться 4.3333(3)

Скорее всего, потому, что надо было вот эту строку

Код

    ci1 := ci*(m-i)/(i+1);

писать как

Код

    ci1 := ci*(m-i)/(i+1.0);

Паскаль такой Паскаль

[Олири]

нет там никакой ошибки i и i+1 просто обозначают следующий элемент а не сами степени..

[iMagus]

нет там никакой ошибки i и i+1 просто обозначают следующий элемент а не сами степени..

Ясен пень, что не степени, а след элемент. На то она и рекуррентная формула эта.

Но я нашел багу (зря я формулу ругал).
Так будет правильно

Код

function calc_binom_series(x: integer, m: integer): integer;
var
  sum: real;
  ci, ci1: real;
  i: integer;
begin
i := 1;
ci := 1;
sum := 1;
while i <= m do
  begin
    ci1 := ci*(m-i+1.0)/(i+1.0-1.0);
    sum := sum + ci1*pow(x,i);
    ci := ci1;
    i := i + 1;
  end;
calc_binom_series := sum;
end;

ЗЫ. Более наглядна формула была бы в виде:
Сⁱ_m=С^i-1_m(m-(i-1))/((i-1)+1)

Сообщение отредактировал iMagus - Dec 29 2010, 11:44

[Firehitman]

Теперь понятно, почему у нас все так криво строят... кирпичи по формулам кладут, да еще и на форумах списанных...

[iMagus]

Да, и никто, главное, не заметил, где я ошибся поначалу. Пришлось даже неработу алгоритма списать на мистические свойства факториала. А ведь ошибка-то была на поверхности - только я ее не видел, потому, что считал, что в данном месте ошибиться невозможно. Стыдоба!

ЗЫ. Фаер, ты о чем?)

[Олири]

ну он имеет ввиду что я тупая...у вас можт и строят криво)) а я не из владивостока