Головна сторінка
Top.Mail.Ru Яндекс.Метрика
Форум: "Основна";
Поточний архів: 2002.01.08;
Завантажити: [xml.tar.bz2];

Вниз

Рівняння третього ступеня. Знайти схожі гілки


mashinist   (2001-12-16 23:16) [0]

Я тут, схоже всіх вже задовбав зі своєю математикою, але все ж. У кого-небудь є шматок коду для вирішення рівняння третього ступеня виду ax ^ 3 + bx ^ 2 + cx + d = 0, причому a, b, c, d: extended;

Допоможіть!



Builder   (2001-12-17 02:35) [1]

А чи не легше ніж довбати, сісти і написати самому?

Берещь будь-яку книжку з чисельних методів - це є в першому розділі.



mashinist   (2001-12-17 11:58) [2]

Не виходить.
Вже пробував.



Romkin   (2001-12-17 12:19) [3]

http://alglib.chat.ru/equat/index.html#cube
По-моєму, вичерпно
Але, взагалі-то, рівняння вище другого ступеня зазвичай вирішують не з точних формулах, а по алгоритмам знаходження коренів поліномів n-го ступеня



savva   (2001-12-17 12:31) [4]

ось що запропонував редактор блок схем
Procedure CubeEquation (a, b, c: real; var nr: byte; x: array [1..3] of real); // струму чому то не описані параметри 8)) сам опишеш begin p: = - a * a / 3 + b; q: = 2 * a / 3 * a / 3 * a / 3-a * b / 3 + c; QH: = (p / 3) * (p / 3) * (p / 3) + (q / 2) * (q / 2); if QH = 0 then begin nr: = 3; AH: = abs (q / 2); if AH <> 0 then begin AH: = sign (q) * exp (ln (AH) / 3) end; x [1]: = 2 * AH-a / 3; x [2]: = - AH-a / 3; x [3]: = - AH-a / 3 end else begin if QH> 0 then begin AH: = - q / 2 + sqrt (QH); AH: = sign (AH) * exp (ln (abs (AH)) / 3); BH: = - q / 2-sqrt (QH); BH: = sign (BH) * exp (ln (abs (BH)) / 3); x [1]: = AH + BH-a / 3; if AH = BH then begin nr: = 3; x [2]: = - (AH + BH) / 2-a / 3; x [3]: = - (AH + BH) / 2-a / 3 end else begin nr: = 1; x [2]: = - (AH + BH) / 2-a / 3; x [3]: = (AH-BH) / 2 * sqrt (3) end; end else begin nr: = 3; u: = - q / 2 / sqrt (-p * p * p / 27); u: = arccos (u); x [1]: = 2 * sqrt (abs (p / 3)) * cos (u / 3) -a / 3; x [2]: = - 2 * sqrt (abs (p / 3)) * cos (u / 3 + Pi / 3) -a / 3; x [3]: = - 2 * sqrt (abs (p / 3)) * cos (u / 3-Pi / 3) -a / 3 end; end; end;

не знаю, чи варто довіряти ...



SergVlad   (2001-12-17 13:58) [5]

To mashinist
Все ще борешся зі своїми регуляторами?



Юрий Зотов   (2001-12-17 14:59) [6]

> Romkin © (17.12.01 12: 19)

Дозволю собі поправку - вище третьої, а не другий. Для кубічних же рівнянь є точні формули - наприклад, формула Кардано. Це набагато простіше, швидше і точніше.



mashinist   (2001-12-17 15:40) [7]

Спасибі!




сторінки: 1 вся гілка

Форум: "Основна";
Поточний архів: 2002.01.08;
Завантажити: [xml.tar.bz2];

наверх









Пам'ять: 0.57 MB
Час: 0.025 c
14-22354
Немезида
2001-10-08 13:02
2002.01.08
мобіла


1-22096
Mikhalyov Dmitry
2001-12-19 15:55
2002.01.08
Про TObject в подіях


3-21942
оборот
2001-12-03 12:19
2002.01.08
Добрі люди допоможіть з збереженої процедурою.


1-22202
T2
2001-12-16 07:12
2002.01.08
До питання про маленьких виконуваних файлах


1-22047
handra
2001-12-17 13:12
2002.01.08
створення довідки





африкаанс албанський арабська вірменин азербайджанець баскський білоруський болгарська каталонський Китайська (спрощене письмо) Китайський традиційний) хорватський чеська данську мову нідерландський Ukranian естонець Філіппінська фінську мову французький
галісійська грузинський німецький грецький гаїтянський креольський давньоєврейську хінді угорський ісландський індонезієць ірландський італійський японський корейський латиська литовець македонець малайський мальтійський норвежець
перс полірування португальська румунський російська сербський словацький словенський іспанська суахілі шведську мову тайський турецька український урду в'єтнамський валлійський ідиш бенгальський боснійський
кебуано есперанто гуджараті хауса хмонг ігбо яванський каннада кхмерская Лао латинь маорі маратхі монгольський непальська панджабі сомалійський тамільська телугу йоруба
зулуський
Англійська Французький Німецький Італійський Португальська Русский Іспанська