Vladimir_S

Вы забыли добавить "милсдарь", а также уведомить меня, что Вы к моим услугам.
А эта сентенция к кому обращена? Ко мне? К разработчикам ЯП? Ежели ко мне, то мне почему-то кажется, что создание библиотек констант, функций, модулей и т.п. определяется моими в оных потребностями, а не абстракциями типа "давно пора".
P.S. Кстати, библиотеки у меня изрядные: спецфункции и т.п. Поскольку нужны реально. А вот особой потребности, скажем, в Ln(π) как-то вот не испытываю. Извините.

Uncle Vasya

Господа, да вы шутите? С очень высокой точностью вы оба получили нуль! Я уверен, что вы знаете как сравниваются вещественные числа, но если вдруг нет... Никогда не сравнивайте вещественные числа напрямую, только если не надо сравнивать на больше или меньше. Равенство вещественных чисел проверяют по вычитанием одного числа из другого, взятия модуля и сравнения с малой константой. Как-то так:

if Abs(x1 - x2) > 0.0000001 then

Иначе можете нарваться на ситуацию, когда

0.1 + 0.2 != 0.3

Простите за оффтоп.
В Паскале с его синтетическим типом Real вообще, на мой скромный взгляд, беда. Он ни аппаратно не поддерживается, да ещё и требует затрат времени на конвертации. А при подключении поддержки сопроцессора так вообще превращается в Extended! На эти грабли я наступил впервые ещё в университете в 90-ые. Чего потом удивляться что значение косинуса вычисляется по-разному? Не зря чуваки из Sun Microsystems очень долго пилили свою Java чтобы обеспечить одинаковые результаты вычислений на всех поддерживаемых платформах.

Viewer

Вы хотите "на шпагах"?
Смешно, с Вашим Паскалем.
Если докажете, что он Ваш.

Viewer

Это Вы извините, что напряг Вас с логарифмом Пи - это особая функция, что помогает входу члена в удлиненное влагалище.

Vladimir_S

Не согласен. Е-19 - да, а вот Е-12 как-то IMHO многовато для нуля.

black88

Брадису расскажите с его 8-значными таблицами

Uncle Vasya

Ну, если рассчитывать траекторию полёта космического аппарата куда-нибудь в пояс Койпера, то вполне возможно, потребуется точность выше 12 знаков, но для большинства применений и 6, а то и 3 хватит. Правда, автор темы не озвучил никаких особых требований.

Точный вещественный нуль можно получить только если задать константу 0 либо вычесть друг из друга два одинаковых вещественных числа (которые, в свою очередь, точно можно задать лишь константами). Результаты любых математических вычислений с вещественными числами может быть лишь приблизительным и очень сильно зависит от самих чисел из-за внутреннего представления.

Vladimir_S

И проронил я скупую мужскую слезу...

Dimafania

Друзья, поскольку я вижу, что здесь собрались действительно профессионалы своего дела, попрошу все таки вас предложить свое решение на мою задачу. В принципе программу я и сам уже написал в некоторых версиях, но вот видится мне, что вы можете предложить куда более оптимальное решение (ввиду большого кол-ва циклов программа считает очень долго).
Итак, у нас есть эксп. значения теплового расширения при ряде температур:
TEC:array [1..20] of Real = (-0.4E-6, -1.9E-6, -4E-6, -4.5E-6, -4.38E-6, -3.9E-6, -2.9E-6, -2.23E-6, -1.8E-6, -1.7E-6, -1.7E-6,
-1.9E-6, -2.1E-6,-2.4E-6, -2.78E-6, -2.9E-6, -2.95E-6, -3.12E-6, -3.25E-6, -3.22E-6);
T:array [1..20] of Integer = (20,28,40,50,60,70,80,90,100,110,120,130,140,150,1 60,170,180,190,200,210);

Задача подогонать модель под экспериментальные данные (т.е. по наименьшему среднеквадратическому отклонению).
Модельное тепловое расширение рассчитывается по формуле (1):

ΔV,V0,R - константы, что у меня есть;
с - фиксированная при определенном Т, варьируется от 0 до 1 (желаемая точность определения 0.01)
ΔS, U - константы, которые следует определить, неизменны при всех Т.
ΔS, U, c(T) и новая константа ΔЕ, которую также надо найти, должны удовлетворять следующим условиям:
Формула (2):

Итак первая производная формулы 2 должна быть равна нулю, вторая производная формулы (2) должна быть положительной. Р везде берется равное 1.
В результате программа должна находить ΔЕ, ΔS, U, фиксированные при всех Т, а также ряд значений с(Т), такими, чтобы среднеквадратическое отклонение параметра B (формула 1) от эксперименитального значения TEC было минимальным.
В случае если аналитического решения нет, то перебор должен быть в следующих диапазонах:
ΔS от -7 до 7 с шагом 0.01
ΔE от от -700 до 700 с шагом 1
U от 0 до 5000 с шагом 1
с от 0 до 1 с шагом 0.01

Dimafania

Забыл добавить, что первая и вторая производные частные по с.