автореферат диссертации по машиностроению и машиноведению, 05.02.13, диссертация на тему:Исследование влияния посадки патронов на шпиндель прядильных веретен с целью улучшения условий их эксплуатации

кандидата технических наук
Ушфун, Тауфик
город
Москва
год
1985
специальность ВАК РФ
05.02.13
цена
450 рублей
Диссертация по машиностроению и машиноведению на тему «Исследование влияния посадки патронов на шпиндель прядильных веретен с целью улучшения условий их эксплуатации»

Оглавление автор диссертации — кандидата технических наук Ушфун, Тауфик

ВВЕДЕНИЕ.

Глава I. ЛИТЕРАТУРНЫЙ ОБЗОР.

1.1. Краткое описание конструкции современных веретен

1.2. Краткий обзор литературы

1.3. Постановка задачи

Глава П. ТЕОРЕТИЧЕСКИЕ ИССЩОВАНИЯ ВЛИЯНИЯ ПРУЖИННЫХ

ФИКСАТОРОВ НА ДИНАМИКУ ВЕРЕТЕНА ТИПА ВИГ 28

С ПАКОВКОЙ

2.1. Динамическая модель

2.2. Динамика шпинделя веретена с фиксаторами вверху и внизу

2.3. Дифференциальные уравнения движения и их решение

2.4. Динамика шпинделя веретена, имеющего пружинные фиксаторы вверху и шарнир внизу

2.4.1. Амплитуды колебаний верхних концов насадка и паковки

2.4.1.1. Амплитуды колебаний верхнего конца насадка

2.4.1.2. Амплитуды колебаний верхнего конца паковки

2.4.2. Реакции в опорах шпинделя веретена

ВИГ 28

2.4.2.1. Реакции в опорах шпинделя веретена в плоскости

2.4.2.2. Реакции в опорах шпинделя веретена в плоскости

2.5. Определение коэффициентов влияния

2.6. Теоретические результаты и их анализ .VI 2.6.1. Анализ влияния различных параметров шпинделя на величину амплитуды колебания его верхнего конца с паковкой

2.6*2. Анализ влияния различных параметров шпинделя на величину динамических нагрузок в верхней опоре

2.7. Критические скорости вращения шпинделя веретена с паковкой

Глава Ш. ЭКСПЕРИМЕНТАЛЬНЫЕ ИССЛЕДОВАНИЯ.

3.1. Экспериментальное определение физических параметров веретена и паковки

3.1.1. Определение моментов инерции шпинделя (насадка) и паковки

3.1.2. Определение упругости опоры веретена

3.1.3. Определение величин динамической неуравновешенности шпинделя веретена и паковки.

3.2. Описание конструкции экспериментальной установки

3.2.1. Методика определения критических скоростей вращения шпинделя веретена с паковкой.

3.2.2. Методика определения амплитуд колебаний верхнего конца шпинделя веретена (паковки) .Ю

3.2.3. Методика определения динамических усилий в верхней опоре веретена

3.3. Проведение многофакторного эксперимента при исследовании влияния пружинных факторов на условия вращения веретена

3.3.1. Выбор плана экспериментальных исследований

3.3.2. Методика обработки результатов эксперимента

3.3.3. Анализ влияния факторов на условия вращения веретена

3.3.3.1. Анализ влияния факторов на величину динамических нагрузок в верхней опоре шпинделя веретена

3.3.3.2. Анализ влияния факторов на величину амплитуд колебаний верхнего конца шпивделя с паковкой

Введение 1985 год, диссертация по машиностроению и машиноведению, Ушфун, Тауфик

После освобождения Алжира от колониального господства в 1962 г., Партия (Фронт национального освобождения) и народ должны были окончательно освободить страцу от неоколониализма -новой формы эксплуатации фирмами и международными корпорациями. Это была нелегкая задача, если учесть, что в молодой республике 95% всего населения было неграмотным.

Национализация богатства страны стала возможной при поддержке СССР и других социалистических стран, которые с первых дней революции были на стороны алжирского народа и продолжали оказывать Алжиру помощь после победы над французскими колонизаторами.

С 1962 года в СССР были подготовлены тысячи алжирских специалистов. Кроме того в Алжире при содействии СССР были построены десятки заводов, институтов, техникумов и центров для подготовки кадров. Среди них можно назвать крупнейший в Африке металлургический комплекс Эль-Хаджар в городе Аннаба, нефтяной институт, институт легкой промышленности в Бумердесе и т.д.

До настоящего времени тысячи советских специалистов работают плечо к плечу с алжирскими специалистами в различных областях культуры и экономики страны.

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

В 1966 году была создана национальная компания текстильной промышленности " SON (ТЕХ которая должна была развивать Ал

- б жирскую текстильную промышленность путем создания новых заводов, фабрик и предприятий торговли. Основными задачами этой компании являлись: создание множества рабочих мест и удовлетворение алжирских потребителей в текстильных изделиях.

В АВДР при содействии СССР, Франции, Японии, Бельгии, Англии, Швейцарии были построены и оборудованы десятки текстильных фабрик и комбинатов.

В настоящее время в Алжире имеется 59 государственных текстильных предприятий, где трудятся около 50000 рабочих. " " производит 92 ООО тонн пряжи и III млн. различных текстильных изделий в год. Это позволяет не только удовлетворить внутренний рынок, но и даже экспортировать часть продукции.

Главная сегодняшняя проблема текстильного производства страны - нехватка ремонтных мастерских и высококвалифицированных кадров по ремонту и уходу за оборудованием, а также отсутствие научно-исследовательских центров. Поэтому данная работа "Исследование влияния посадки патронов на шпиндель прядильных веретен с целью улучшения условий их эксплуатации" является актуальной темой для развивающегося отечественного машиностроения.

Актуальность темы: несмотря на внедрение пневмомеханических прядильных машин в текстильную промышленность, кольцепрядильные машины остаются основными машинами прядильного производства. Поэтому повышение их производительности и производительности труда работающих на них рабочих являются актуальными вопросами текстильной промышленности.

Одним из основных механизмов кольцепрядильных машин является крутильно-мотальный механизм, состоящий из веретена, кольца, бегунка и раскладочного механизма.

В прядильном производстве динамика веретен сильно влияет на качество продукции и производительность труда обслуживающего персонала.

С повышением скорости вращения шпинделя веретена и объема паковки, наблюдается увеличение амплитуд колебаний шпинделя, вредно влияющих на ход технологического процесса и работоспособность веретена. Как показывает практика, выход веретена из строя в большинстве случаев происходит из-за износа опор шпинделя веретена, который зависит от динамических нагрузок в его опорах.

Динамические нагрузки в опорах в свою очередь зависят от многих факторов:

I) от точности изготовления и сборки вращающихся частей шпинделя веретена; от дисбаланса веретена и паковки;

3) от собственной неуравновешенности патрона;

4) от неправильной формы намотки ;

5) от частоты вращения шпинделя веретена;

6) от характера посадки патрона на шпиндель веретена и других факторов.

Отметим, что с увеличением вибрации шпинделя веретена растет обрывность нити.

Оуществуют несколько способов улучшения динамических характеристик веретена. Основными из них являются увеличение точности изготовления шпинделей и патронов, обеспечение правильной посадки патрона на шпиндель веретена, а также применение упругих элементов в посадочных частях насадка.

В данной работе рассматривается влияние посадки паковки на шпиндель веретена, а также влияние дисбалансов паковки в двух плоскостях приведения, угла мезду векторами дисбалансов и частоты вращения веретена на амплитуды его колебаний и реакций в верхней опоре. Для экспериментальных исследований было выбрано наиболее массовое веретено ВНГ 28-63.

Цель и задачи исследования: целью работы является исследование работоспособности веретена ВНГ 28-63 с паковкой в зависимости от его характеристик. Выяснить в процессе исследования влияние жесткости пружинных фиксаторов, статического и момент-ног о дисбалансов паковки, частоты вращения шпинделя веретена на амплитуды колебаний верхнего конца шпинделя веретена и динамических усилий в его верхней опоре.

Для решения поставленной задачи было выполнено следующее:

1. Экспериментально определены массово-геометрические параметры исследуемого веретена и паковки.

2. Экспериментально определены статическая и моментная неуравновешенности шпинделя веретена и паковки.

3. Экспериментально найдены жесткости упругой втулки шпинделя и пружинных фиксаторов.

4. Теоретически и экспериментально определены критические частоты вращения шпинделя веретена ВНГ 28-63 с паковкой.

5. Теоретически подсчитаны амплитуды колебаний верхнего конца шпинделя веретена и динамические усилия в его верхней опоре.

6. Проведен эксперимент по определению амплитуд колебаний верхнего конца шпивделя веретена и динамических усилий в его верхней опоре.

7. Амплитуды колебаний верхнего конца шпинделя веретена и динамические усилия в его верхней опоре определялись в зависимости от 5-ти факторов: а) жесткости пружинных фиксаторов б) дисбаланса в верхней плоскости приведения паковки в) дисбаланса в нижней плоскости приведения паковки г) угла между векторами дисбалансов д) частоты вращения шпинделя веретена

8. Проведен анализ полученных результатов и даны выводы и рекомендации по работе.

Общая методика исследования. В работе проведены теоретические и экспериментальные исследования. В теоретических исследованиях использован метод расчета динамики веретена по схеме полужесткого шпинделя с применением уравнений Лагранжа 2-го рода для его решения.

Теоретические расчеты проведены с использованием ЭВМ.

Достоверность данных теоретических исследований проверена путем сопоставления результатов расчета с экспериментальными данными. При этом использовались современные метода планирования эксперимента и статистической обработки результатов опытов.

Экспериментальные исследования проведены в лабораторных условиях на специальном стенде на кафедре проектирования текстильных машин Московского текстильного института.

Научная новизна работы. При проведении теоретических и экспериментальных исследований автором самостоятельно выполнено следующее:

- Разработана методика расчета амплитуд колебаний шпинделя и динамических нагрузок в его опорах с учетом жесткости пружинных фиксаторов, величин дисбалансов в верхней и нижней плоскостях приведения паковки, угла между векторами дисбалансов» податливости упругой втулки и ее приведенной массы, а также частоты вращения шпинделя веретена, составлены программы расчета на ЭЦВМ.

- Использован математический эксперимент для теоретических расчетов амплитуд колебаний верхнего конца шпинделя и динамических нагрузок в его опорах на базе матрицы плана РВ1В-схемы.

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

Практическая ценность. Предложенная методика аналитического определения влияния жесткости пружинных фиксаторов, статического и моментного дисбалансов на динамику шпинделя веретена ВНТ 28-63 при различных частотах вращения с помощью разработанных программ позволяет быстро определить критические частоты вращения, амплитуды колебаний верхнего конца шпинделя и паковки и динамические усилия в верхней опоре любого веретена с упругой втулкой,имеющего полужесткий шпиндель при разных значениях жесткости пружинных фиксаторов и дисбалансов насадка и паковки.

В работе приведены графики результатов теоретических и экспериментальных исследований динамики веретена для различных параметров системы. Эти графики представляют собой большой практический интерес для конструкторов.

Использование предложенных рекомендаций по изучению жесткости пружинных фиксаторов в посадочных частях насадка шпинделя веретена позволит улучшить работоспособность веретена.

Объем работы: диссертационная работа состоит из введения, трех глав, выводов и рекомендаций, имеет 94 странищмашино-писного текста, 52 рисунка на 52 страницах, 6 таблиц на 10

Заключение диссертация на тему "Исследование влияния посадки патронов на шпиндель прядильных веретен с целью улучшения условий их эксплуатации"

- 153 -ОБЩИЕ ВЫВОДЫ

1. Полученные формулы для расчета веретена ВНГ 28-63 дают хорошее совпадение с экспериментальными данными и могут быть рекомендованы для определения критических частот веретен с аналогичной конструктивной схемой.

2. Для анализа условий вращения веретена ВНГ 28-63 с паковкой разработана методика расчета амплитуд колебаний шпинделя и динамических нагрузок в его опорах с учетом жесткости пружинных фиксаторов, величин дисбалансов в верхней и нижней плоскостях приведения паковки, угла между векторами дисбалансов, податливости упругой втулки и ее приведенной массы, а также частоты вращения шпинделя веретена. Методика ориентирована на использование ЭВМ и может быть применена для анализа динамики аналогичных веретен.

3. Теоретически и экспериментально установлено, что существенное влияние на амплитуды колебаний шпинделя веретена и динамические нагрузки в верхней опоре шпиццеля оказывает жесткость пружинных фиксаторов. При переходе величины жесткости пружинных фиксаторов от 700 н/м до 4350 н/м, амплитуды колебаний верхнего конца шпинделя веретена с паковкой увеличились в среднем на 14,45$, а динамические усилия в верхней опоре на 12,60$, что говорит о целесообразности уменьшения жесткости фиксаторов.

4. Как показали теоретические и экспериментальные исследования, основное влияние на амплитуды колебаний верхнего конца шпинделя с паковкой и реакции в верхней опоре оказывает дисбаланс в верхней плоскости приведения паковки. Для уменьшения этого влияния желательно уменьшить жесткость пружинных фиксаторов.

5. С увеличением частоты вращения шпинделя значительно растут динамические усилия в верхней опоре шпинделя веретена. Амплитуды колебаний верхнего конца шпивделя, кроме зон критических частот, изменяются мало.

6. Теоретические значения амплитуд колебаний верхнего конца шпивделя веретена с паковкой и реакций в его верхней опоре мало отличаются от экспериментальных значений, что позволяет рекомендовать полученные теоретические зависимости для использования при проектировании веретен типа ВН.

7. Выбранный план эксперимента РВ1В -схема является эффективным планом для проведения экспериментов. Он позволяет из малого числа опытов (43) получить достаточно надежную информацию о эффектах взаимодействий пяти важнейших факторов на динамику веретена. Этот план целесообразно использовать при наличии 4-х и более факторов.

РЕКОМЩЦАЦИИ

1. Для обеспечения стабильной работы прядильных веретен типа ВН, работающих на повышенных частотах вращения, желательно применять пружинные фиксаторы в местах их посадки на шпиндели веретен.

2. Уменьшить жесткость пружинных фиксаторов, находящихся в веретенах ВНГ 28-63 от значения 1700 н/м до 5004.700 н/м и обеспечить гарантированный зазор между насадком и паковкой вне зоны фиксаторов.

3. Улучшить качество бумажных патронов, которые после относительно короткого времени использования деформируются, садятся по значительной части конуса насадка и делают, тем самым, бесполезным установку фиксаторов.

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

4. Методики определения критических частот вращения шпинделя веретена, амплитуд колебаний его врхнего конца и динамических нагрузок верхней опоры могут быть рекомендованы для использования при проектировании и исследовании веретен типа ВН.

Библиография Ушфун, Тауфик, диссертация по теме Машины, агрегаты и процессы (по отраслям)

1. Официально-документальные документы

2. I. SONiT£* -DONN££S OFTfcißtLfiS -49841. П. Книги

3. Анурьев В.И. Справочник конструктора-машиностроителя. В 3-х тт., т. I. 5-е изд., перераб. и доп. - М.: Машиностроение, 1979. - 565 с.

4. Бидерман В.П. Теория механических колебаний. М.: Высшая школа, 1980. - 407 с.

5. Бронштейн И.Н., Семевдяев К.А. Справочник по математике. -М.: Наука, 1980. 974 с.

6. Гернет М.М., Ратобылский В.Ф. Определение моментов инерции. М.: Машиностроение, 1969. - 247 с.

7. Кельзон A.C. Расчет и конструирование роторных машин. -Л.: Машиностроение, 1977. 286 с.

8. Коритысский Я.И. Исследование динамики и конструкций веретен текстильных машин. М.: Машгиз, 1963. - 643 с.

9. Коритысский Я.И. Колебания в текстильных машинах. М.: Машиностроение, 1973. - 319 с.

10. Коритысский Я.И. Вибрация и шум в текстильной и легкой промышленности. М.: Легк. индустрия, 1974. - 326 с.

11. Коритысский Я.И. Динамика упругих систем текстильных машин. М.: Легкая и пищевая промышленность, 1982. - 265 с.

12. Корнев И.В., Щукин А.й. и др. Веретена, центрифуги, прядильные камеры текстильных машин. М.: Легк. индустрия, 1978. - 128 с.

13. Макаров А.И. и др. Расчет и конструирование машин прядильного производства. М.: Машиностроение, 1969. - 511 с.

14. Макаров А.И. и др. Основы проектирования текстильных машин. М.: Машиностроение, 1976. - 416 с.

15. Макаров А.И. Расчет и конструирование машин прядильного производства. М.: Машиностроение, 1981. - 384 с.

16. Малышев А.П. Веретено. М.: Гизлегпром, 1950. - 239 с.

17. Мальмберг К.Ё. Влияние главнейших параметров паковок на работу веретен в прядении и кручении. М.: Легк. индустрия, 1965. - 134 с.

18. Маркова Ё.В., Лисенков А.Н. Планирование эксперимента в условиях неоднородностей. М.: Наука, 1973. - 218 с.

19. Маркова Е.В., Лисенков А.Н. Комбинаторные планы в задачах многофакторного эксперимента. М.: Наука, 1979. - 348 с.

20. Маслов Г.С. Расчеты колебаний валов. Справочник. М.: Машиностроение, 1980. - 151 с.

21. Моменты инерций тел. Справочник. М.: Машиностроение, 1977. - 511 с.

22. Николаи Е.Л. Теория гироскопов. М.: Гостехиздат, 1948. -142 с.

23. Основы балансировочной техники / Под ред. В.А. Щепетильни-кова. М.: Машиностроение, 1975. T. I. Уравновешивание жестких роторов и механизмов. 1975. - 526 с.

24. Пановко Я.Г. Основы прикладной теории колебаний и удара. -Л.: Машиностроение, 1976. 319 с.

25. Пановко Я.Г. Введение в теорию механических колебаний. -М.: Наука, 1980. 270 с.

26. Писаренко Г.С. и Др. Справочник по сопротивлению материалов. К.: Наукова думка, 1975. - 703 с.

27. Пономарев С.Д. и др. Расчеты на прочность в машиностроении. В 3-х тт., т. 3. М.: Машиностроение, 1959. - III8 с.

28. Попов Э.А. и др. Основы теории конструкции и расчета текстильных машин. М.: Машиностроение, 1975. - 388 с.

29. Севостьянов А.Г. Методы и средства исследования механико-технологических процессов текстильной промышленности. -М.: Легкая индустрия, 1980. 392 с.

30. Тимошенко С.П. Колебания в инженерном деле. Перевод под рук. Пановко Я.Г. М.: Наука, 1976. - 444 с.

31. Тихомиров В.Б. Планирование и анализ эксперимента. М.: Легкая индустрия, 1974. - 263 с.

32. Туричин A.M. Электрические измерения неэлектрических величин. М.: Энергия, 1975. - 600 с.

33. Феодосьев В.И. Сопротивление материалов. М.: Наука, 1972. - 544 с.

34. Филиппов А.П. Колебания деформируемых систем. М.: Машиностроение, 1970. - 728 с.

35. САсд/еибед, /., билеte, сИм deífcVioitaAM. mdUwtUel ,А ЫЬ,еЙ1/1АлЛл€ HadMte .235. §)J!u¿ Р. (Ый&нли. ote* таХ/1мош . lecWcyae Jr1. Ш. Статьи

36. Жуковский Н.Е. Об упругой оси турбины Лаваля и об осях с качающимися подшипниками. М., 1899. - 22 с.

37. Кущуль М.Я. Приближенный метод определения критических скоростей многоопорного вала переменного сечения. Сб. № 2. "Поперечные колебания и критические скорости". Изд. АН СССР, 1953.

38. Макаров А.И. Значение работ проф. Н.Е.Жуковского для изучения вопросов устойчивости вращения веретен при рабочих скоростях. Сб. научно-исслед. тр. / МТИ, т. 15. - М.: Гиз-легпром, 1955. - 21-32 с.

39. Степанов С.Г. Исследование и проектирование фрикционных намоточных механизмов машин для производства химических волокон.

40. Фрид И.А. Исследование нестационарных колебаний веретен прядильно-крутильных машин при анизотропных характеристиках жесткостей опор. Технология текстильной промышленности, 1979, № 2, с. 125-131.

41. Шнайдер Г.З. и др. Исследование колебаний и данамические расчеты новых конструкций электроверетен ЗВН-2 (обзор). -М., 1971.

42. Щукин А.И. Влияние упругости пружинных фиксаторов на динамику веретен текстильных машин. ЦНИИТЭИЛегпищемаш. М. Машиностроение для текстильной промышленности. Реферетивный- 160 сборник. Вып. П, 1978. 8-II с. 1У. Диссертации

43. Акимов A.A. Исследование влияния неуравновешенности бобино-держателя маятникового типа и биения поверхности тела намотки на динамические характеристики. Канд. диссертация. М., 1983.

44. Артюнин А.И. Исследование динамики высокоскоростных электроверетен тяжелого типа. Канд. диссертация. М., 1979.

45. Квартин Л.М. Исследование и расчет основных конструктивных параметров веретен прядильных и крутильных машин. Канд. диссертация. М., 1981.

46. Кучеров Б.К. Динамика кольцепрядильного веретена для больших паковок с невращающимся шпинделем и гасителем вибрации. Докторская диссертация. М., 1958.

47. Макаров А.И. Вибрация шпиеделей кольцепрядильных и кольце-крутильных веретен. Канд. диссертация. М., 1936.

48. Макаров А.И. Применение центрифуг с опрокинутыми кружками для выработки пряжи из хлопка и шерсти. Докторская диссертация. М., 1952.

49. Папоян А.Р. Разработка метода проектирования и расчета опор катушек для больших паковок крутильных машин. Канд. диссертация. М., 1982.

50. Полыциков В.И. Влияние собственной неуравновешенности бумажных армированных патронов на амплитуду колебаний шпинделя веретен типа ВНГ при рабочей скорости. Канд. диссертация. -М., 1970.

51. Попов Э.А. Нестационарные колебания веретен с гибким шпинделем и тяжелыми насадками. Канд. диссертация. М., 1964.

52. Попов Э.А. Динамика шпинделей крутильных механизмов текстильных машин. Докторская диссертация. М., 1977.

53. Эль-Бакри М.С.И. Некоторые вопросы динамики веретен коль-цепрядильных машин. Канд. диссертация. М., 1965.

54. Эль-Хавари И.А. Исследование влияния посадки шпуль на вибрацию крутильного веретена. Канд. диссертация. М., 1978.- 162 -IIpEJioseHHe № I

55. DZ1<D=DZ1<9> DZlU7>=l.E-05 DO 14 1=24/31

56. DZia>=DZK17> DO 15 1=36/43

57. DZKI>=DZia7> DZ1<20>=DZK1> DZK21)=DZl(i> DZ1<22>=DZ1<9> DZi <23>=DZ1<9> DZ1< 32 >=DZ1<1) DZK33>=DZia> DZI<34>=DZ1<9> DZI<35>=DZ1C9> READ(1/101> DZ2 READ<1/202> ALFA READa/202> OM31 FORMAT 20F3.1>

58. FORMATS14F5.0) DO 22 1 = 1/431. KF2< i >=k:fi a >«1000.

59. ALFA i I > =3.14«ALFA (. I > 180. DO 2 1 = 1/43 DZ2 CI > =DZ2 < I > 100000.

60. DZ=SQRTCDZ1 < I >»=s-:2+DZ2( I >**:2+2. sDZl < I >«DZ2< I >:*COS<ALFA< I > > >1—1a 1-1 1—1 c"j 1.

61. A LUл. ж' 1—4

62. И™'. Ы OJ •w- +1—1 Ж C'J C'J•X Ж Ll Ж1.. ьн p- *-.;'' ж1.. —I V J Ж •X-J •X •X + C'J

63. Z 1—1 Ll жto —I 1—1 ж KliTl .u. О •:X iTi TT*о о 1-4 J ■+1 о в. ж —r* Ll + C'J x¡1 ж 1—1 JUL. to l C'J

64. ГГ\ "M 1—1 Ь-1 to 1—1 ж 1—4 Cl ж to

65. СМ P-J ID P-J Jú W: C'J Ж C'J C'J P- 1—4 C'J -•■ C' J ••-•• CO J i—t Ж !i -í- ö я i—i i + P-- Ll P- LL Ж J Ll CM Ll C'J J + Ж LU

66. CM е. O 1— © i—i J -y*' l Ж 1 1-4 JJL, Ll •jiilj. 1 C'J SIL. ■M" Pж CL О I- 00 — 1 Ж 1 ж •-- --• Ll Ж ;-,••• Ж .--■• Ж l~l J Ж

67. C i—i ь- !i i—i 1—I t—l + iTi p- 1—11—1 P-- Ж iTi Ж i—1 Ж ж C'J

68. CM 1 •■„■■ 1 © C'J i—i i—i J •X i—i j ■■-•• -- n + l iTl J со •-• C'J + ж

69. J —! ••-•• О "tf- h-Ч 1-4 1—¡ PO Kl Li. 1 i 1 1 1 1 -•• 1 -- 1 1 1 J © Ж i l 1 "-• © О © ©T Cl

70. PM II 1! 0"' <r 1— ll J •X •X O G" (Ti О -X -X 1-1 -X Ll »•. a. a*. a. a. в. a. a. a. a. a. a. a. a. a., a., a., a. a. a. a., a., .•-••

71. Pi i-i O-J II I- —т* •:Х -Œ •X •x II II 1— II II h- и и h-s: Ii 1-4 1—1 1—1 1—1 1—1 ■t-i C'J C'J C' J C'J C'J C' J Kl Ki Kl Ki Ki Ki -nj- Ti- -if -íf rj" Tj- LO Ll"i Li~i Li~i Li"> C'J K' "tf

72. N P-J P-J J I—1 Ll. 1—1 C'J 1-й C'J -■•-i C'J 2: ií: to ••-■• •■„■' ••„•' V V ••«.♦ •w* ■-„•• V V V -—•* •*

73. DO 3 12=1/5 3 R2CI1/I2)=RKI1/12)

74. CALL GEL6 < 61R1 / 5/ 1 / 1. 8E-86/1 PR ) IF':! I PR. NE. 0 ) PR I NT 182/ I PR 32 FORM AT a e^-' IPR=-'/I2)

75. CALL 6EL6<62/ R2/ 5,-1. 8E-86/ I PR) IF<IF'R.NE.8) PRINT 182/IPR

76. FORMAT (18!::!/ A1 = •' / 612.5/ 2;:i/ •' A3= •' / G12.5/ 2ft' •' A5= / 612.5/ 2X/ ' A7=' / «612.5/ 2X> •'A9=''612.5)

77. FORMATi 18^ Bl-f / 612.5/ 2Xr ' B3=-' /612.5/ 2X>' B5=-' / 612.5 «/2!K/ -'B7=J/612.5/ 2'yif ''B9=-'/612.5)

78. A1=61<1) A3=61<2) A5=61 < 3) Bl=62<1) B3=62 <2) B5=62i'3':> B7=62<4) A7=61 c'4> A9=61(5) B9=62 <5)

79. R1 R= < < -B1 +L4'«B9-"A) ( M1 *L4-M2 > + ■:.' -B5+B9,-'A) * < -M2*L4+M3 ) s+OM < I) ma2&B?x < < AF'-CP ) -MP*:L8» (L8+B ) ) ) --'A RR1<I)=SQRT <R1N$$2+R1R*s2)

80. R2H= C < A1-L4'« A9--'A ) s (-M1 £ < L4+A ) +M2 ) + < A5-A9/A ) s < M2*: < L4+A ) :*-M3 ) + A7:*0M < I) ( < AF'-CF') -MPsLS® ( L8+B+A > ) ) .-'A

81. R2R= (< -B1 +L4$B9--'A ) ( -M1 ( L4+A ) +M2 ) + < -B5+B9.--A )< M2# < L4+A ) s-M3 ) +OM < I) $:*2$B7$ < MF'KL8$ < L8+B+A) (AF'-CF') ) ) •••"A RR2< I )=Si3RT<R2N»:s2+R2R»;'*2)

82. FORMAT (18IK/ •' KF=/ F5.8/ DZ1 = / 618.3/ •' DZ2= •' / 618.3/ * . ALFA= •' / F8.5/' OM=' / F5.8 >

83. PRINT 185/AK/AF/RR1(I)/RR2(I) 2 CONTINUE

84. FORMAT< 1AK= / 612.5/ •' AF=-"/612.5/RR1 = ' / 612.5/1. RR2=-'/G12.5) STOP END

85. SUBROUTI NE GEL6 < R/ A/ M/ N/ EPS/1ER ) DIMENSION A<1)/R<1) IF<M)23/23/l 1 IER=8 F'IU=8. MM=M*M NM=N*M DO 3 L=l/MM TB=hBS(A(L >)1.aB-FIU>3.'3,22 PIU=TB I=L

86. CONTINUE TOL=EPS'«PIU LST=í

87. DO 17 K=1.'M I FC PH.i-3 27,23'44 IF<IER>7.'5.'75 IF <PIU-TOL>6.' б." 76 IER=K-1

88. PIUI=1./h(I) J=< I~1 >.--м I=I-J*M-K J=J+1-K1. DO 8 L=K.f NM/ M LL=L+I

89. TB=PIUI$R<LL> R<LL::'=R<L>8 R<L>=TB lF<K-m9r 18.' 18

90. LEND=LST+M-K IF< J> 12.' 12.'1010 II=J*:M

91. DO 11 L=LST.' LEND TB=R<L> LL=L+II n<L>=h(LL>11 A(LL>=TB12 DO 13 L=LST.'MMj-M LL=L+I1. TB=PIUI-H(LL) n(LL>=H<L)

92. A<L>=TB A<LST>=J PIU=0. LST=L8T+1 , l=ñ

93. DO 16 I I=LST.' LEND PIUI=-A<11> I8T=II+M J=J+11. DO 15 L=IST.'MM.'M LL=L-J

94. A(L>=h<L> +pIUIжA i LL) TB=ABS<A<L)> IF(TB-PIU) 15.' 15.' 1414 PIU=TB I=L15 CONTINUE1. DO 16 L=K/ NM/ M LL=L+.J

95. R('. LL > =R (. LL > +PIUI sR <L >17 LST=LST+M18 IF < M-1') 23/ 22/ 1919 IST=MM+M LST=M+11. DO 21 1=2/M1.=LST-I1.T=IST-LST1.IST-M1.A (. L > +0.51. DO 21 J=II/NM/M1. TB=R<J>1.=J1. DO 20 K=IST/MM/M LL=LL+1

96. TB=TB-ft < K > MR < LL > K=J+L ROJ>=R<K>21 R':'K":'=TB22 RETURN23 IER=-1 RETURN END

97. KF=255@. nZl= 8.158E-04 BZ2= 0.150E-04 ALFñ= 1,55256 OM=1050.

98. Al= 0.29593Е-06 АЗ=-8.14214Е-03 А5= 0.26707Е-05 А7=-0.119S7E-82 А9= 0.39112Е-07

99. В1= 0.11125Е-05 ВЗ=-0.12017Е-04 В5= 0.14879Е-04 В7=-0.98861Е-04 В9=-0.67460Е-06 АК= 0.40811Е~05 AF= 0.30502Е-03 RR1= 20.816 RR2= 38.301 KF=2550. DZ1= 0.150E-04 DZ2= 0.150E-04 ALFA= 1.55256 0M= 630.

100. Al= 0.58211E-05 A3=-0.15618E-03 A5= 0.75142E-84 A7=-0.13288E-82 A9= 0.85736E-05

101. Bl= 0.16530E-05 B3=-0.13314E-04 B5= 0.21622E-04 B7=-0.11527E-83 B9= 0.25695E-05 AK= 0.21220E-04 AF= 0.33671E-03 RR1= 8.3027 RR2= 15.322 KF=2550. DZ1= 0.150E-04 HZ2= 0.500E-05 ALFA= 1.55256 0M= 630.

102. Al= 0.43232E-05 A3=-0.11600E-03 A5= 0.55806E-04 A7=-0.98630E-03 A9= 0.63675E-05

103. Bl= 0.12277E-05 B3=-0.98878E-05 B5= 0.16058E-04 B7=-0.85609E-04 B9= 0.19083E-05 нК= 0.15760E-04 AF= 0.25007E-03 RR1= 6.1663 RR2= 11.380 KF=2550. DZ1= 0.150E-04 HZ2= 0.500E-05 ALFA= 1.55256 0M=1858.

104. Al= 0.21978E-06 A3=-0.10556E-03 A5= 8.19834E-85 A7=-0.89022E-03 A9= 0.29048E-07

105. В1= 0.82627E-06 B3=-0.89247E-05 B5= 0.11851E-04 B7=-0.73422E-84 B9=-0.50101E-06 AK= 0.30310E-05 AF= 0.22653E-83 RR1= 15.460 RR2= 28.445 KF=2550. DZ1= 0.150E-04 DZ2= 0.100E-04 ALFA= 3.12256 0M=1858.

106. Al= 0.69199E-07 A3=-0.33237E-04 A5= 0.62450E-06 A7=-8.28829E-83 A9= 0.91457E-08

107. Bl= 0.37032E-06 B3=-0.39999E-05 B5= 0.49526E-05 B7=-0.32906E-04 B9=-8.22454E-86 AK= 0.13447E-05 AF= 0.71576E-04 RR1= 4.8843 RR2= 9.0243 KF=2550. DZ1= 0.150E-04 DZ2= 8.108E-04 ALFA= 3.12256 0M= 630.

108. Al= 0.13612E-05 A3=-0.36522E-04 A5= 0.17571E-04 A7=-0.31054E-83 A9= 0.28048E-05

109. Bl= 0.55022E-06 B3=-0.44316E-05 B5= 0.71971E-05 B7=-0.38368E-04 B9= 0.85526E-06 AK= 0.51518E-05 AF= 0.79031E-04 RR1= 1.9489 RR2= 3.6103 KF=2550. DZ1= 0.150E-04 DZ2= 6.188E-04 ALFA= 0.01744 0M= 630.

110. Al= 0.67982E-05 A3=-0.18240E-03 A5= 0.87755E-04 A7=-8.15518E-82 A9= 8.18813E-84

111. Bl= 8.19775E-86 B3=-8.15927E-85 B5= 8.25866E-85 B7=-8.13789E-84 B9= 8.38738E-86 AK= 8.23833E-84 AF= 8.39188E-83 RR1= 9.6686 RR2= 17.761 KF=2558. DZ1= 8.150E-04 CZ2= 0.100E-04 ALFA= 0.01744 OM=1050.

112. Al= 8.34561E-06 A3=-0.16600E-03 A5= 8.31198E-05 A7=-8.13999E-82 A9= 8.45677E-07

113. Bl= 0.13309E-06 B3=-0.14375E-05 B5= 0.17800E-05 B7=-0.11826E-04 B9=-0.80700E-07 AK'= 0.10643E-05 AF= 0.35500E-03 RR1= 24.230 RR2= 44.400 KF=2550. DZ1= 0.500E-05 DZ2= 8.158E-84 ALFA= 1.55256 0M=1858.

114. Al= 8.21978h~86 A3=-8.18556E-83 A5= 8.19834E-85 A7=-8.89822E-83 A9= 8.29848E-87

115. В1= 8.82627E-86 B3=-8.89247E-85 B5= 8.11851E-84 B7=-8.73422E-84 B9=-8.58181E-86 AK= 8.38318E-85 AF= 8.22653E-83 RR1= 15.468 RR2= 28.445 KF=2550. DZ1= 0.500E-05 DZ2= 0.150E-04 ALFA= 1.55256 0M= 630.

116. Al= 0.43232E-05 A3=-0.11600E-03 A5= 0.55806E-04 A7=-0.98630E-03 A9= 0.63675E-05

117. Bl= 0.12277E-05 B3=-0.98878E-05 B5= 0.16058E-04 B7=-0.85609E-04 B9= 0.19083E-85 AK= 0.15760E-04 AF= 8.25887E-83 RR1= 6.1663 RR2- 11.380 KF=2550. DZ1= 0.500E-05 DZ2= 0.508E-05 ALFA= 1.55256 OM=1050.

118. Al= 0.98644E-07 A3=-0.47380E-04 A5= 0.89022E-06 A7=-0.39955E-03 A9= 0.13037E-07

119. Bl= 0.37085E-06 B3=-0.40056E-05 B5= 0.49597E-05 B7=-0.32954E-04 B9=-8.22487E-86 AK= 8.13684E-85 AF= 8.18167E-83 RR1= 6.9388 RR2= 12.767 KF=2558. DZ1= 8.588E-85 DZ2= 8.588E-05 ALFA= 1.55256 0M= 630.

120. Al= 0.19404E-05 A3=-0.52062E-04 A5= 0.25047E-04 A7=-0.44268E-03 A9= 0.28579E-05

121. Bl= 0.55101E-06 B3=-0.44379E-05 B5= 0.72074E-05 B7=-0.38423E-84 B9= 0.85649E-06 AK= 0.70734E-05 AF= 0.11224E-03 RR1= 2.7676 RR2= 5.1874 KF-2558. DZ1= 8.588E-85 DZ2= 8.180E-04 ALFA= 3.12256 OM=1050.

122. Al= 0.69149E-07 A3=-0.33213E-04 A5= 0.62404E-06 A7=-0.28009E-03 A9= 0.91392E-08

123. Bl= 0.37005E-06 B3=-0.39970E-05 B5= 0.49491E-05 B7=-0.32883E-04 B9=-0.22438E-06 AK= 0.13437E-05 AF= 0.71524E-04 RR1= 4.8808 RR2= 9.0178 KF=2550. HZ1= 0.500E-05 DZ2= 8.188E-84 ALFA= 3.12256 0M= 638.

124. R= IPR= Al = Bl = AK= KF= IPRS I PR= Al = Bl = AK= KF= Al=-Bl = AK= KF= Al=-Bl = AK= KF= Al=-Bl = AK= KF= Al=-Bl = AK= KF= Al=-Bl = AK= KF= Al=-Bl = AK= KF= Al=-Bl =4 = 4 8.88888 8.08800 0.0ЙЙ00 =2550. DZ44 0.88880 8 81 =

125. A3= 8.88888 B3= 8.ñññññ AF= 8.88888 5.800 DZ2=1. A5=1. B5= RR1 = 0.800Eи.инииУ A7= и.йииии0.88888 В7= 8.880800.00008 RR2= 8.88888 •84 ALFA= 8.81744 ÜM= 8.1. А5=

126. ALFA= 8.81744 -8.18479Е-84 А7^8.16122Е-85 В7' 15.255 RR2=

127. А9= 8.37878Е-85 В9=-0.30i36Е-86

128. А9= 8.74185Е-86 В9=-8.83792Е-86

129. I'- ■ I'- * X'- ■ I:- " I'- ■ " X- * I> " X- " • X " I:- ■ X> ■

130. II II II 11 II II II II II II II II 11 11 11 11 II 11 11 11 11 11 II II II

131. B3=-0.81548E-05 B5= 0.17380E5.0843 05 ALFA=8.15185E-04 A7

132. Bl= 0.68842E-06 B3=-0.65575E-05 B5= 0.91630E-05 B7=-0.55540E-04 B9= AK= 8.48482E-05 AF= 0.16636E-03 RR1= 11.355 RR2= 28.914-0.40996E-06- 171 -npHJiOKeHHe IS 21.PLICIT REALOG Lr M)

133. R1 CI/1 >=-MSH»:0M*:s2+Ml+KFl+KF2 R1 <1/2 > =-KF1-KF2

134. R1 (1/ 3>=-M2+Li 1*KF1H.7*KF2

135. CALL LADATF ( R1, R2/ 5/ 5/ 5/ D1 / D2/ G1 / G2/ G3/1ER > IF<IER.NE.0> WRITE<5/200> IER :00 FORMATS 10H/ •' IER=-"/ II > D='::2.s®D2>*Dl WRITE(5/102) D1. < ABS < D . GT. 1. E+12 > D=D-" 1. E+12 WRITEC5/182) D 02 FORMAT < 10!::!/ D= •' / G12.5 >1. IF IF

136. ABS < D>„LE„1.06E-06> GO TO 7 •J. EQ, 1 > GO TO 3•::'::D«DPRED>/<ABSCDPRED>>>.LE.0. > GO TO 69 ;om=om+domega

137. FORMAT <i 8Х/'DOMEGA=',618.3) WRITERS/184) DOMEGA DPRED=D GO TO 1

138. OM=OM~DOMEGA DOMEG A=D0MEGA-"'2.1.<DOMEGA.LE.8.85) GO TO 71. DPRED=D1. GO TO 17 PRINT 183/OM

139. FORMAT < 18л/ •' ЧАСТОТА СОБСТВЕННЫХ КОЛЕБАНИЙ ж СТЕР1+1НЯ OMEGA=' / G 15.8 > STOP END

140. SUBROUTINE LADATF < A.- Lib Н/ I А/1DGT/ D1 / D2/ I PUT/ EQIJIL* U А/ IER >

141. DIMENSION REAL DATA IER=8 RN=N

142. WREL=ZERO Dl=ONE D2=ZER0 BIGA=ZERO DO 18 1=1/N BIG=ZERO DO 5 J=l/N P=A< I/ J) LIKI/J)=P P=ABS<P>

143. CP .GT. BIG) BIG=P CONTINUE IF (BIG .GT.BIGA) BIGA=BIG IF <BIG .EQ.ZERO) GO TO 118 EQUILC I )= ONE.--'BIG 18 CONTINUE

144. A(1А/1)/LIK1А/1)/IPUTCD/EQUILCl) LU

145. ZERO/ ONE/ FOUR/ SIXTH/ SIXTH^'8.8/ 1. / 4. /16. / . 8625.-'51. DO 185 J=l/N •JM1=J-1

146. CJM1 .LT. 1) GO TO 48 ВЫЧИСЛЕНИЕ U •:: I / J ) / 1 = 1/./ J-1 DO 35 1=1/JM1 SUM = LU<I/J) IM1 =1-11. aDOT .EQ. 8) GO TO 25 ПРОВЕРКА НА ТОЧНОСТЬ РАЗЛОЖЕНИЯ AI =ABS(SUM) bJI=ZERO

147. CIMl .LT. 1) GO TO 28 DO 15 К=1/IM1

148. T = LiJC I/ K):#L¡ К К/ • I)1.ti'ÍHscT)15 CONTINUE

149. CI.-J > =SUM 20 WI =WI+ABSCSUM)1. <AI .EQ. ZERO) AI=BIGA1. TEST = WI.--AI

150. СTEST .ST. WREL) UREL=TEST GO TO 35

151. БЕЗ ПРОВЕРКИ НА ТОЧНОСТЬ РАЗЛОЖЕНИЯ 25 IF СIMl .LT.1) GO TO 35

152. DO 30 K=i.' IM1 Si IM=SUM-LU С Ь К ) *LU С К/ J ) 30 CONTINUE1.CI.'J)=SUM 35 CONTINUE

153. ВЫЧИСЛЕНИЕ U С J, J ) И L СI , J ).' I = J+1.'.' N 40 P= ZERO

154. DO 70 I=J.'N SUM =LUC I.' J)

155. CIDGT .EQ. 0) GO TO 55 ПРОВЕРКА НА ТОЧНОСТЬ AI =ABSCSU!i)1. Ш =ZERO

156. CJM1 .LT. 1) GO ТЛ m DO 45 K=1.'JM1 T=LUCI.'K)»:LUCK'.I) SUM=SUM-T1. WI=WI+ABSCT) '5 CONTINUE

157. C I.' J)=SUM >0 I.J I =i.,.ii +ABSCSUM)

158. CAI .EQ.ZÉRO) AI=BIGA TEST=WI,-AI

159. С TEST . GT.WREL) I.IREL =ТЕЯТ GO ТП

160. БЕЗ ПРОВЕРКИ НА ТОЧНОСТЬ 55 IF (JM1 .LT. 1) GO ТП1. DO 68 K=1.'.JM1

161. SUM =SUM -LUC I.' K)$Li iCK' ■ i) 60 CONTINUE

162. C I.'J) =SUM Q = EQUIL СI) :+:ABS < SUM ) IF CP .GE. Q) GO Tfl 70 P=Q IMAX=I P0 CONTINUE

163. ПРОВЕРКА НА АЛГОРИТМИЧЕСКУЮ СИНГУЛЯРНОСТЬ IF CRN+P .EQ. RN) GO TO 110 " '1. <J .EQ. I МАЮ GO TO SO ПЕРЕСТАНОВКА СТРОКИ J С IМАИ D1=-D1 DO 75 К=ЬН

164. P= LUCIМАХ/К) LLKIMAK.-Ю =LU<J/K) LLK J/ K)=P 75 CONTINUE

165. EQiJIL< IMAX>= EQUILCJ) 88 IF'UTOJ) = IMAX

166. DI = Dl:*L!J(J/J> 85 IF САВЗШ> .LE. ONE) GO TO 98 DI =D1*SIXTH D2 =D2+F0UR GO TO 8598 ÍF <ABS<D1) .GE. SIXTH)GO TO 95 Di = DlstíSIXTN D2 = D2-F0UR GO ТП 98 95 CONTINUE JPl =J+11. <JP1 .GT. N) GO TO 185 ДЕЛЕНИЕ НА ВЕДУЩИЙ ЭЛЕМЕНТ IKJz J)

167. F-LLK.ii" J> DO 188 I=JF'l.' N1.a/J)=LUc:i/J).-P 88 CONTINUE 85 CONTINUE1. <IDGT .EQ. 8) GO TO 9885p=3*N+3l.JA = P*WREL1. < 1 у „ ЙЖЖ < -1 DGT ) .NE. ЫА) GO TO 9885 IER-34 GO TO 988'й 18 IER =129 D1 =ZERO D2 =ZERO38 CONTINUE 35 RETURN END

168. ЧАСТОТА СОБСТВЕННЫХ КОЛЕБАНИЙ СИСТЕМЫ "ВЕРЕТЕНО-ПАКОВКА" ОМЕГА=119.21875

169. ЧАСТОТА СОБСТВЕННЫХ КОЛЕБАНИЙ СИСТЕМЫ "ВЕРЕТЕНО-ПАКОВКА" ОМЕГА=234.53125

170. ЧАСТОТА СОБСТВЕННЫХ КОЛЕБАНИЙ СИСТЕМЫ "ВЕРЕТЕНО-ПАКОВКА" ОМЕГА=8882.8313- 176 -EtpHJiojKeHHe № 3

171. F'L ICIT REAL < J/ K/ L/ M>

172. KF1<I>=KFi <20> DO 111 1=33/43

173. KF1<I>=KF1<32> DZKl>=1.5E-05 DO 12 1=2/8

174. DZKI>=DZ1<1> DZ1< 9 > =0.5E-05 DO 13 1=10/16

175. DZ1<I>=DZ1<9> DZ1<17>=1.E-05 DO 14 1=24/31

176. DZKI>=DZ1<17> DO 15 1=36/43

177. DZ1<I>=DZ1<17> DZ1<20>=DZ1<1> DZK21>=DZ1<1> DZ1<22>=DZ1<9> DZ1<23>=DZ1< 9 > DZ1<32>=DZ1<1> DZ1<33>=DZ1<1> DZ1<34>=DZ1<9> DZ1<35>=DZ1<9> READ<1/101> DZ2 READ<1/202) ALFA READ<1/202) OM

178. FORMAT<20F3.1) i82 FORMAT<14F5.0> DO 22 1 = 1/43 KF2<I)=KFKI> 22 ALFA<I>=3.14«ALFA<IV180. DO 2 1 = 1/43 DZ2 < I > =DZ2 (I VI 00000.

179. DZ=SQRT < DZ 1 < I) »:.*'2+DZ2 < I > ««2+2. «DZ 1 < I) *DZ2 < I) «COS < ALFA < I) ) )1. E2=DZ.-'MP1. J0Z1=DZ«L11. J0Z2=DZ«L2

180. JOZ=SQRT< J0Zl«:*2+J0Z2:*:«2-2. «J0Z1«J0Z2*C0S<ALFA< I > ) >

181. DELTA2=J0Z/ С 9.81 * < CP-AP ) )

182. PR I NT 104/ KF ICI)DZ KI) / DZ2 < I ) ALFA С I > , ОМ С I )1.(ALFA <I).EQ.8.) GO TO 41. СALFAСI).EQ.90.) GO TO 51. <ALFA <I).EQ.1SO.> GO TO 4

183. ALFA I > =3« 14$ALFA С I )/180.

184. F11=ÄTÄN(-TAN <ALFA(I))/С DZ1<I)/< DZ2СI)*COS<ALFA(I))) +1.)) F12=ATAN С 1. /TAN < ALFA ( I ) ) JOZ1 / С J0Z2®S IN С ALFAСI ) ) ) > GO TO 64 FI1=0.

185. FI2=90.$3.14/180. hfi ТП f.

186. F11=ATAN С-DZ2СI)/DZ КI )) F12=ATAN С-JOZ1/J0Z2 >6 CONTINUE

187. FORMATOS 11610.3) EF'S=F 11 —FI

188. R К bl ) =-MSH*OM < I ) ®*2+M 1+KFKI) +KF2 С I ) RКЬ2)=-KFКI)-KF2<I) RК U 3) =-M2+Lll:*KFKI)-L7&KF2СI ) R1 С 1, 4)=-L12*:KFK I )+L8i*KF2C I ) R К1 5 )=< -M 1 *L4+M2 ) / А GlС1)=0.

189. R К 2/ 1 ) =-M2+L 11 *KF KI) -L7«KF2 С I )

190. RК2/2)=-L11*KFКI)+L7*KF2<I>

191. RI ( 2? 3>=<r:SH-ASH>'«OMí I )'*"«2+M3+Ll 1**2*KF1 С I )+L7'«'*2'«KF2( I > RI С 2.' 4 > =-L 11 12-«KF Kl) -L7®L8sKF2 С I ) R1С2r5) = СM2$L4-M3)/А GK2)=0.

192. R1 С 3.' 1 ) =-KF Kl) -KF2 С I )

193. R1 < 3.' 2 ) =-MF':*OM С I ) #:#2+KF KI) +KF2 ( I )

194. R К 3/ 3 ) =-L 11 *KF КI ) +L7*KF2 С I )

195. R1 С3r 4) =L 12*KF 1 СI)-L8*KF2C I )1. R1C3.' 5)=0.1. Gl С 3) =МРаС2*ЮМСI )

196. RI С4.' 1>=-L12*KFKI>+L8*KF2CI)

197. RI С4.' 2)=L12:*KF1 СI )-L8*KF2C I )

198. R1С 4.' 3 ) =-L 12*:L 11 *KF ICI) -L8*L7:*KF2 СI >

199. RI С 4.» 4 )=С CP-АР > '*OM СI )**2+L12**2*KFl СI >+L8**2:*KF2C I )1. R1 C4.f 5)=0.

200. Gl С 4 > =- С CP-AP ) *DELTA2*0M СI ) IN С EPS2 )

201. R1С 5.' 1 ) = С -L4:*M 1+М2 ) /А1. R1С 5.' 2)=8.

202. R1С 5.' 3 )=С M2$L4-M3 ) / А1. R1C5.'4)=0.

203. R1С 5.' 5 ) =-МЕ::*ОМ СI ) С L4**2*M 1 -2. *M2*L4+M3)/A5 (Ш2+КВ1. G1C5)=0.1. G2C1)=0.62С2>=0.1. G2C3)=0.

204. G2 С 4 ) = С CP-АР ) жОМ СI ) Ä»2:*DELTA2:*C0S С EPS2 )

205. G2<5>=8. DO 3 11=1/5 DO 3 12=1/5 3 R2CI1/I2>=Riai/ 12)

206. CALL GELG C G1 / R1 / 5/ 1 / 1. £€-86/1 PR ) IFaPR.HE.8> PRINT 182/ I PR 82 FORMAT< 18!K/ ■' IF'R="' /12)

207. CALL GELG<G2/R2/5/l/1.8E-86/IPR) IFCIPR.NE.8) PRINT 182/IPR

208. FORMAT < 18>-,/ A1 -' / G12.5/ 2X/ ' A3=/ G12.5/ 2K- A5= •' / G12.5/ 2X/ ' A7= / «G12.5/ 2Xr A9="' / G12.5)

209. FORMAT-; 18X/ •' Bl = ' / G12.5/ 2X/ ' B3='/ G12.5/ 2K/ ' B5=-' / G12.5 :*/2X/ B7=-'/ G12.5/ 2X/ 'B9=' / G12.5)

210. A1=G1C1> A3=G1<2) A5=G1C3> B1=G2C1) B3=G2<2) B5=G2<3) B7=G2<4> A7=G1<4) A9=G1C5> B9=G2C5>

211. R1 R= < C -B1 +L4:*B9.-'A ) < M1 ®L4-M2 )+C -B5+B9/A ) * C -M2«L4+M3 ) s+OM < I > wm2&B7'# C < AF'-CP > -MP*:L8:* < L8+B ) > ) .--'A RR1 < I) =S0RT C R1 N«:*2+R 1 R*:«2 >

212. R2N= C C A1 -L4«A9--'A ) ® C -M1 * (L4+A ) +M2 > + C A5-A9.--A > $ C M2« C L4+A > s-M3 ) + A7«0M < I) '««2« < C AF'-CF') -MP*L8$ < L8+B+A ) ) ) / A

213. R2R= < C -B1 +L4«B9/A ) $ C -M 1:*CL4+A > +M2 ) + C -B5+B9/A > C M2$ C L4+A ) ;«-M3 ) +OM CI) :*:.*2«B7:* CMF'sLSfc" C L8+B+A > C AF'-CF') ) ) / A RR2 C I > =SQRT < R2N$*2+R2Rs$2 )

214. FORMAT < 18X/ KF=' / F5.8/ •' DZ 1 -' / G18.3/ DZ2= ' / 618.3/ •' ALFA= / F8.5/ •' OM=•' / F5.8 >

215. PRINT 185/AK/AF/RR1CI)/RR2CI> 2 CONTINUE

216. FORMAT< 18>i/ ' AK=' / G12.5/ AF=' / G12.5/ ' RR 1 -' / G12.5/1. RR2=-'/G12.5> STOP END

217. SUBROUTINE GELG <R/A/M/ N/EPS/IER > DIMENSION AC 1)/RC1> IFCM>23/23/1 1 IER=8 F'IU=8. MM=M:*M NM=N*M DO 3 L=1/MM

218. TB=ABS<A(L>) IF<TB-PIU>3.r 3,22 PIU=TB I=L

219. CONTINUE TOL=EPS*PIU . LST=1

220. DO 17 К=ЬМ IF<PIU>23/23.'44 IF<IER>7.'5/75 IFCPIU-TODb.' б/ 76 IER-K-1

221. PIUI = 1.--A<I> J=<I-1VM I=I-J*M-K J=J+1-K1. DO 8 L=K.' NM/ M LL=L+I1. TB=PIUI*R<LL> R<LL>=R<L>

222. R(L>=TB IF<K-M>9.' 18.' 181. Э LEND=LST+M-K1.(J> 12.' 12.' 18 10 П=.ВД

223. DO 11 L=LST.' LEND ТВ=н < L) LL=L+II A<L>=A<LL>и a«::ll>=tb12 DO 13 L=L8T.'MM.'M LL=L+I

224. TB=PIUr*:ACLL> AÍ:LL>=A(L>13 A<L>=TBa<:lst>=.j1. PIU=0. LST=LST+1 , 1=0

225. DO 16 I I=LST.' LEND PIUI=-A<II> I8T=II+M J=.J+1

226. DO 15 L=IST.' MM.' M LL=L-J

227. A < L >=A <L)+PIUIжA < LL > TB=ABS<A<L>> IF<TB-PIU> 15.' 15.' 1414 PIU=TB I=L15 CONTINUE1. DO 16 L=K.'NM.'M1.=L+J

228. R < LL > =R < LL >+PIUI :*R < L >17 LST=LST+M18 IF < M-1 > 23."- 22.' 1919 IST=MM+M LST=M+11. DO 21 I =2.i-M1.=LST-I1.T=IST-LST1.IST-M1.ACL>+8.51. DO 21 J=I I.' NM.i- M1. TB=R(J)1.=J

229. DO 2@ K=IST/MM/M LL=LL+1 28 TB-TB-A <K>$R <LL > K=J+L RC,D=p <K>21 R<K)=TB22 RETURN23 IER=-1 RETURN END

230. KF=2550. DZi= S.158E-84 DZ2= 0.150E-04 ALFA= 1.55256 ОМ=1050.

231. Al = 8.18822Е-85 АЗ=-8.17888Е-83 А5= 8.24466Е-84 А7= 8.34185Е-84 АЗ=-8.18214Е-8о

232. В1= 8.43841Е-86 В3= 8.34271Е-86 В5= 8.55865Е-85 В7=-8.62988Е-83 В9=-8.21856Е-86 АК= 8.67968Е-85 AF= 8.18654Е-83 RR1= 18.981 RR2= 21.585 KF=2558. D21= 8.158E-04 HZ2= 8.158E-84 ALFA= 1.55256 0M= 638.

233. Al= 8.57386E-85 A3=-8.18487E-83 A5= 8.74579E-84 A7= 8.11824E-83 A9= 8.86982E-85

234. Bl= 8.16959E-85 B3= 8.11853E-85 B5= 8.21915E-84 B7=-8.86726E-83 B9= 8.25888E-85 AK= 8.21064E-04 AF= 0.20419E-03 RR1= 5.4533 RR2= 10.545 KF=2550. DZ1= 0.150E-04 HZ2= 0.500E-05 ALFA= 1.55256 0M= 630.

235. Al= 8.42619E-05 A3=-0.13671E-03 A5= 0.55388E-04 A7= 0.87812E-04 A9= 0.64541E-05

236. Bl= 0.12595E-05 B3= 0.88034E-06 B5= 8.16276E-84 B7=~8.64418E-83 B9= 8.18633E-05 AK= 0.15644E-04 AF= 0.15165E-03 RR1= 4.0500 RR2= 7.8318 KF=2550. DZ1= У.150E-04 DZ2= 0.500E-05 ALFA= 1.55256 0M=1858.

237. Al= 0.13979E-05 A3=-0.12691E-03 A5= 0.18171E-04 A7= 0.25388E-04 A9=-0.75856E-06

238. Bl= 0.14326E-06 B3= 0.11407E-06 B5= 0.18329E-05 B7=-8.28937E-83 B9=-0.72749E-07 AK= 0.16283E-05 AF= 0.48000E-04 RR1= 3.6255 RR2= 7.1535 KF=2558. D21= 8.158E-84 DZ2= 8.188E-84 ALFA= 3.12256 0M= 638.

239. Al= 8.13419E-85 A3=-8.43842E-84 A5= 8.17439E-84 A7= 8.27648E-84 A9= 0.20321E-05

240. Bl= 0.56449E-06 B3= 0.39455E-06 B5= 0.72945E-05 B7=-8.28867E-03 B9= 0.83508E-06 AK= 0.51230E-05 AF= 0.55076E-04 RR1= 1.8067 RR2= 3.4963 KF=2550. DZ1= 0.150E-04 DZ2= 0.100E-04 ALFA= 0.01744 0M= 630.

241. Al= 0.67019E-05 A3=-0.21497E-03 A5= 0.87098E-04 A7= 0.13808E-03 A9= 8.18149E-84

242. Bl= 8.28287E-86 B3= 8.14188E-86 B5= 8.26216E-85 E¡7=-8.18375E-83 B9= 8.38812E-86 AK= 8.23618E-84 AF= 8.19682E-83 RR1= 1.8789 RR2= 1.9753 KF=2558. DZ1= 8.158E-84 DZ2= 8.188E-84 ALFA= 8.81744 0M=1858.

243. Al= 8.21982E-85 A3=-8.19957E-83 A5= 8.28573E-84 A7= 0.39923E-04 A9=-0.11928E-05

244. Bl= 0.51488E-07 B3= 0.40997E-07 B5= 0.65872E-06 B7=-8.75245E-84 B9=-0.26146E-07 AK= 8.77435E-85 AF= 8.19444E-83 RR1= 1.4888 RR2= 2.8548 KF=2558. DZ1= 8.588E-85 DZ2= 8.158E-84 ALFA= 1.55256 0M=1858.

245. Al= 0.13979E-05 A3=-0.12691E-03 A5= 0.18171E-04 A7= 0.25388E-04 A9=-0.75856E-06

246. Bl= 8.31966E-86 B3= 8.25452E-86 B5= 8.48895E-85 B7=-8.46714E-83 B9=-8.16232E-86 AK= 8.50472E-05 AF= 0.13854E-03 RR1= 8.0957 RR2= 15.971 KF=2550. DZ1= 0.500E-05 DZ2= 0.150E-04 ALFA= 1.55256 0M= 630.

247. Al= 8.42619E-85 A3=-ñ.13671E-03 A5= 0.55388E-84 A7= 0.87812E-04 A9= 8.64541E-85

248. Bl= 0.12595E-05 B3= 0.88034E-06 B5= 0.16276E-04 B7=-8.64410E-03 B9= 0.18633E-05 AK= 0.15644E-04 AF= 0.15165E-03 RR1= 4.0500 RR2= 7.8318 KF=2550. DZ1= 0.500E-05 DZ2= 0.500E-05 ALFA= 1.55256 OM=1050.

249. Al= 0.62740E-06 A3=-0.56960E-04 A5= 0.81555E-05 A7= 0.11395E-04 A9=-0.34046E-06

250. Bl= 0.14347E-06 B3= 0.11424E-06 B5= 0.18355E-05 B7=-8.28967E-83 B9=-8.72853E-07 AK= 0.22653E-05 AF= 0.62180E-04 RR1= 3.6335 RR2= 7.1682 KF=2550. DZ1= 0.500E-05 DZ2= 0.500E-05 ALFA= 1.55256 0M= 630.

251. Al= 0.19129E-05 A3=-ñ.61357E-04 A5= 8.24868E-84 A7= 8.39412E-84 A9= 8.28967E-85

252. Bl= ñ.!=¡6530E-86 B3= 8.39512E-86 B5= 8.73049E-85 B7=-8.28989Е-83 B9= 0.83628E-06 AK= 0.70213E-05 AF= 0.68062E-04 RR1= 1.8178 RR2= 3.5151 KF=2550. DZ1= 0.500E-05 DZ2= 0.100E-04 ALFA= 3.12256 OM=Í050.

253. Al= 0.43981E-06 A3=-0.39929E-04 A5= 0.57170E-05 A7= 0.79879E-05 A9=-0.23866E-06

254. В1= 0.14316E-06 B3= 0.11399E-06 B5= 0.18315E-05 B7=-8.28921E-83 B9=-8.72696E-87 AK= 0.16271E-05 AF= 0.47965E-04 RR1= 3.6229 RR2= 7.1483

255. KF=2550. D21= 0.500E-05 DZ2= 8.180E-84 ALFA= 3.12256 0M= 630.

256. Al = 8. 13409E- •05 A3= -8.4381IE -84 A5: = 0. 17427E-84 A7= • 8. 27628E- -84 A9= 8. 28386E- -05

257. Bl = 0. 564Й8Е- •06 B3= 8.39427E -86 B5= = 0. 72892E-85 B7= :-8. 28847E- -83 B9= 8. 83447E- -86

258. AK= 0. 51193E- •05 AF= 8.55836E-I 84 RR1 = 1. 8854 RR2= 3. 4938

259. KF= =2550. DZ1 = i 0 .588E-85 DZ2= 8. 188b -04 ALFA= 0.01 744 0M= 1050.

260. Al = 0. 13189E- 85 A3= -8.11974E' -ЙЗ A5= = 0. 17144E-04 A7= ' 8. 23954E- -84 A9=- -8. 71578E- -ЙЙ

261. Bl = 0. 30893E- ■07 B3= 8.24599E -87 B5: = 0. 39524E-06 B7= ;-8. 45147E- -84 B9—■ -8. 15687E- -07

262. AK= 0. 46461E- •05 AF= 8.11666E-I 03 RR1 = 0. S 18851 RR2= 1. 7124

263. KF: =2550. DZ1 = ¡ 0 . 588 E-85 DZ2= 8. 100E- -04 ALFA= 0.01 744 0M= : 63Й.

264. Al = 0. 40212E- •05 A3= -8.12898E' -ЙЗ A5= = 0. 52259E-04 A7= î 8. 82851E- -84 A9= 8. 68894E- ■05

265. Bl = 0. Г2173Е- •06 B3= 8.85888E- -87 E:5: = 0. 15730E-05 B7= -0. 62249E- •84 B9= 8. 18888E- •06

266. AK= 0. 14166E- •04 AF= 8.11889E-1 83 RR1 = 0. 6 4736 RR2= 1. 1852

267. KF= =2550. DZ1 8 . 188 E-84 DZ2- 8. 100E- "04 ALFA= 1.55 256 0M= ■ 840.

268. Al = 8. 28772E- •85 A3= -8. 11649E' -83 A5: = 8. 26960E-04 A7= : 0. 37738E- •84 A9=- -8. 47448E- ■05

269. Bl = 0. 58797E- •86 B3= 8.37825E -86 B5: = 8. 65283E-05 B7= -0. 45887E- -83 B9=- -8. 11137E- ■05

270. AK= 0. 75198E- ■05 AF= 8.12728E-I 93 RR1 = 5. 8964 RR2= 9. 9921

271. KF: =2550. DZ1 = i 8 . 888 DZ2= 8. 800E- -84 ALFA= 8.81 744 0M= 8.1.F'R: = 4 1.PR: = 4

272. Al = 0. 88888 A3= 8.88888 A5: = 8. 88888 Hi' = • 0. 88800 A9= 8. 88888

273. Bl = 8. 88888 B3= 8.88888 B5: = 8. 88888 B7= ^ 0. £10088 B9= 8. 88888

274. AK= 0. 88888 AF= й.88888 RR1 = 8.88888 RR2= 0.00000

275. KF: =2550. DZ1 = 0 . 888 DZ2- 8. 800E- -84 ALFA= 8.81 744 0M= 8.1.PR' = 4 1.PR: = 4

276. Al = 0. 88888 A3= 8.88888 A5: = 8. 88888 A7= : 8. 88888 A9= 8. 88888

277. Bl = 0. 88888 B3= 8.00000 B5: = 8. 88888 B7= : 8. 00000 B9= 8. 88888

278. AK= 8. 88888 AF= 8. £10880 RR1 = 8.88088 RR2= 8.88888

279. KF: i-' '58. DZ1 = 1 0 .158E-84 DZ2= 8. 100E- -84 ALFA= 3.12 256 0M= : 840.

280. Al = 0. 28188E- •06 A3= -0.39150E -84 A5= = 8. 26896E-85 A7= ■ 8. 34775E- ■85 A9=- -8. 45914E- "06

281. Bl = 0. 61378E- •07 B3= 8.49626E' -87 B5: = 8. 78868E-86 B7= -8. 19349E- ■83 B9=- -8. 13448E- "06

282. AK= 0. 73911E- ■06 AF= 8.46647E-! 84 RR1 = 1438 RR2= 4. 2526

283. KF: = 750. DZ1 0 .158E-84 DZ2= 8. 180E- -84 ALFA= 8.81 744 0M= ^ 848.

284. Al = 0. 18843E- •05 A3= -й.19553E -83 A5: = 8. 13833E-84 A7= ^ 8. 17368E- ■84 A9=- -8. 22931E- -й5

285. Bl = 0. 22859E- ■07 B3= 8.17835E' -87 B5: = 8. 28345E-86 B7= -8. 69541E- ■84 B9=- -8. 48330E- -07

286. AK= 8. 35336E- ■05 AF= 8.19341E-I 83 RR1 = 8.7 '9665 RR2= 1. 5699

287. KF: = 750. DZ1 0 .588E-85 DZ2= 8. 100E- -84 ALFA= 3.12 25F. 0M= ■ 848.

288. Al = 8. 28893E- 06 A3= -8.39122E' -84 A5= = 8. 26877E-85 A7= • 8. 34758E- •85 A9=- -8. 45881E- -06

289. Bl = 0. 61334E- ■07 B3= 8.49598E' -87 B5: = 8. 78811E-86 B7= •-8. 19335E- •83 B9=- -8. 13438E- "06

290. AK= 0. 73858E- •06 AF= 8.46614E-! 84 RR1 = 2. 1415 RR2= 4. 2496

291. KF: = 750. DZ1 0 .588E-85 DZ2= 0. 100E- -84 ALFA= 8.81 744 on= : 848.

292. Al = 0. 60255E- •0 f, A3= -8.11732E -йЗ A5: = 8. 78288E-85 A7= ^ 8. 10421E- ■84 A9=- -8. 13759E- -05

293. Bl = 0. 13235E- •87 B3= 8.1878IE -87 B5: = 8. 17887E-86 B7= '~0. 41725E- -84 B9=- -8. 28998E- -07

294. AK= 0. 21201E- •05 AF= 8.11685E-! 03 RR1 = 8. A ■7799 RR2= 0.94194

295. KF: — I-' '58. DZ1 0 !. 100E-04 DZ2= 0. 150E- -84 ALFA= 3.12 256 on= : 848.

296. Al = 0. 28108E- 'fifi A3= -8.39158E -04 A5: = 8. 26896E-85 Л —1 Hf = ^ 0. 34775E- -85 A9=- -8. 45914E- "06

297. Bl = 0. 61378E- ■07 B3= 8.49626E -07 B5: = 8. 78868E-86 B7= •-0. 19349E- -83 B9=' -8. 13448E- "06

298. AK= 0. 73911E- •06 AF= 8.46647E-! 04 RR1 = iL ■ 1438 RR2- 4. 2526

299. KF= = 758. DZ1 = 8 i. 100E-04 DZ2= 0. 150E' -84 ALFA= 8.81 744 0M= = 848.

300. Al = 0. 10043E- •05 A3= -8. Г9553Е -03 A5; = 8. 13833E-84 A7= : 0. 17368E- -84 A9=' -8. 22931E- "05

301. Bl = 0. 22059E- •07 B3= 8.17835E -07 B5: = 8. 28345E-86 B7= ■-0. 69541E- -84 B9=- -8. 48330h- -07

302. AK= 0. 35336E- •05 AF= 8.19341E- 03 RR1 = 8.79665 RR2= 1. 5699

303. KF: — г* '50. DZ1 8 .188E-84 DZ2= 0. 5Й0Е- -85 ALFA= 3.12 256 0M= : 848.

304. Al = 8. 28093E- •06 A3= -8.39122b -Й4 A5: = 8. 26877E-85 A7= ■ £1. 34758E- -85 A9=' -8. 45881E- -06

305. Bl = 0. 61334E- ■07 B3= : 8.49598E -07 B5; = 8. 78811E-86 B7= 0. 19335E- -83 B9=' -8. 13438E- "06

306. АК= 8.73858Е-06 AF= 8.46614Е-84 RR1= 2.1415 RR2= 4.2496 KF= 75@. DZ1= 8.188E-04 DZ2= 0.588E-85 ALFA= 8.81744 0M= 848.

307. Al= 8, 6Й255Е-86 Bl= 8.13235E-87

308. A3=-8.11732E-03 B3= 8.18781E-871. A5= B5=8.78288E-85 8.17887E-861. A7= 8.18421E-84

309. AK= 8.21281E-85 AF= 8.11685E-83 RR1= 8.47799

310. B7=-@.41725E-84 RR2= 8.94194

311. KF= 758. DZ1 = 8.188E-84 DZ2= 8.158E-84 ALFA= 1.55256 0M=1858. A1= 8.45191E-86 A3=-8.14133E-83 A5= 8.58748E-85 A7= 8.79635E-85 Bl= 8.95618E-87 B3= 8.79835E-87 B5= 8.12238E-85 B7=-8.48254E-83 AK= 8.16259E-85 AF= 8.15461E-83 RR1= 8.3523 RR2= 16.682

312. A9=-8.13759E-85 B9=-8.28998E-87

313. Й9=-0.24519E-86 B9=-0.48519E-87

314. KF= Al = Bl = AK= KF= Al= Bl = AK= KF= Al = Bl = AK= KF= Al = Bl = AK= K'F-Al = Bl = AK= KF= Al = Bl = AK= KF= Al = Bl = AK= KF= Al = Bl = AK= KF= Al = Bl = AK= KF= Al = Bl = AK= KF= Al = Bl = AK= KF= Al =

315. B3= 8.49457E-87 B5= 85 AF= 8.95778E-84 RR1= = 8.188E-84 DZ2- 8.588E85 8. 8.5. 85 8.

316. B3= 8.14467E-06 B5= 0. 05 AF= 8.97478E-84 RR1= 2. = 8.158E-84 DZ2= 8.188E 85 A3=-8.42834E-8486 A3=-8.89272E-84 A5=1. A5=

317. B3= 8.36981 E-86 B5= 85 AF= 0.54484E-04 RR1= = 8.150E-84 DZ2= 8. 199E05 A3=-0.21393E-03 A5:06 B3= 0.13262E-06 B5:

318. AF= 0.19696E-03 RR1= = 8.588E-05 DZ2= 0.100E05 A3=-0.12836E-031. A5=-04 : 0. = 8.3. -84 = 8. = 8. 1. -84 = 8.86 B3= 8.79572E-87 B5= 8.1. 84

319. AF= 8.11818E-83 RR1= = 8.588E-85 DZ2= 8.188E85 A3=-8.42883E-84 A5

320. B3= 8.36874E-86 B5 85 hF= 8.54444E-84 RR1= = 0.100E-04 DZ2- 0.150E0.05 A3=-0.42834E-04 A5=04 0.

321. B3= 0.36901E-06 B5= 0. 05 AF= 0.54484E-04 RR1= 3. 0.100E-04 DZ2- 0.150E-0405 A3=-0.21393E-03 A5= 0.

322. B3= 0.13262E-06 B5= 8. @4 hiF= 8.19696E-83 RR1 = 8.188E-84 DZ2= 8.588E85 A;

323. A9= 8.19369E-85 B9= 0.42144E-06

324. A9=-0.15189E-06 B9=-0.30057E-07

325. A9= 0.11999E-05 B9= 0.26108E-06

326. A9=-0.30348E-85 B9=-8.12821E-85

327. A9=-8.15153E-84 В9=-й.43283E-86

328. A9=-8.98917E-85 B9=~8.25922E-86

329. A9=-8.38318E-85 B9=-8.12812E-85

330. A9=-8.38348E-85 B9=-8.12821E-85

331. A9=-8.15153E-04 B9=-0.43203E-06

332. A9=-0.90917E-85 B9=-0.25922E-86

333. A9=~8.30318E-85 B9=-8.12812E-851. A9=-8.15414E-85

334. Bi= 8.71884Е-86 В3= 8.53655Е-86 В5= 8.98956Е-85 В7=-8.59S97E-83 В9=-8.36121Е-8ь АК= 8.18384Е-84 AF= 8.16318Е-83 RRi= 18.485 RR2= 28.367 KF=4358. HZi= 8.188E-84 HZ2= 8.158E-84 ñLFft= 1.55256 0M= 638.

335. Al= 8.92915E-85 A3=-8.17243E-83 A5= 8.12877E-83 A7= 8.25157E-83 A9= 8.14881E-84

336. Bl= 8.45781E-85 B3= 8.25228E-85 B5= 8.59868E-84 B7=-8.12942E-82 B9= 8.67666Е-Й5 AK= 8.36437E-84 AF= 8.22897E-83 RR1= 8.2134 RR2= 15.533 KF=4358. DZ1= 8.188E-84 DZ2- 8.588E-85 ALFA= 1.55256 0M= 638.

337. A1= 8.57568E-85 A3=-8.18682E-83 A5= 8.74818E-84 A7= 8.15584E-83 A9= 8.87227E-85

338. Bl= 8.28311E-85 B3= 8.15624E-85 B5= 8.36592E-84 B7=-8.88174E-83 B9= 8.41918E-85 AK= 8.22572E-84 AF= 8.13689E-83 RR1= 5.8881 RR2= 9.6222 KF=4358. DZ1= 8.188E-84 DZ2= 8.588E-85 ALFA= 1.55256 0M=1858.

339. Al= 8.17593E-85 A3=-8.92467E-84 A5= 8.22878E-84 A7= 8.33155E-84 A9=-8.95498E-86

340. Bl= 8.44836E-86 E¡3= 8.33239E-86 B5= 8.56346E-85 B7=-8.37185E-83 B9=-8.22377E-86 AK= 8.63831E-85 AF= 8.18189E-83 RR1= 6.4468 RR2= 12.617