# =====================================================================
#   СКРИПТ-ЗРАЗОК
#
#   Тут показано, як змусити R порахувати покоління одне за одним.
#   Прочитайте коментарі, запустіть, подивіться на результат.
#
#   Як запускати: поставте курсор на рядок і натисніть Ctrl+Enter.
#   R виконає цей рядок і перейде на наступний. Тисніть далі.
# =====================================================================


# ---------------------------------------------------------------------
#   МОДЕЛЬ 1.  Популяція, що щороку подвоюється
#
#        N     =  N  * 2
#         i+1      i
# ---------------------------------------------------------------------

n <- 10
# Стрілка <- читається "поклади в". У коробочку з іменем n "кладемо" 10.
# Ми створили об'єкт. Тут будемо позначати кількість поколінь у розрахунку.

N <- numeric(n)
# Заготовка: низка з n чисел, поки що всі нулі.
# Така низка називається вектор. Її місця пронумеровані: 1, 2, 3, ...

N[1] <- 5
# Квадратні дужки означають "яке саме місце".
# У перше місце об'єкту N "кладемо" початкову чисельність — 5 особин.

for (i in 1:(n - 1)) {
  N[i + 1] <- N[i] * 2
}
# for означає "повтори", це команда-цикл. Літера i по черзі стає 1, 2, 3, ... , n-1.
# Усередині дужок { } записано, що саме повторювати.
# Рядок, що виконується в циклі, читається так: у місце i+1 "поклади" те, що в місці i, помножене на 2.

N
# Просто напишіть ім'я об'єкта — і RStudio у консолі (знизу зліва) покаже, що в ньому.

plot(N, type = "b")
# Намалювати. type = "b" означає both: і точки, і лінії між ними.


# ---------------------------------------------------------------------
#   МОДЕЛЬ 2.  Популяція, до якої щороку додають 10 особин
#
#        N     =  N  + 10
#         i+1      i
#
#   Порівняйте цей блок із попереднім.
#   Змінився РІВНО ОДИН рядок — той, у якому записана формула.
#   Решта команд ті самі.
# ---------------------------------------------------------------------

n <- 10
N <- numeric(n)
N[1] <- 5

for (i in 1:(n - 1)) {
  N[i + 1] <- N[i] + 10
}

N
plot(N, type = "b")


# =====================================================================
#   ПІДСУМОК: у цьому скрипті всього п'ять дій
#
#     n <- 10              скільки кроків рахуємо
#     N <- numeric(n)      заготовка під результат
#     N[1] <- 5            з чого починаємо
#     for (...) { ... }    сама формула, повторена n-1 раз
#     plot(N, type = "b")  подивитися, що вийшло
#
#   Щоб порахувати будь-яку іншу модель, треба змінити тільки формулу.
# =====================================================================