ДИНАМИКА УРОЖАЙНОСТИ УЛУЧШЕННЫХ СЕНОКОСОВ

Введение
Для функционального описания поведения территории нами вводятся новые понятия - активность и интенсивность растительного покрова. Причем территория понимается как простейшее геодезическое изображение ландшафта. А сам ландшафт, в свою очередь, является первым компонентом динамической геотриады «ландшафт + население + хозяйство».
Активность учитывается по доле площади растительного покрова (леса и древесно-кустарниковая растительность, луга и пастбища, особо охраняемые территории и болота) и этот экологический параметр позволяет хаpaктеризовать фактически образовавшиеся отклонения от территориального экологического равновесия [1] на конкретной территории.
Однако пока отсутствуют данные земельных кадастров по конкретным землепользователям. Поэтому приведены значения параметров растительного покрова по административным районам и городам субъекта федерации. В частности рассмотрены районы и города Республики Марий Эл (РМЭ) по состоянию распределения земель на 01.01.07 г.
В наиболее общем случае интенсивность проявляется как активность во времени. Физически интенсивность - это скорость изменений. А активность - это сами изменения в природной, природно‒техногенной или технической среде (по площади, урожайности растений, продуктивности почвы и пр.) в некотором срезе времени.
Ритмика жизнедеятельности растений четко предопределена циклами вращения Земли вокруг Солнца, поэтому единицей измерения любого показателя (параметр - это показатель, хаpaктеризующий исследуемую систему) растительного покрова и его компонент станет год. Активность и интенсивность поведения человека по отношению к процессу уничтожения ландшафта быстротечны. Например, планировалось атомными взрывами построить русло крупного канала для переброски воды северных рек на юг.
Растительный покров - это территориально распластанное живое существо, обладающее многими свойствами поведения живого вещества. Индикатором экологического равновесия на данной территории становится его активность. Она в простейшем случае исчисляется по доле занятой площади, причем в ходе конкуренции за плодородные почвы, с человеком. А интенсивность поведения (скорость реакции на внешние раздражения, в том числе и на антропогенные воздействия) растительного покрова за более чем 350‒450 млн. лет эволюции наземных растений пока не изменилась. За скоростью эволюции человека растения явно не успевают. Но они это возместят высокой своей активностью роста и развития сорняков, а также ликвидацией нужных человеку видов культур болезнями.
Травяной покров. Особое место в растительном покрове занимает травяной покров, и он более живуч, устойчив и стабилен в своем поведении по сравнению с древесными пологами. По данным американских ученых, трава появилась около 100 млн. лет назад. А деревья возникли 350‒450 млн. лет назад, то есть трава появилась как результат отклика древесных растений на поведение наземных животных. Поэтому трава биологически устойчивее деревьев, кустарников и кустарничков. При снижении плодородия почвы и ухудшении климата деревья исчезают, а их место обитания занимают степные травы.
Интенсивность воспроизводства травяного покрова. Годичная цикличность определяет урожайность травы в естественных условиях роста и развития, какими являются естественные луга, на которых заготовляется сено. Улучшенные луга уже меняют свое поведение, повышая урожайность на позитивное воздействие человека культурными и техническими мерами и внесением на такие луга удобрений.
Из статистического сборника [6, с. 181] были взяты сводные данные по сенокосным лугам Республики Марий Эл (РМЭ), приведенные в таблице 1.
Таблица 1
Урожайность сена в сельских хозяйствах на территории РМЭ, ц/га
|
Годы учета* |
в хозяйствах всех категорий |
в крупных и средних хозяйствах |
|||||
|
Время t, лет |
Естественные и улучшен. сенокосы |
Однолетние травы |
Много-летние травы |
Время t, лет |
Естественные сенокосы |
Улучшенные сенокосы |
|
|
1970 1975 1980 1985 1990 1995 1996 1997 1998 1999 2000 2001 2002 2003 |
0 5 10 15 20 25 26 27 28 29 30 31 32 33 |
14.5 8.0 12.0 17.1 19.0 3.0 12.2 16.0 11.3 9.9 13.6 14.6 9.1 12.5 |
12.4 7.9 20.5 29.9 26.9 16.1 30.4 18.8 21.4 11.1 33.7 20.6 17.7 21.1 |
19.9 10.1 18.6 35.1 33.6 16.8 18.8 21.3 16.1 15.5 18.2 19.0 12.1 18.6 |
- - 0 5 10 15 16 17 18 19 20 21 22 23 |
- - 11.0 14.7 16.7 12.1 10.0 13.0 10.2 9.5 12.0 14.1 8.7 12.8 |
- - 16.3 21.8 26.2 15.4 15.2 21.1 14.6 11.3 18.8 15.5 10.4 11.2 |
Примечание: * Недостатком динамического ряда является отсутствие данных по всем годам в промежутках между пятилетками, что не позволяет точно выявить циклы с полупериодом менее пяти лет.
Вначале рассмотрим статистические закономерности динамики урожайности травы после естественной сушки в сено на сенокосных лугах сельских хозяйств всех категорий. При этом будем выполнять все этапы статистического моделирования:
‒ эвристическая идентификация, когда выясняется логика каждого динамического ряда исходных количественных данных;
‒ структурная идентификация, когда по эвристической модели понимание процессов изменения значений изучаемого показателя выбирается основная (естественная составляющая) закономерность, то есть тренд изменения во времени исследуемого показателя;
‒ параметрическая идентификация, когда в решающей математической среде типа CurveExpert‒1.3 подбираются значения параметров статистической модели исходной конструкции, полученной на предыдущем этапе структурной идентификации.
После выполнения всех процедур и этапов идентификации биотехнического закона, предложенного проф. П.М. Мазуркиным [2‒5], сформируется готовая статистическая модель. Урожайность сена с естественных и улучшенных сенокосов со всей площади по РМЭ определяется рядом значений по точечному графику на рис. 1.
Рис. 1. Изменение урожайности сена естественных и улучшенных сенокосов
Республики Марий Эл, ц/га
Из роя точек резко выделяется одна (самая нижняя) точка.
Полнота динамического ряда. Пусть исходные данные [6, с.181] количественно достоверны и статистики адекватно описали табличными моделями типа табл. 1 картину динамики исследуемого явления или процесса. Тогда можно поставить задачу оценки полноты статистического ряда динамики значений изучаемого показателя.
Критерий полноты значений показателя. Представим динамический ряд значений урожайности сена с естественных и улучшенных сенокосов со всей площади по РМЭ в виде последовательности, показанной в таблице 2.
Весь этот ряд динамики можно разделить на две части:
а) динамика по пятилетним циклам, причем четко по хронологическому времени 1970, 1975, 1980 и т.д.;
б) динамика по годичным циклам, что соответствует действительности, то есть ритму вращения Земли вокруг Солнца.
Таблица 2
Полный динамический ряд урожайности сена, ц/га
|
Время t, лет |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
|
|
Урож. q, ц/га |
14.5 |
- |
- |
- |
- |
8.0 |
- |
- |
- |
- |
|
|
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
21 |
|
12.0 |
- |
- |
- |
- |
17.1 |
- |
- |
- |
- |
19.0 |
- |
|
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
33 |
|
- |
- |
- |
3.0 |
12.2 |
16.0 |
11.3 |
9.9 |
13.6 |
14.6 |
9.1 |
12.5 |
Введем новый показатель оценки добротности статистических данных динамического ряда - это коэффициент полноты ряда. Этот коэффициент k будет вычисляться как отношение количества m имеющихся значений показателя в ряду к общему количеству возможных члeнов n динамического ряда, то есть по формуле k=m/n .
Из данных табл. 2 получаем m=14 и n=34. Поэтому коэффициент полноты ряда будет k=14/34=0,412. Динамический ряд оказался сильно деформированным. Специальными уловками группировки исходных данных вовсю занимались советские экономисты, чтобы доказывать недоказуемое. Нужно статистическим службам запретить публикации неполных динамических рядов, так как это почти всегда приводит к ложной идентификации.
Картина распределения точек, показанная роем на рисунке 1, оказалась ложной из-за отсутствия 20 точек из требуемых 34 значений изучаемого показателя. При этом не ясно, какими по расположению относительно принятой системы координат были точки в промежутках между пятилетками. Поэтому может оказаться, что полный статистический ряд даст совершенно иную модель динамики (если урожайности 1970, 1975, 1980, 1985, 1990, 1995 годов окажутся не максимальными, а иными). Вот почему надо обращать серьезное внимание на полноту статистического ряда динамики значений изучаемого показателя.
Границы доверительного интервала статистического ряда. Применим известный прием выделения верхней и нижней границ доверительного интервала. Это обезопасит от ложной идентификации и неправильной тpaктовки найденной формулы, так как верхняя и нижняя границы показывают только возможное поле распределения значений показателя вне зависимости от полноты статистического ряда.
Перепишем первую часть показателей таблицы 1 с учетом границ доверительного интервала и результаты приведем в таблице 3.
Верхняя граница урожайности сена естественных и улучшенных сенокосов определится уравнением (рис. 2) вида
qmax=14,5092exp(-0,0091099t)+4,7361×10-8t9,23434exp(-0,44353t). (1)
Максимальная относительная погрешность уравнения (1) равна 4,59 %. Коэффициент корреляции по рис. 2 равен 0,9861.
Нижняя граница урожайности сена (рис. 3) определится статистическим уравнением
qmin = 7.64985exp(-0,0052882t)+0,000403119 t 5,45669 exp(-0,31879t). (2)
Максимальная относительная погрешность уравнения (2) равна 3,84 %, а коэффициент корреляции по рисунку 3 равен 1,0000, то есть формула (2) становится однозначно функциональной математической зависимостью.
Таблица 3
Урожайность сена в сельских хозяйствах всех категорий на территории РМЭ, ц/га
|
Годы учета |
Время t, лет |
Естественные и улучшенные сенокосы |
Однолетние травы |
Многолетние травы |
||||||
|
ряд |
max |
min |
ряд |
max |
min |
ряд |
max |
min |
||
|
1970 1975 1980 1985 1990 1995 1996 1997 1998 1999 2000 2001 2002 2003 |
0 5 10 15 20 25 26 27 28 29 30 31 32 33 |
14.5 8.0 12.0 17.1 19.0 3.0* 12.2 16.0 11.3 9.9 13.6 14.6 9.1 12.5 |
14.6
17.1 19.0
16.0
14.6
12.5 |
8.0 12.0
12.2
9.9
9.1 |
12.4 7.9 20.5 29.9 26.9 16.1 30.4 18.8 21.4 11.1 33.7 20.6 17.7 21.1 |
12.4
20.5 29.9
30.4
33.4
21.1 |
7.9
16.1
18.8
11.1
17.7 |
19.9 10.1 18.6 35.1 33.6 16.8 18.8 21.3 16.1 15.5 18.2 19.0 12.1 18.6 |
19.9
35.1 33.6
21.3
19.0
18.6 |
10.1
16.8 18.8
16.1 15.5
12.1 |
Примечание: * Такие резко отклоняющиеся точки исключаются из статистических рядов.
Рис. 2. Верхняя граница доверительного интервала урожайности сена естественных
и улучшенных сенокосов РМЭ, ц/га
Закономерности динамики урожайности сена. Конструкция уравнений для описания верхней и нижней границы доверительного интервал одинакова. Она содержит две составляющие:
а) закон гибели (спада) в упрощенной конструктивной форме, когда третий параметр статистической модели будет равен единице;
б) биотехнический закон проф. П.М. Мазуркина [2‒5] в упрощенной форме (когда четвертый параметр модели равен единице).
Когда известна общая конструкция статистической закономерности, то можем её применить ко всему динамическому ряду новых исходных статистических данных.
Рис. 3. Нижняя граница доверительного интервала урожайности сена естественных и улучшенных сенокосов РМЭ, ц/га
В результате идентификации был получен тренд (рис. 4) из двух составляющих в виде устойчивых законов, в котором знак экспоненциального закона изменился с отрицательного на положительный, и он стал показывать рост значений урожайности сена. При этом 1995 год не был исключен.
q=10,9017exp(0,00083530t)+1,1128×10-10 t 14,0597 exp(-0,88382t), (4)
После исключения резко отклоняющейся от других значений показателя точки, отмеченной в данных табл. 3 звёздочкой, была получена формула (рис. 5) статистической закономерности вида
q=11,4401exp(-0,0010562t)+4,9549×10-12 t 14,74586 exp(-0,80755t). (5)
Коэффициент корреляции увеличился до 0,7266, а закономерность стала одинаковой с границами доверительного интервала.
Рис. 4. Тренд урожайности сена естественных и улучшенных сенокосов РМЭ, ц/га
Рис. 5. Тренд урожайности сена естественных и улучшенных сенокосов РМЭ после исключения резко отклоняющейся точки, ц/га
При этом попытка наращивания первой составляющей до закона гибели в полной форме не дает результата. Остатки приведены на рис. 6.
Чередование знаков остатков ε показывает, что можно дополнить составляющей в виде закона колебательного движения.
Структурно-параметрическая идентификация выполняется относительно вейвлет‒функции проф. П.М. Мазуркина [2‒5] вида
q3 = Acos(πt/p0.5±a10), (6)
где A - амплитуда (половина) колебательного движения значений изучаемого показателя во времени по биотехническому закону
(6а)
p0.5 - половина периода колебательного возмущения показателя по биотехническому закону (с учетом предыстории в виде постоянного члeна)
, (6б)
a1...a10- параметры исходной конструкции статистической модели (6), являющейся обобщенным законом колебательного изменения значений изучаемого показателя.
Волновую закономерность нельзя искать по статистическим выборкам с пропусками данных. Поэтому принимается интервал времени с 1995 по 2003 гг. (причем результат 1995 года исключается).
После нескольких сеансов поиска значений параметров модели (6) в математической среде CurveExpert‒1.3 была получена конкретная закономерность в виде сложного уравнения
q=5,64537exp(0,022908t)+6,0178×10-13 t 14,8210 exp(-0,74130t)-A1cos(πt/p1+3,29024)+
+A2cos(πt/p2+1,31367), (7)
A1 = 2,09541 t 0,15870 exp(-0,0037692t), p1 = 2,79259 - 0,019856t,
A2 = 4,1925×10-9 t 8,82490 exp(-0,38095t), p2 = 1,12563 + 0,00014490t.
Рис. 7. Динамика урожайности сена естественных и улучшенных сенокосов РМЭ за период 1970‒2003 гг., ц/га
В статистической модели (7) из четырех составляющих две последние являются колебательными возмущениями. Из-за отрицательного знака первая волновая зависимость является кризисным возмущением сенокосных лугов на внешние раздражители (естественные и антропогенные), причем частота возмущения растет, а вторая волна с начальным периодом в 2×1,12563 = 2,25 лет показывает позитивное успокоение травяного покрова со снижением частоты адаптивного возмущения.
Сокращение ряда до 1996‒2003 гг. сужает прогнозные возможности готовой статистической модели. Исходя из высокой точности, когда максимальная относительная погрешность равна всего 0,47 %, уравнение (7) дает возможность принять горизонт прогноза, равный основанию прогноза, то есть в восемь лет до
2003 + 8 = 2011 года.
А при ряде со всеми 34 точками получили бы возможность прогноза до
2003 + 34 = 2037 года.
При этом остатки располагаются так, как это показано на рис. 8.
Расположение остатков показывает, что существует еще одна волновая составляющая с примерным периодом в пять лет. Погрешность измерений урожайности сена недостаточно высока, чтобы искать вейвлет‒функцию с амплитудой менее 0,08 ц/га (запись в справочнике [6, с.181] выполнена с округлением до одной десятой, поэтому точность измерений составляет ±0,05 ц/га).
Рис. 8. Остатки готовой статистической модели тренда (4.11) урожайности сена естественных и улучшенных сенокосов РМЭ, ц/га
Изменение шкалы времени. Если исключить пятилетние данные, а также резко отклоняющуюся точку 1995 г., то получим динамический ряд с новой шкалой абсцисс. Переход выполняется по формуле t1=t‒26.
Для новой шкалы получаем условие t1 = 0 для 1996 года. Такое сокращение времени приводит и к уменьшению количества составляющих модели.
После повторной параметрической идентификации модели (4.11) была получена формула (рис. 9) с тремя составляющими
q=12,3327exp(-0,00036540t) -A1cos(πt/p1+1,53433)+A2cos(πt/p2 -2,93294), (8)
A1 = 3,63948exp(0,0089534t), p1 = 1,91645-0,028705t,
A2 = 0,034478 t 8,51786 exp(-1,80140t), p2 = 1,16822 - 0,0086746t
Рис. 9. Динамика урожайности сена естественных и улучшенных сенокосов РМЭ за период 1996‒2003 гг., ц/га
Рис. 10. Остатки готовой статистической модели тренда (4.12) урожайности сена естественных и улучшенных сенокосов РМЭ за период 1996‒2003 гг., ц/га
По коэффициенту корреляции 1,000 и нулевым остаткам (рис. 10) модель (8) стала однозначной.
Даже по малым остаткам, находящихся в интервале пренебрежимо малых для пpaктики измерения массы сена чисел ǀɛǀ=7·10-12...7·10-10 на рис. 10 наблюдается новая волновая закономерность. Этот факт означает, что повышение точности измерения массы травы и сена с единицы площади позволит находить статистические закономерности динамики урожайности.
Список литературы
- Реймерс, Н.Ф. Природопользование: Словарь‒справочник / Н.Ф. Реймерс. - М.: Мысль, 1990. - 637 с.
- Мазуркин, П.М. Геоэкология: Закономерности современного естествознания: Научное изд. / П.М. Мазуркин. - Йошкар-Ола: МарГТУ, 2006. - 336 с.
- Мазуркин, П.М. Закономерности устойчивого развития / П.М. Мазуркин. ‒ Научное издание. - Йошкар-Ола: МарГТУ, 2002. - 302 с.
- Мазуркин, П.М. Рациональное природопользование: учебное пособие. В 3‒х ч. Ч. 1: Экологически ответственное землепользование / П.М. Мазуркин, С.Е. Анисимов, С.И. Михайлова; под ред. П.М. Мазуркина. - Йошкар-Ола: МарГТУ, 2006. - 176 с.
- Мазуркин, П.М. Математическое моделирование. Идентификация однофакторных статистических закономерностей: Учебное пособие / П.М. Мазуркин, А.С. Филонов. - Йошкар-Ола: МарГТУ, 2006. - 292 с.
- Республика Марий Эл: Статистический ежегодник / А.В. Целищев и др. - Часть 1. - Йошкар-Ола: ФСГС, 2004. - 342 с.
Статья подготовлена и опубликована при поддержке гранта 3.2.3/4603 МОН РФ
Статья в формате PDF
108 KB...
21 07 2026 15:46:15
Статья в формате PDF
113 KB...
20 07 2026 13:33:52
Статья в формате PDF
121 KB...
19 07 2026 3:54:32
Статья в формате PDF
101 KB...
18 07 2026 0:12:24
17 07 2026 2:49:17
В статье проанализирован опыт лечения больных острым аппендицитом за последние 10 лет. Из 1073 поступивших в приемное отделение, 229 больных отправлены в другие отделения, у 730 диагноз подтвержден и выполнена операция аппендэктомия. Гистологическое исследование отростков показало, что у 353 (48,4%) больных отросток был флегмонозный, у 87 (11,9%) – гангренозный, в том числе у 15 (2%) – гангренозно-перфоративный, у 290 (39,7%) – катаральный. Большой процент катаральных форм автор связывает с гипердиагностикой. 24 (3,2%) больных был диагностирован разлитой перитонит. В комплексном лечении больных наряду с антибактериальными средствами, последнее время широко стали применяться современные методики (дренирование брюшной полости силиконовыми трубками д 0,5-1,0 см, назогастральное дренирование, гемосорбция, УФО крови, химическая детоксикация гипохлоритом натрия). ппендикулярный инфильтрат был диагностирован у 14 (1,9%) больных. Тактика при этом осложнении была традиционной. У 35 (4,79%) больных развились послеоперационные осложнения: нагноение подкожно-жировой основы у 19 (2,66%), инфильтраты послеоперационного шва – у 9 (1,2%), гематомы подкожной клетчатки – у 7 (0,9%), в том числе у 7 (0,9%) больных с нагноением подкожно-жировой основы, развились дополнительно послеоперационные пневмонии. а эти годы серьезных полостных послеоперационных осложнений не отмечалось, также не было послеоперационной летальности. лучшение результатов лечения автор связывает с повышением профессионального роста врачей, продуманной взвешенной хирургической тактикой. Также имеет значение и возраст больных. У 88% он равнялся 1822 годам. При поступлении больные были физически крепкими и тренированными (военнослужащие), что позволило им значительно лучше справиться в послеоперационном периоде даже с перитонитом.
...
16 07 2026 12:42:33
15 07 2026 16:58:36
14 07 2026 19:15:32
Статья в формате PDF
255 KB...
13 07 2026 4:39:52
Статья в формате PDF
126 KB...
12 07 2026 23:35:34
Статья в формате PDF
311 KB...
11 07 2026 12:12:20
Статья в формате PDF
258 KB...
10 07 2026 6:13:12
Статья в формате PDF
260 KB...
09 07 2026 12:53:20
08 07 2026 3:53:42
Статья в формате PDF
206 KB...
07 07 2026 4:46:20
Статья в формате PDF
242 KB...
06 07 2026 22:29:50
Статья в формате PDF
267 KB...
05 07 2026 5:41:16
Статья в формате PDF
737 KB...
04 07 2026 3:25:58
Статья в формате PDF
109 KB...
03 07 2026 15:40:19
Статья в формате PDF
241 KB...
02 07 2026 13:17:57
01 07 2026 12:13:44
Статья в формате PDF
147 KB...
30 06 2026 0:12:51
Статья в формате PDF
127 KB...
29 06 2026 7:20:46
Статья в формате PDF
101 KB...
27 06 2026 19:17:31
Статья в формате PDF
114 KB...
26 06 2026 18:34:49
Статья в формате PDF
124 KB...
24 06 2026 14:36:13
Статья в формате PDF
108 KB...
23 06 2026 20:35:37
Статья в формате PDF
275 KB...
22 06 2026 10:47:45
В статье постовариоэктомический синдром рассматривается как предиктор метаболического синдрома у женщин различных возрастных групп. На результатах анализа разнообразного клинико-диагностического материала показано, что женщин с постовариоэктомическим синдромом в возрасте после 40 лет достоверно чаще наблюдаются метаболические нарушения.
...
21 06 2026 14:32:21
Статья в формате PDF
112 KB...
20 06 2026 6:23:34
Статья в формате PDF
137 KB...
19 06 2026 8:58:35
Статья в формате PDF
376 KB...
18 06 2026 23:52:40
Статья в формате PDF
104 KB...
17 06 2026 22:37:43
Статья в формате PDF
120 KB...
16 06 2026 23:19:40
Статья в формате PDF
105 KB...
15 06 2026 5:12:47
Статья в формате PDF
314 KB...
14 06 2026 2:18:57
Статья в формате PDF
105 KB...
13 06 2026 18:10:27
Современный этап развития мирового и отечественного языкознания хаpaктеризуется антропоцентрической направленностью лингвистических исследований. Антропоцентризм является одним из фундаментальных свойств человеческого языка, так как взаимосвязь и взаимообусловленность языка и человека очевидна и не может вызывать никаких сомнений. «Идею антропоцентричности языка в настоящее время можно считать общепризнанной: для многих языковых построений представление о человеке выступает в качестве естественной точки отсчета» [1, 5]. Антропоцентрический подход в изучении языка или антропоцентрическая парадигма предполагает анализ человека в языке и языка в человеке. В.А. Маслова пишет, что «…антропоцентрическая парадигма выводит на первое место человека, а язык считается конституирующий хаpaктеристикой человека, его важнейшей составляющей. Человеческий интеллект, как и сам человек, немыслим вне языка и языковой способности как способности к порождению и восприятию речи. Если бы язык не вторгался во все мыслительные процессы, если бы он не был способен создавать новые ментальные прострaнcтва, то человек не вышел бы за рамки непосредственно наблюдаемого. Текст, создаваемый человеком, отражает движении человеческой мысли, строит возможные миры, запечатлевая в себе динамику мысли и способы ее представления с помощью средств языка» [1, 8].
...
12 06 2026 8:19:11
Еще:
Поддержать себя -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 ::