File:Population of ancient tehuacan valley 1.svg
Original file (SVG file, nominally 914 × 454 pixels, file size: 31 KB)
Captions
Captions
Summary
[edit]
| Description |
English: Ancient population growth of ancient Tehuacán valley afret McNeish |
| Date | |
| Source | Own work |
| Author | Merikanto |
Sources:
Richard McNeish: The Prehistory of the Tehuacan Valley Vol. 1-5, 1967–1972
https://www.academia.edu/36910072/The_Prehistory_of_the_Tehuacan_Valley_Vol_1_MacNeish_Byers_y_otros
https://users.pop.umn.edu/~rmccaa/mxpoprev/cambridg3.htm
The Peopling of Mexico from Origins to Revolution
Robert McCaa
rmccaa@tc.umn.edu
preliminary draft (notes omitted from web version) for:
Richard Steckel & Michael Haines (eds.),
The Population History of North America
(Cambridge University Press).
December 8, 1997
- tehuacan valley after mcneish and also mccaa
import numpy as np
import math
import matplotlib.pyplot as plt
times1=np.array([-11000,-6800,-5000,-3000,-1500,-200, 700])
times2=np.array([-6800,-5000,-3000,-1500,-200,700,1521])
times3=(times1+times2)/2
popula1=np.array([12,48,120,480,1800,12000,60000])
popula2=2*popula1
popula3=(popula1+popula2)/2
populacoeff1=np.array([1,4,10,40,150,1000,5000])
densitytimes1=np.array([-7000,-3400,-1000,-500,-100,700])
densitytimes2=np.array([-5000,-2300,-1000,-500,700, 700])
densitytimes3=(densitytimes1+densitytimes2)/2
denisties_per_100km2=np.array([2.2,14, 43, 165, 1100,3600])
print(times3)
p1=popula3[0] ## population at time zero
p2=popula3[6] ## population at time colonization
ta1=times3[0] ## presumed time of arrival paleoindians
ta2=times3[6] ## time just before arrival of Columbus
deltat=100 ## time step
- print (math.log(2.72))
- quit(-1)
td1=abs(ta2-ta1)
pd1=p2/p1
rangelen1=int(td1/(deltat))
r1=math.log(pd1)/td1
t_two=math.log(2)/r1
r2=int(r1*100000)/1000
- print (pd1)
- print (td1)
print(" td1 ",td1)
print(" r %",r2)
print(" T2 ",int(t_two))
- quit(-1)
- r1=0.000656
tees=[]
pops=[]
tb1=0
taa1=ta1
print(ta1,p1)
for n in range(0,rangelen1):
tb1=tb1+deltat
taa1=taa1+deltat
p2=p1*math.exp(tb1*r1)
p2=int(p2/10)*10
tees.append(taa1)
pops.append(p2)
#print(taa1,tb1,p2)
print(taa1,p2)
plt.plot(tees,pops, ":")
plt.plot(times3, popula3, "r-")
plt.plot(times3, popula1, "g-")
plt.plot(times3, popula2, "b-")
MEDIUM_SIZE = 12
plt.rcParams.update({'font.size': MEDIUM_SIZE,})
plt.tick_params(labelsize=MEDIUM_SIZE)
- plt.rc('font', size=SMALL_SIZE) # controls default text sizes
plt.rc('xtick', labelsize=MEDIUM_SIZE) # fontsize of the tick labels
plt.rc('ytick', labelsize=MEDIUM_SIZE)
plt.suptitle("Population of Tehuacán valley")
plt.title("estimated")
plt.xlabel("Year BC/AD", size=MEDIUM_SIZE)
plt.ylabel("Population", size=MEDIUM_SIZE)
- annot1="Initial population at "+str(ta1)+ " BC "+str(p1)
- plt.text(-10000, 4.0e6, "Assumptions", fontsize=12)
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.
- share alike – If you remix, transform, or build upon the material, you must distribute your contributions under the same or compatible license as the original.
File history
Click on a date/time to view the file as it appeared at that time.
| Date/Time | Thumbnail | Dimensions | User | Comment | |
|---|---|---|---|---|---|
| current | 14:27, 10 November 2022 | | 914 × 454 (31 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.
