Главная Промышленная автоматика.

for m:=2 step 1 until n do

begin s:= (kXr-dr+ (m-1) Xs) Xda; k:=-к end m; bs:=s; fin: end bs;

Свидетельство к алгоритрлу 1096

Алгоритм 1096 получен в результатесущественного усовершенствования алгоритма 109а в следующих отношениях.

1. Исключение из алгоритма 109а локального массива Ь[1:п] и четырех локальных переменных типа real (с ц&лью экономии машинной памяти).

2. Объединение двух циклов, начинающихся метками even и odd, в один цикл, начинающийся меткой iter (с целью сокращения записи процедуры bs).

3. Включение глобальных идентификаторов alim и deltal09 в список формальных параметров (с целью повышения удобств пользования процедурой в библиотеках трансляторов, подобных архиву БЭСМ-АЛГОЛ).

После внесения вышеуказанных изменений алгоритм был транслирован в системе 4-70, и с его помощью были повторены расчеты для значений п-1, 2, ..., 16, а-0.25 и 24.0. Полученные результаты во всех выданных на печать восьми цифрах совпали со значениями Bn-iia), приведенными в табл. 3 (по данным Миллера и др.).

Свидетельство к алгоритму 109а

Алгоритм 109а получен в результате модификации, сокращения и ординарной переработки алгоритма 109 (Кругляк Ю. А., Whitman D. R. «САСМ», 1962, № 7).

Модификация алгоритма 109 заключалась в следующем.

1. В алгоритме 109 использовались локализованные переменные alim и gamma, которые нужно было заменять в теле процедуры числами в каладохм конкретном случае применения алгоритма 109. Б алгоритме 109а вместо трех переменных alim, delta и gamma использовались две глобальные переменные alim и deltalOQ, в результате чего тело процедуры bs стало независимым от конкретного применения.

2. Отношение mod{n-1,2):=0 с процедурой-кодом mod было заменено более простым и универсальным отношением й-2Х2~п.

Алгоритм транслирован с исходными данными из «Подтверждения к алгоритму 109», приведенного а.вторами, для п=1, 2, 9. Результаты совпали с прелними с точностью до семи значащих цифр (см. табл. 3).

Подтверждение к алгоритму 109 -И

Интегоалы 5„(а)= Г л:" ехр (-ax)dx встречаются в задачах физк-

ки, использующих сфероидальные бицентрические координаты, в частности в квантовохимических расчетах.



Алгоритм был запрограммирован для вычислительной машины Burroughs 220, использующей транслятор Burroughs Algebraic Compiler. Программа была применена для вычисления таблиц Вп{а) при п=0, 1, 15 и 0=.00, 0.01, 32.54. Например, для л=0, 1, 15; а=0.25 и а=24.0 были получены результаты, приведенные в табл. 3. Они сравнивались с результатами (столбцы 3 и 5) из работы Миллера и др. [9i].

Таблица 3

а=..0.25

а=24.0

Алгоритм 109

По данным Миллера и др.

Алгоритм 109

По даннЕШ Миллера и др

2 3 4 5 6 7 8 9 10 II 12 13 14 15 16

О -О О

-О, О, -О,

о, -о,

0. -0.

0. -0.

0. -0.

0. -0.

.20208984X10

.16771064X10»

.67921322X10»

,10074584X10»

.40896479X10»

,72008754X10-»

,29268836X10»

,56030292X10-"

22792911X10»

45856272X10-»

,18664760X10»

388097 8X10-»

15803198X10»

33640Б62Х10-»

13702696X10»

29686662X10-»

0.20208985344653X10» -0.16771066П7520Х10»

0.67921324506375X10» -0.10074585827159X10»

0.40 96480211998X10» -0.7200875663692 ХЮ"»

0.292688375 7905x10» -O.5603O294023I7 ХЮ"

0.22792912573392X10" -0. 5856272975462X10-»

0. 8664761544688X10» -0.38809719373731x10-»

0.15803200452627X10" -О.33640 63670387ХШ-»

0.13702696892367X10» -0.29686663616401X10-»

0.11037134x10" -0.10577253X10»

0.10155696X10" -0.97676725X10»

0.94091887X10» -0.90768866X10"

О 76791 9X10" -0.847982 2X10"

О 8210 25 ХЮ" -0.79581 7 ХЮ"

0.7721 22 ХЮ" -0.7498 404x10"

0.728 141X10" -0.70894600X10»

0.69016158X10» -0.67236245ХЮ"

0.110371342208ХЮ» -0.105772536282X10"

0.101ББ6964184ХЮ» -0.976767216847x10"

0.940918885936x10» -0.907688634174x10"

0.876921258533X10" -0.847982638338ХЮ"

0.821052542631x10» -0.795818718590X10»

0.772122289331X10" -0.74S824039467X10"

0.728801402343X10» 0.7П"94--907х10"

О 901 1 91 89 Ю» -0.0/2362427583X10"

В упомянутых интервалах точность не менее шести значащих цифр. Такая точность достаточна для большинства квантовохимических расчетов.

Дополнительные сведения об алгоритмах 108 и 109

В личной переписке автор алгоритмов 108 и 109 Ю. А. Кругляк сообщил редактору настоящего выпуска и просил опубликовать следующие дополнительные сведения.

1. Подробные блок-схемы, технические характеристики алгоритмов 108 и 109 и сами алгоритмы на языке АЛГОЛ-58 онубликоваиы в «Журнале структурной химии», 1962, № 5, с. 569.

2. Оба алгоритма были попользованы для вычисления подробных таблиц интегралов А и В, изданных в работе: Кругляк Ю. А.,

Уитмен д. р. Таблицы интегралов Л„(1, а)= x"e-dx и В„{а) -

г yzg axj-g Хаблицы интегралов квантовой химии. Т. L

М., ВЦ АН СССР, 1965.

Свидетельство к алгоритму 1106 [S22]

Алгоритм 1106 «Квзнтовомеханичеокие интегралы для орбит электронов» не публикуется здесь потому, что соответствующий алгоритм ПО (Кругля.к Ю. А., Whitman D. R. «САСМ», 1962, № 7) содержит ряд ошибок, из которых можно отметить следующие, бросающиеся в глаза при начале раэбора алгоритма.

1. Описание процедуры INTSOLI начинается с описателя real, хотя в остальном эта процедура не является функцией.



2. Описания массивов а, b, G и be находятся среди спецификаций, а не в начале тела процедуры,

3. Описания процедур as и bs находятся среди операторов тела процедуры INTSOLI, а не в начале тела последней.

4. Делается попытка вычислить оператор

fo:=(/2) f (п+З) X (2XZSG) f («+l/2) Xzsb (5/2) X G[2Xn] f (1 2)

с неопределенным значением переменной t.

5. Алгоритм 110 содержит много опечаток. Например, в первых же операторах вместо переменной q напечатана переменная g, вместо переменной Ьс[1,1] -переменная Ьс[1, 1) и т. д.

В личной переписке автор алгоритма ПО Ю, А. Кругляк сообщил редактору данного выпуска и просил опубликовать следующие сведения.

1. На языке АЛГОЛ-58 вместе с подробным изложением математики алгоритм опубликован в «Журнале структурной химии», 1963, № 1, с. 91.

2- Алгоритм использоващ для составления таблиц интегралов, изданных в работе: Кругляк Ю. А. Таблицы некоторых молекулярных двухцентровых одноэлектронных интегралов.- «Таблицы интегралов квантовой химии», т. II. Харьковск. ГУ, 1966.

Свидетельство к алгоритму f 116 [S22]

Алгоритм 1116 «Молекулярно-орбитальные расчеты молекулярных взаимодействий» (Кругляк Ю. А., Whitman D. R. «САСМ», 1962, № 7) не публикуется здесь по тем же причинам, что и алгоритм ПО.

В личной переписке автор алгоритма 111 Ю. А. Кругляк сообщил редактору данного выпуска и просил опубликовать следующие сведения об алгоритме: подробное математическое описание и сам алгоритм 111 на языке АЛГОЛ-58 опубликован в «Журнале структурной химии», 1963, № 2, с. 254. Об использовании алгоритма даются сведения в «Трудах химического факультета и научно-исследовательского института химии Харьковск. ГУ», 1963, № 19, с. 179.

АЛГОРИТМ 1126

Положение точки относительно многоугольника [Z]

Если точки {x[i], y[i\) (t=l, 2, п) являются верщинами простого замкнутого многоугольника, пронумерованными в циклическом порядке, а {хО, г/0) -точка, не лежащая ни на одной из сторон многоугольника, то процедура-функция pointpol (сокращение от point - точка, ро/г/gon-многоугольник) определяет (принимая при этом значение true), лежит ли точка (хО, г/0) внутри многоугольника. При обращении к процедуре pointpol массивы х я у должны иметь размерность [1:п+1].

Boolean procedure pointpol (n,x,y,xO,yO);

value n,xO,yO; real xO,yO; integer n; array x,y; begin integer i; Boolean b;





0 1 2 3 [4] 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41

0.0017