Библиотека математических функций, (разделено) Опции

[ProtasSoft]

Коллеги!
Мне нужна библиотека математических функций для Pascalя.

Подскажите где можно что-то подобное посмотреть.

Сообщение отредактировано: ProtasSoft - 18.08.2005 14:24

[volvo]

Что именно тебе нужно? Математические функции - это слишком расплывчатое понятие, конкретнее: может работа с матрицами, может - тригонометрия, может - комплексные числа?

+ к этому придерживаемся правила - "новый вопрос - новая тема"... ;)

Сообщение разделено из темы "Обработка прямоугольных матриц"

[ProtasSoft]

Мне нужна программа, реализующая хи-квадрат-распределение с M степенями свободы.

[volvo]

Теорию можешь посмотреть здесь (есть примеры программ, только на С):
http://www.geo.pu.ru/ecobez/edu/books/eis/mstat/test3.htm

и здесь:
http://www.is.isea.ru/~i_gestrin/common/in...has&i=118&r=118

[ProtasSoft]

теория мне не нужна
Формула функции, которую нужно реализовать на Паскале, следующая:

F(M,x)={1/2^(M/2)*Г(М/2)}*Integral (от 0 до x) {t^{M/2-1}exp(-t/2)dt}, где

Г(y)=Integral (от 0 до бесконечности) {exp(-t)*t^(y-1)*dt} -- гамма-функция Эйлера.

Кто-нибудь может что-то предложить

[volvo]

Давай все-таки, наоборот: ТЫ предлагаешь свой вариант, и говоришь, где у тебя затруднения, ибо решать за тебя никому не нужно... Готового решения не жди...

P.S. Я бы на твоем месте вначале прочитал FAQ по процедурным типам, там есть кое-что полезное для тебя...

Сообщение отредактировано: volvo - 18.08.2005 15:25

[volvo]

Кстати, Хи-квадрат распределение затабулировано для 1 <= M <= 30. У меня есть PAS-файл содержащий таблицу значений этого распределения при вероятностях от 0.95 до 0.05, возможно этого тебе хватит. Если нужно - скажи...

[ProtasSoft]

у меня M может быть очень большим (>>30).

[ProtasSoft]

Кстати нашел интересную мат библиотечку. Если кому интересно см. аттач.

[ProtasSoft]

Коллеги. подскажите как вышепрведенную библиотеку подсоединить к компилятору FPC?

Сообщение отредактировано: ProtasSoft - 18.08.2005 16:08

[volvo]

как вышепрведенную библиотеку подсоединить к компилятору FPC?

Распакуй содержимое библиотеки в рабочую директорию и подключай нужные модули директивой Uses (так же, как и в Turbo/Borland Pascal)

Uses m_llong;
...

Или ты имеешь в виду, что есть проблемы с компиляцией кода в FPC?

[volvo]

Кстати, в связи с этим:

Г(y)=Integral (от 0 до бесконечности) {exp(-t)*t^(y-1)*dt} -- гамма-функция Эйлера.
Кто-нибудь может что-то предложить

возник вопрос - тебе обязательно находить гамма-функцию методом Эйлера, или можно воспользоваться тем, что

Г(n+1) = n!

, и соответсвенно

Г(n) = (n - 1)!

? Возможно, через факториалы будет проще...

[ProtasSoft]

Подскажите как стереть созданный ранее файл.
Такая запись верна

file1:='matrixdat.txt';
    erase(file1);
    file2:='dat.txt';
    erase(file2);

[klem4]

uses crt;
var
   f:text;

Begin

   Assign(f, 'c:\test.txt');
   erase(f);

end.

ps FAQ Файлы