Модераторы: Poseidon
  

Поиск:

Ответ в темуСоздание новой темы Создание опроса
> [Pascal] Реализация сложной дроби 
V
    Опции темы
InviZible
Дата 23.9.2006, 15:23 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 63
Регистрация: 13.9.2006

Репутация: нет
Всего: нет
Есть задача, вроде всё понятно в ней, но
как можно реализовать вот эту дробь в паcкале:


(-1)^i *x^(4i+1)
-----------------
(4i+1)*(2i)!


"^" - "в степени"
"*" - умножить
"!" - факториал

"i" в ходе программы увеличивается (Inc(x);)
PM MAIL   Вверх
Kuvaldis
Дата 23.9.2006, 15:35 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


механик-вредитель
***


Профиль
Группа: Участник Клуба
Сообщений: 1189
Регистрация: 16.6.2006
Где: Минск

Репутация: 32
Всего: 61
InviZible,
Код

 if (i mod 2 = 0) then
   sign := 1
else
   sign := -1;
 sum := sign * exp(ln(x) * (4 * i + 1)) / ((4 * i + 1) * fact(2 * i) )

В Паскале нет функции возведения в степень, но есть экспонента (е в степени х)
 и есть логарифм 
 smile 

Факториал можно считать отдельной функцией, но если у тебя идет накапливание суммы вычислений, то лучше сделать так:
Код

  fact := 2;
  sign := -1;
  for i := 1 to n do
  begin
     temp := sign * exp( ln(x) * (4 * i + 1) ) / ( (4 * i + 1) * fact  );
     sign := sign * (-1);
     fact := fact * (2 * i  + 1) * (2 * i + 2);
     sum := sum + temp;
  end;

Т. е. переход к очередному факториалу производится путем простого домножения старого факториала на 2 недостающих множителя

Это сообщение отредактировал(а) Kuvaldis - 23.9.2006, 15:43


--------------------
Помни - когда ты спишь, враг не дремлет
Спи чаще и дольше, изматывай врага бессоницей
PM MAIL ICQ   Вверх
InviZible
Дата 23.9.2006, 21:10 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 63
Регистрация: 13.9.2006

Репутация: нет
Всего: нет
Большое спасибо, буду пробовать!
PM MAIL   Вверх
volvo877
Дата 24.9.2006, 01:09 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Эксперт
****


Профиль
Группа: Комодератор
Сообщений: 2073
Регистрация: 15.11.2004

Репутация: 34
Всего: 116
Цитата(Kuvaldis @  23.9.2006,  15:35 Найти цитируемый пост)
В Паскале нет функции возведения в степень, но есть экспонента (е в степени х) и есть логарифм 

И ни то, ни другое для решения данной задачи не нужно:

Код
...
sign := -1; fact := 2; { для первой итерации - 2!}
newX := X * sqr(sqr(X)); { для первой итерации - X^5 } 

for i := 1 to n do begin
  T := (sign * newX) / ((4 * i + 1) * fact);
  newX := X * sqr(sqr(X)); { увеличиваем степень на 4 }
  fact := fact * (2 * i + 1) * (2 * i + 2);
  sign := - sign;
  s := s + T;
end;
...
PM MAIL   Вверх
InviZible
Дата 30.9.2006, 09:30 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 63
Регистрация: 13.9.2006

Репутация: нет
Всего: нет
Большое всем спасибо. Всё реализовал. Попозже выложу
PM MAIL   Вверх
InviZible
Дата 7.10.2006, 12:37 (ссылка) | (нет голосов) Загрузка ... Загрузка ... Быстрая цитата Цитата


Шустрый
*


Профиль
Группа: Участник
Сообщений: 63
Регистрация: 13.9.2006

Репутация: нет
Всего: нет
Вот та часть этой дроби:
Код


 .....
   if (i mod 2) = 0 then // (-1)^i
    st1:=1
   else st1:=-1;

   vd:=st1*vd*sqr(x)*sqr(x); // числитель

   fact:=1;            // факториал
    for j:=1 to 2*i do
     fact:=fact*j;

   nd:=(4*i+1)*fact; // знаменатель
      
   slag:=(vd/nd); // вычисление дроби
 ....  
 
   

Всем спасибо
PM MAIL   Вверх
  
 Ответ в темуСоздание новой темы Создание опроса
Правила форума "Центр помощи"

ВНИМАНИЕ! Прежде чем создавать темы, или писать сообщения в данный раздел, ознакомьтесь, пожалуйста, с Правилами форума и конкретно этого раздела.
Несоблюдение правил может повлечь за собой самые строгие меры от закрытия/удаления темы до бана пользователя!


  • Название темы должно отражать её суть! (Не следует добавлять туда слова "помогите", "срочно" и т.п.)
  • При создании темы, первым делом в квадратных скобках укажите область, из которой исходит вопрос (язык, дисциплина, диплом). Пример: [C++].
  • В названии темы не нужно указывать происхождение задачи (например "школьная задача", "задача из учебника" и т.п.), не нужно указывать ее сложность ("простая задача", "легкий вопрос" и т.п.). Все это можно писать в тексте самой задачи.
  • Если Вы ошиблись при вводе названия темы, отправьте письмо любому из модераторов раздела (через личные сообщения или report).
  • Для подсветки кода пользуйтесь тегами [code][/code] (выделяйте код и нажимаете на кнопку "Код"). Не забывайте выбирать при этом соответствующий язык.
  • Помните: один топик - один вопрос!
  • В данном разделе запрещено поднимать темы, т.е. при отсутствии ответов на Ваш вопрос добавлять новые ответы к теме, тем самым поднимая тему на верх списка.
  • Если вы хотите, чтобы вашу проблему решили при помощи определенного алгоритма, то не забудьте описать его!
  • Если вопрос решён, то воспользуйтесь ссылкой "Пометить как решённый", которая находится под кнопками создания темы или специальным флажком при ответе.

Более подробно с правилами данного раздела Вы можете ознакомится в этой теме.

Если Вам помогли и атмосфера форума Вам понравилась, то заходите к нам чаще! С уважением, Poseidon, Rodman
 
1 Пользователей читают эту тему (1 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема | Центр помощи | Следующая тема »

Подписаться на тему | Подписка на этот форум | Скачать/Распечатать тему


 




[ Время генерации скрипта: 0.1160 ]   [ Использовано запросов: 21 ]   [ GZIP включён ]


Реклама на сайте     Информационное спонсорство

 
 По вопросам размещения рекламы пишите на vladimir(sobaka)vingrad.ru
Отказ от ответственности     Powered by Invision Power Board(R) 1.3 © 2003  IPS, Inc.