# =====================================================================
#   ДОМАШНЄ ЗАВДАННЯ
#
#   У цьому скрипті бракує одного рядка. Допишіть його — за зразком
#   із файлу Sample.R. Усі потрібні команди там уже є; змінюється
#   тільки формула.
#
#   Формула Мея:
#
#        x     =  lambda * x  * (1 - x )
#         i+1                i         i
#
# =====================================================================


# ---------------------------------------------------------------------
#   ЧАСТИНА 1.  Допишіть рядок
# ---------------------------------------------------------------------

lambda <- 2.8          # коефіцієнт розмноження
n      <- 40           # скільки поколінь рахуємо

x <- numeric(n)
x[1] <- 0.15           # початкова чисельність

for (i in 1:(n - 1)) {

  # <<<<<<<<<<  СЮДИ ВПИШІТЬ ОДИН РЯДОК  >>>>>>>>>>
  # Зліва має бути x[i + 1], потім стрілка <-,
  # а праворуч — формула Мея, записана мовою R.
  #
  # Пам'ятайте: знак множення в R не можна пропускати.
  # Там, де в математиці пишуть lambda·x, у R пишуть lambda * x

}

plot(x, type = "b", main = paste("lambda =", lambda))
# paste склеює слова й числа в один підпис,
# щоб на малюнку було видно, з якою lambda він побудований


#   ПЕРЕВІРКА. Якщо все правильно, то при lambda = 2.8
#   чисельність вийде на сталий рівень. Наберіть:

x[40]

#   Має бути приблизно 0.6429
#   Якщо у вас нуль — рядок ще не вписано.
#   Якщо інше число — шукайте помилку у формулі.


# ---------------------------------------------------------------------
#   ЧАСТИНА 2.  Знайдіть три режими
#
#   Тепер міняйте ТІЛЬКИ значення lambda в першому рядку
#   і щоразу запускайте скрипт заново.
#
#   Знайдіть такі значення lambda, за яких популяція поводиться так:
#
#     A) чисельність виходить на стале значення й більше не міняється
#     B) чисельність стрибає між двома значеннями, з року в рік
#     C) чисельність не повторюється ніколи
#
#   Пробуйте lambda від 2 до 4.
# ---------------------------------------------------------------------

#   ВІДПОВІДІ:
#
#   A) lambda =
#   B) lambda =
#   C) lambda =