svg("Population of Germany.svg", width = 15, height = 9)

year <- seq(1901, 2021)
population <- data.frame(year)
################################################################################
rownames(population) <- year
population[c(50:121), "humans"] <- c(
    (18.388 + 50.958125), # 1950
    (18.351 + 51.434777),
    (18.328 + 51.863761),
    (18.178 + 52.453806),
    (18.059 + 52.943295),
    (17.944 + 53.517683),
    (17.718 + 53.339626),
    (17.517 + 54.064365),
    (17.355 + 54.719159),
    (17.298 + 55.257088),
    (17.241 + 55.958321), # 1960
    (17.125 + 56.589148),
    (17.102 + 57.247246),
    (17.155 + 57.864509),
    (16.983 + 58.587451),
    (17.020 + 59.296591),
    (17.058 + 59.792934),
    (17.082 + 59.948474),
    (17.084 + 60.463033),
    (17.076 + 61.194591),
    (17.058 + 61.001164), # 1970
    (17.061 + 61.502503),
    (17.043 + 61.809378),
    (16.980 + 62.101369),
    (16.925 + 61.991475),
    (16.850 + 61.644624),
    (16.786 + 61.441996),
    (16.765 + 61.352745),
    (16.756 + 61.321663),
    (16.745 + 61.439342),
    (16.737 + 61.657945), # 1980
    (16.736 + 61.712689),
    (16.697 + 61.546101),
    (16.699 + 61.306669),
    (16.671 + 61.049256),
    (16.644 + 61.020474),
    (16.624 + 61.140461),
    (16.641 + 61.238079),
    (16.666 + 61.715103),
    (16.614 + 62.679035),
    79.753227,  # 1990
    80.274564,
    80.974632,
    81.338093,
    81.538603,
    81.817499,
    82.012162,
    82.057379,
    82.037011,
    82.163475,
    82.259540, # 2000
    82.440309,
    82.536680,
    82.531671,
    82.500849,
    82.437995,
    82.314906,
    82.217837,
    82.002356,
    81.802257,
    81.751602, # 2010
    80.327900,
    80.523746,
    80.767463,
    81.197537,
    82.175684,
    82.521653,
    82.792351,
    83.019213,
    83.166711,
    83.155031, # 2020
    83.237124
)
comment(population$humans) <- "Sources:
    East German population: -Special series with contributions for the territory of the former GDR- Booklet 3, Population Statistical Surveys 1946 to 1989, released May 1993 https://www.statistischebibliothek.de/mir/receive/DEMonografie_mods_00002289 (page 24);
    West German population since 1950 and combined populations since 1990: Federal Statistical Office Germany Table 12411-0001"

########################################
years <- 1950:2021
par(mar = c(2.7, 4.5, 2.5, 0)) #bottom, left, top, right
plot(
    years,
    population$humans[50:121],
    ylim = c(65, 85),
    bty = "l",
    type = "b",
    yaxs = "i",
    xlab = "",
    ylab = "",
    pch = 15,
    lab = c(10, 10, 1),
    col = "red",
    main = "Population of Germany",
    cex.main = 2,
    sub = "",
    lwd = 2,
    yaxt = "n",
    xaxt = "n"
)

axis(
    1,
    at = seq(1950, 2021, by = 5),
    las = 0,
    tck = 1,
    lty = 1,
    col = "#c9c9c9",
    lwd = 1
)
axis(
    2,
    at = seq(65, 85, by = 1),
    las = 2,
    tck = 1,
    lty = 1,
    col = "#c9c9c9",
    lwd = 1
)
mtext( # Y Axis
    side = 2,
    text = "million",
    cex = 1.3,
    line = 2.5
)
dev.off()