svg("Bev\u00F6lkerungsentwicklung Deutschland.svg", width = 15, height = 9)

year <- seq(1901, 2022)
popdevelopment <- data.frame(year)
########################################
rownames(popdevelopment) <- year
popdevelopment[50:122, "born"] <- c(1116701, 1106380, 1105084, 1095029, 1109743, 1113408, 1137169, 1165555, 1175870, 1243922, 1261614, 1313505, 1316534, 1355595, 1357304, 1325386, 1318303, 1272276, 1214968, 1142366, 1047737, 1013396, 901657, 815969, 805500, 782310, 798334, 805496, 808619, 817217, 865789, 862100, 861275, 827933, 812292, 813803, 848232, 867969, 892993, 880459, 905675, 830019, 809114, 798447, 769603, 765221, 796013, 812173, 785034, 770744, 766999, 734475, 719250, 706721, 705622, 685795, 672724, 684862, 682514, 665126, 677947, 662685, 673544, 682069, 714927, 737575, 792141, 784901, 787523, 778090, 773144, 795492, 738856)
comment(popdevelopment$born) <- "Sources: 1950-  Federal Statistical Office Germany Tables 12612-0001"
########################################
popdevelopment[50:122, "died"] <- c(748329, 752697, 767639, 790654, 775291, 795938, 812111, 840195, 818418, 835402, 876721, 850300, 878814, 895070, 870319, 907882, 911984, 914417, 976521, 988092, 975664, 965623, 965689, 962988, 956573, 989649, 966873, 931155, 955550, 944474, 952371, 954436, 943832, 941032, 917299, 929649, 925426, 901291, 900627, 903441, 921445, 911245, 885443, 897270, 884661, 884588, 882843, 860389, 852382, 846330, 838797, 828541, 841686, 853946, 818271, 830227, 821627, 827155, 844439, 854544, 858768, 852328, 869582, 893825, 868356, 925200, 910902, 932272, 954874, 939520, 985572, 1023687, 1066317)
comment(popdevelopment$died) <- "Sources: 1950-  Federal Statistical Office Germany Tables 12613-0002"
########################################


years <- 1950:2022
par(mar = c(2.5, 6.5, 4.0, 0)) #bottom, left, top, right
plot(years, popdevelopment$born[50:122], ylim = c(500000, 1400000), bty = "l", type = "l", yaxs = "i", xlab = "", ylab = "", pch = 16, cex = 2, lab = c(10, 10, 1), col = "forestgreen", main = "Bev\u00F6lkerungsentwicklung in Deutschland\n(1950-2022)", cex.main = 1.5, sub = "", lwd = 4, yaxt = "n", xaxt = "n", lty = 1)
lines(years, popdevelopment$died[50:122], lty = 1, type = "l", pch = 15, cex = 2, lwd = 4, col = "navyblue")

axis(1, at = seq(1950, 2022, by = 5), las = 0, tck = 1, lty = 2, col = "#c9c9c9", lwd = 1, cex.axis = 1.5)
axis(2, at = seq(500000, 1400000, by = 100000), las = 2, tck = 1, lty = 2, col = "#c9c9c9", lwd = 1, cex.axis = 1.5)
axis(2, at = seq(0, 0, by = 1), tck = 1, las = 2, lty = 2, col = "#000000", cex.axis = 1.5)

legend(x = "topright",      # Position
    inset = 0.02,
    legend = c("Gestorben", "Geboren"),  # Legend texts
    lty = c(1, 1),
    col = c("navyblue", "forestgreen"),      # Line colors
    lwd = 4,
    cex = 1.5)
dev.off()