VPF::решение кубического уравнения

DeMoNN

Дата 17.2.2009, 01:19 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

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

Репутация: нет
Всего: 2

в инете нашел несколько алгоритмов на си, яваскрипте и матлабе, но в паскале не могу найти некоторых функций-з, извлечение куб корня(а лучше-возведение в рациональную степень), например, или acos() и pow()есть ли их аналоги на паскале?

SirGrey

Дата 17.2.2009, 09:51 (ссылка)

(голосов:1)

Загрузка ...

Новичок

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

Репутация: нет
Всего: нет

Решил быть программистом, Учи математику!!! Логарифмы рулят)))
Найти кубический корень - то же, что возвести в степень 1/3.

Следовательно Q:=Exp((Y*Ln(X)), а для кубического корня соответственно Q:=Exp(Ln(x)/3)

DeMoNN

Дата 17.2.2009, 10:20 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

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

Репутация: нет
Всего: 2

спасибо!
З.Ы. С математикой у меня вроде плохо небыло(на уровне моего 8го класса)

volvo877

Дата 17.2.2009, 10:47 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

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

Репутация: 2
Всего: 116

Цитата(SirGrey @ 17.2.2009, 08:51

)

Решил быть программистом, Учи математику!!!

Решил давать советы на форуме - учи математику: вычисли своим методом кубический корень от (-3)... Что получилось, не расскажешь?

DeMoNN

Дата 17.2.2009, 11:32 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

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

Репутация: нет
Всего: 2

А можно так?

Код


function cbrt(base: extended): extended;
begin
  if base >= 0 then
    result:=power(base, 1/3)
  else
    result:=-power(-base, 1/3);
end;

а потом юзать эту функцию, где нужно

Это сообщение отредактировал(а) DeMoNN - 17.2.2009, 11:43

DeMoNN

Дата 17.2.2009, 12:06 (ссылка)

(нет голосов)

Загрузка ...

Бывалый

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

Репутация: нет
Всего: 2

вот, вроде работает:

Код


procedure cburan(a, b, c, d: extended);
var
  r, s, t, p, q, v, u, discrim: real;
  x1, x2, x3: string;
begin
  r := b / a;
  s := c / a;
  t := d / a;
  p := (3*s - r*r) / 3;
  q := (2*r*r*r)/27 - (r*s)/3 + t;
  discrim := (p*p*p)/27 + (q*q)/4;
  v := cbrt(-(q/2) - sqrt(abs(discrim)));
  u := cbrt(-(q/2) + sqrt(abs(discrim)));
  x1 := floattostr(u + v - r/3);
  x2 := floattostr(-(u + v)/2-r/3) + '+' + floattostr(((u-v)/2)*sqrt(3)) + 'i';
  x3 := floattostr(-(u+v)/2-r/3) + '-' + floattostr(((u-v)/2)*sqrt(3)) + 'i';
  showmessage('X1=' + x1 + #13 + 'X2=' + x2 + #13 + 'X3=' + x3);
end;

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| Object Pascal: кроссплатформенные технологии \| Следующая тема »