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

rm (list = ls ())          # очищення робочого простору

# Маса тіла (кг) шести ссавців — від миші до слона
M <- c(0.02, 0.3, 4, 70, 300, 4000)

# Інтенсивність обміну (Вт) за законом Клайбера: B = 3.4 * M^0.75
B <- 3.4 * M^0.75
#   Знак ^ означає «у степені». M^0.75 — маса в степені 0.75.
#   R рахує це одразу для всіх шести тварин.

B
#   Подивіться, що вийшло: шість чисел, по одному на кожну тварину.


# ---------------------------------------------------------------------
#   Два графіки поруч
# ---------------------------------------------------------------------

par(mfrow = c(1, 2))
#   Ділимо вікно графіків на дві частини: 1 рядок, 2 стовпці.

plot(M, B, pch = 16)
#   Звичайні координати. pch = 16 — суцільні точки.
#   Виходить крива.

plot(M, B, pch = 16, log = "xy")
#   Те саме, але з log = "xy": логарифмічна шкала на обох осях.
#   Виходить пряма.

par(mfrow = c(1, 1))
#   Повертаємо параметри графіки за умовчанням.

# ---------------------------------------------------------------------
#   Нахил прямої — по першій і останній точках
# ---------------------------------------------------------------------

(log10(B[6]) - log10(B[1])) / (log10(M[6]) - log10(M[1]))

#   B[1] — перше число у векторі B (миша), B[6] — шосте (слон).
#   log10 — десятковий логарифм.
#
#   Нахил = на скільки піднялися / на скільки просунулися вправо.
#   Угорі — різниця по вертикалі, унизу — різниця по горизонталі.
#
#   Має вийти рівно 0.75 — той показник степеня, який ми самі заклали.