File:Population of south mesopotamia towns 5000-2900 bc 1.png
Original file (1,071 × 704 pixels, file size: 100 KB, MIME type: image/png)
Captions
Captions
Summary
[edit]
| Description |
English: Population of South Mesopotamia towns 5000-2900 BC |
| Date | |
| Source | Own work |
| Author | Merikanto |
Source of data of population of towns are maps in
- I made a map-thing
- https://www.nma-fallout.com/threads/i-made-a-map-thing.218744/
Python3 script to create this graphics
import math
import numpy as np
import matplotlib.pyplot as plt
- I made a map-thing
- https://www.nma-fallout.com/threads/i-made-a-map-thing.218744/
- year sigmapop biggestpop numcenters
- 5000 8500 6500 3
- 4700 19500 8500 4
- 3600 81000 14000 12
- 3300 126000 40000 12
- 3100 131000 45000 12
- 2900 264000 80000 12
years=np.array([-5000,-4700,-3600,-3300,-3100,-2900])
populations=np.array([8500,19500,81000,126000,131000,264000])
biggestplaces=np.array([6500,8500,14000,40000,45000,80000])
numplaces=np.array([3,4,12,12,12,12])
population0=populations[0]
year0=years[0]
n1=0
n2=5
t1=years[n1]
t2=years[n2]
p1=populations[n1]
p2=populations[n2]
dt1=t2-t1
dp1=p2/p1
r1=math.log(dp1)/dt1
doublingtime1=math.log(2)/r1
- r1=0.0015
begintime1=-5000
endtime1=-2900
numtimes1=100
print(r1, doublingtime1)
annot2="Growth rate is "+str(round((r1*100),2)) + " %\n"+ "Doubling time is "+ str( round((doublingtime1)) ) + " years"
annot1="Initial population "+str(round((population0),2)) + " \n at year "+str(round((abs(year0)),0))+ " BC"
time1=np.linspace(begintime1, endtime1, numtimes1)
reltime1=np.copy(time1)-time1[0]
pop1=populations[0]*np.exp(r1*reltime1)
plt.title("Population growth of South Mesopotamia", size=22)
plt.xlabel("Years (negative is BC)" , fontsize=15)
plt.ylabel("Population of South Mesopotamia" , fontsize=15 )
plt.xticks(fontsize=15)
plt.yticks(fontsize=15)
plt.text(-5000, 160000, annot1, color="#00003f", alpha=0.7, fontsize=17)
plt.text(-5000, 120000, annot2, color="#003f00", alpha=0.7, fontsize=17)
plt.plot(years, populations,lw=5,linestyle="-" ,alpha=0.9, label="population data")
plt.plot(years, biggestplaces,lw=7,linestyle="--" ,color="#003f00", alpha=0.9, label="population of biggest town")
plt.plot(time1, pop1, lw=6,linestyle="--" ,alpha=0.6,label="exponent curve")
plt.legend(fontsize=16)
plt.grid()
plt.show()
Licensing
[edit]
- You are free:
- to share – to copy, distribute and transmit the work
- to remix – to adapt the work
- Under the following conditions:
- attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
File history
Click on a date/time to view the file as it appeared at that time.
| Date/Time | Thumbnail | Dimensions | User | Comment | |
|---|---|---|---|---|---|
| current | 08:22, 8 July 2024 | | 1,071 × 704 (100 KB) | Merikanto (talk | contribs) | Uploaded own work with UploadWizard |
You cannot overwrite this file.
File usage on Commons
There are no pages that use this file.
