(* Content-type: application/mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 6.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       145,          7]
NotebookDataLength[     29048,        856]
NotebookOptionsPosition[     26271,        765]
NotebookOutlinePosition[     26607,        780]
CellTagsIndexPosition[     26564,        777]
WindowFrame->Normal
ContainsDynamic->False*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{
   RowBox[{
   "The", " ", "following", " ", "data", " ", "was", " ", "taken", " ", 
    "from", " ", "Livi"}], "-", 
   RowBox[{"Bacci", " ", "Massimo", 
    RowBox[{".2007", ".", " ", "A"}], " ", "Concise", " ", "History", " ", 
    "of", " ", "World", " ", "Population", " ", "4", "th", " ", 
    RowBox[{"ed", ".", "Blackwll"}], " ", 
    RowBox[{"Publishing", ".", "Malden"}], " ", "MA", " ", "Table", " ", 
    "1.3", " ", 
    RowBox[{"pp", ".", " ", "26."}]}]}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.4240242885528584`*^9, 3.424024313685807*^9}, {
   3.424024349030258*^9, 3.424024360211664*^9}, 3.4240244009821777`*^9, {
   3.429087423182743*^9, 3.429087507544731*^9}, {3.429088353426104*^9, 
   3.429088354783078*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.429087417585534*^9, 3.429087417611125*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"Asian", " ", "population", " ", "data"}], " ", "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Asia", " ", "=", " ", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "400.`"}], ",", "95.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.`", ",", "170.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"200.`", ",", "158.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"600.`", ",", "134.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1000.`", ",", "152.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1200.`", ",", "258.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1340.`", ",", "238.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1400.`", ",", "201.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1500.`", ",", "245.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1600.`", ",", "338.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1700.`", ",", "433.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1750.`", ",", "500.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1800.`", ",", "631.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1850.`", ",", "790.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1900.`", ",", "903.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1950.`", ",", "1376.`"}], "}"}]}], "}"}], "}"}]}], 
   ";"}]}]], "Input",
 CellChangeTimes->{{3.4240244145938673`*^9, 3.42402449417234*^9}, {
   3.4240246098534327`*^9, 3.4240246295537853`*^9}, {3.42402468162136*^9, 
   3.424024686834099*^9}, 3.4240325295000153`*^9, {3.427289256238428*^9, 
   3.427289284726506*^9}, {3.427289318999168*^9, 3.427289321058813*^9}, {
   3.429087541649441*^9, 3.429087551961969*^9}, {3.4290876684832*^9, 
   3.429087684349862*^9}, 3.4290877184320602`*^9}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.4290876797231007`*^9, 3.429087679742828*^9}}],

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"(*", " ", 
   RowBox[{"European", " ", "Population"}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.4290875986469927`*^9, 3.429087611445505*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Euro", "=", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "400.`"}], ",", "19.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0.`", ",", "31.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"200.`", ",", "44.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"600.`", ",", "22.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1000.`", ",", "30.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1200.`", ",", "49.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1340.`", ",", "74.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1400.`", ",", "52.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1500.`", ",", "67.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1600.`", ",", "89.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1700.`", ",", "95.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1750.`", ",", "111.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1800.`", ",", "146.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1850.`", ",", "209.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1900.`", ",", "295.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1950.`", ",", "393.`"}], "}"}]}], "}"}], "}"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.424032545825543*^9, 3.424032549871756*^9}, {
  3.4240325828270206`*^9, 3.424032590703824*^9}, {3.427289366662715*^9, 
  3.427289368083634*^9}, {3.429087734064307*^9, 3.429087739471101*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.429087731412033*^9, 3.429087731441736*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"Eurasion", " ", "Population"}], " ", "*)"}], 
  "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.429087615137693*^9, 3.42908762845052*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"USSR", "=", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "400.`"}], ",", "13.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0.`", ",", "12.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"200.`", ",", "13.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"600.`", ",", "11.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1000.`", ",", "13.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1200.`", ",", "17.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1340.`", ",", "16.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1400.`", ",", "13.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1500.`", ",", "17.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1600.`", ",", "22.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1700.`", ",", "30.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1750.`", ",", "35.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1800.`", ",", "49.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1850.`", ",", "79.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1900.`", ",", "127.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1950.`", ",", "182.`"}], "}"}]}], "}"}], "}"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.4240331565090103`*^9, 3.424033165810562*^9}, 
   3.427289379472474*^9, {3.429087752433784*^9, 3.429087756839623*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.429087749820283*^9, 3.4290877498530617`*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"African", " ", "Population"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.429087795050838*^9, 3.429087812005619*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Africa", "=", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "400.`"}], ",", "17.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0.`", ",", "26.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"200.`", ",", "30.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"600.`", ",", "24.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1000.`", ",", "39.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1200.`", ",", "48.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1340.`", ",", "80.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1400.`", ",", "68.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1500.`", ",", "87.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1600.`", ",", "113.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1700.`", ",", "107.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1750.`", ",", "104.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1800.`", ",", "102.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1850.`", ",", "102.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1900.`", ",", "138.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1950.`", ",", "224.`"}], "}"}]}], "}"}], "}"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.424033404255427*^9, 3.424033432302999*^9}, 
   3.42728938535118*^9, {3.429087772265606*^9, 3.429087777545369*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.427290337012286*^9, 3.427290413214183*^9}, {
   3.427290455161539*^9, 3.4272904837836533`*^9}, 3.427290526188349*^9, {
   3.4272905578490334`*^9, 3.427290596506034*^9}, 3.427290663971134*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"America", "=", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "400.`"}], ",", "8.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0.`", ",", "12.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"200.`", ",", "11.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"600.`", ",", "16.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1000.`", ",", "18.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1200.`", ",", "26.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1340.`", ",", "32.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1400.`", ",", "39.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1500.`", ",", "42.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1600.`", ",", "13.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1700.`", ",", "12.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1750.`", ",", "18.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1800.`", ",", "24.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1850.`", ",", "59.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1900.`", ",", "165.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1950.`", ",", "332.`"}], "}"}]}], "}"}], "}"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.4240335780729723`*^9, 3.4240335871829233`*^9}, 
   3.427289388552115*^9, {3.4290878231461983`*^9, 3.429087829041749*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.4290878209895163`*^9, 3.429087821024727*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"Oceania", " ", "Population"}], " ", "*)"}], "\[IndentingNewLine]",
   
  RowBox[{
   RowBox[{"Ocean", "=", 
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{
         RowBox[{"-", "400.`"}], ",", "1.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"0.`", ",", "1.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"200.`", ",", "1.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"600.`", ",", "1.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1000.`", ",", "1.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1200.`", ",", "2.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1340.`", ",", "2.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1400.`", ",", "2.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1500.`", ",", "3.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1600.`", ",", "3.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1700.`", ",", "3.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1750.`", ",", "3.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1800.`", ",", "2.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1850.`", ",", "2.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1900.`", ",", "6.`"}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"1950.`", ",", "13.`"}], "}"}]}], "}"}], "}"}]}], 
   ";"}]}]], "Input",
 CellChangeTimes->{{3.424033688753791*^9, 3.42403370035121*^9}, {
   3.424033874524856*^9, 3.424033877834593*^9}, 3.427289391192255*^9, {
   3.429087847962233*^9, 3.4290878865298157`*^9}}],

Cell[BoxData[
 RowBox[{"(*", " ", 
  RowBox[{"Total", " ", "Population"}], " ", "*)"}]], "Input",
 CellChangeTimes->{{3.429087923999133*^9, 3.429087939536084*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Tot", "=", 
   RowBox[{"{", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "400.`"}], ",", "153.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"0.`", ",", "252.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"200.`", ",", "257.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"600.`", ",", "208.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1000.`", ",", "253.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1200.`", ",", "400.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1340.`", ",", "442.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1400.`", ",", "375.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1500.`", ",", "461.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1600.`", ",", "578.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1700.`", ",", "680.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1750.`", ",", "771.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1800.`", ",", "954.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1850.`", ",", "1241.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1900.`", ",", "1634.`"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"1950.`", ",", "2520.`"}], "}"}]}], "}"}], "}"}]}], 
  ";"}]], "Input",
 CellChangeTimes->{{3.424033889235499*^9, 3.424033898289261*^9}, 
   3.4272893938562117`*^9, {3.429087904771534*^9, 3.429087911530466*^9}}],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.429087901791236*^9, 3.429087901824355*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Tot", "[", 
  RowBox[{"[", "1", "]"}], "]"}]], "Input",
 CellChangeTimes->{{3.424033909419714*^9, 3.424033918201599*^9}, {
  3.429087972553999*^9, 3.429088014202675*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "400.`"}], ",", "153.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.`", ",", "252.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"200.`", ",", "257.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"600.`", ",", "208.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1000.`", ",", "253.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1200.`", ",", "400.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1340.`", ",", "442.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1400.`", ",", "375.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1500.`", ",", "461.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1600.`", ",", "578.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1700.`", ",", "680.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1750.`", ",", "771.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1800.`", ",", "954.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1850.`", ",", "1241.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1900.`", ",", "1634.`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"1950.`", ",", "2520.`"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{{3.42908799960433*^9, 3.429088014751315*^9}, {
   3.4290883461663427`*^9, 3.429088371687114*^9}, 3.429088603070915*^9}]
}, Open  ]],

Cell[BoxData[""], "Input",
 CellChangeTimes->{{3.429070397147241*^9, 3.429070415641631*^9}, 
   3.429087972555079*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{
   "ftot", " ", "is", " ", "an", " ", "interpolation", " ", "of", " ", "the", 
    " ", "total", " ", "population", " ", "curve", " ", "above"}], " ", 
   "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", 
  RowBox[{"ftot", "=", 
   RowBox[{"Interpolation", "[", 
    RowBox[{"Tot", "[", 
     RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.4290882642165747`*^9, 3.429088301925995*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"InterpolatingFunction", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "400.`"}], ",", "1950.`"}], "}"}], "}"}], 
    ",", "\<\"<>\"\>"}], "]"}],
  False,
  Editable->False]], "Output",
 CellChangeTimes->{{3.4290883462705517`*^9, 3.429088371751124*^9}, 
   3.429088603114689*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
  "try", " ", "some", " ", "values", " ", "here", " ", "for", " ", "ftot", 
   " ", "to", " ", "see", " ", "how", " ", "good", " ", "of", " ", "an", " ", 
   "iterpolation", " ", "it", " ", "is"}], "*)"}]], "Input",
 CellChangeTimes->{{3.429088436610135*^9, 3.42908845303552*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{"ftot", "[", "1825", "]"}]}]], "Input",
 CellChangeTimes->{{3.429088385439399*^9, 3.429088426887827*^9}, {
  3.429088458569296*^9, 3.429088473876491*^9}}],

Cell[BoxData["1084.375`"], "Output",
 CellChangeTimes->{{3.429088400352502*^9, 3.429088427591873*^9}, {
   3.4290884621095743`*^9, 3.4290884755809717`*^9}, 3.42908860316182*^9}]
}, Open  ]],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"Numberofbirthshigh", "=", "7"}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "Birth", " ", "rate", " ", "of", " ", "females", " ", "historically", " ", 
    "High", " ", "value"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Numberofbirthslow", "=", "4"}], ";"}], 
  RowBox[{"(*", " ", 
   RowBox[{
   "Birth", " ", "rate", " ", "of", " ", "females", " ", "historically", " ", 
    "Low", " ", "value"}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"deathratefemales", "=", 
    RowBox[{"1", "/", "4"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"Infant", " ", "mortality", " ", "rate", " ", "females"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"deathratemales", "=", 
    RowBox[{"(", 
     RowBox[{"deathratefemales", "*", "1.2"}], ")"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
   "compariative", " ", "rate", " ", "of", " ", "male", " ", "mortality"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"sexratiofemale", "=", 
    RowBox[{"103", "/", "107"}]}], ";"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"Historical", " ", "sex", " ", "ratio", " ", "at", " ", "birth"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"livetobreadingage", "=", 
   RowBox[{"600", "/", "1000"}]}], ";", " ", 
  RowBox[{"(*", 
   RowBox[{
   "proportion", " ", "that", " ", "live", " ", "to", " ", "childbearing", 
    " ", "age"}], " ", "*)"}], "\[IndentingNewLine]"}]}], "Input",
 CellChangeTimes->{{3.4290678067241583`*^9, 3.429067880693791*^9}, {
   3.4290679498552713`*^9, 3.429067978919648*^9}, {3.429068583045089*^9, 
   3.4290686020578938`*^9}, {3.429088054458716*^9, 3.429088170373909*^9}, 
   3.4290883673565598`*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Highbirthratemale", "=", 
   RowBox[{
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{"ftot", "[", "i", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"-", "400"}], ",", "1950", ",", "25"}], "}"}]}], "]"}], "*", 
    "sexratiofemale", "*", "Numberofbirthshigh", "*", "deathratemales", "*", 
    "livetobreadingage"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.4290676757074537`*^9, 3.429067798951816*^9}, {
  3.429067894343157*^9, 3.42906794545778*^9}, {3.429067983921072*^9, 
  3.4290680028973494`*^9}, {3.429068048365182*^9, 3.429068060690133*^9}, {
  3.4290681023898983`*^9, 3.429068106667297*^9}, {3.4290686126859293`*^9, 
  3.429068614123827*^9}}],

Cell[CellGroupData[{

Cell[BoxData["Highbirthratemale"], "Input",
 CellChangeTimes->{{3.429068022613187*^9, 3.429068026253207*^9}}],

Cell[BoxData["45846.35134715302`"], "Output",
 CellChangeTimes->{{3.4290680273180227`*^9, 3.429068067050695*^9}, {
   3.429068114187462*^9, 3.42906812049366*^9}, 3.4290683678472223`*^9, 
   3.429068640296093*^9, 3.429069040656645*^9, {3.429088346405748*^9, 
   3.4290883718455667`*^9}, 3.429088603235285*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"Lowbirthratemale", "=", 
   RowBox[{
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{"ftot", "[", "i", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"-", "400"}], ",", "1950", ",", "25"}], "}"}]}], "]"}], "*", 
    "sexratiofemale", "*", "Numberofbirthslow", "*", "deathratemales", "*", 
    "livetobreadingage"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.429068146337947*^9, 3.429068189204238*^9}, 
   3.429068625336152*^9}],

Cell[CellGroupData[{

Cell[BoxData["Lowbirthratemale"], "Input",
 CellChangeTimes->{{3.429068194003677*^9, 3.4290682023874617`*^9}}],

Cell[BoxData["26197.91505551601`"], "Output",
 CellChangeTimes->{
  3.429068203081582*^9, 3.429068367950231*^9, 3.429068640402734*^9, 
   3.429069040769474*^9, {3.429088346492078*^9, 3.429088371917233*^9}, 
   3.429088603289366*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"Highbirthratefemale", "=", 
   RowBox[{
    RowBox[{"Sum", "[", 
     RowBox[{
      RowBox[{"ftot", "[", "i", "]"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"-", "400"}], ",", "1950", ",", "25"}], "}"}]}], "]"}], "*", 
    "sexratiofemale", "*", "Numberofbirthshigh", "*", "deathratefemales", "*",
     "livetobreadingage"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.429068232039505*^9, 3.42906825615728*^9}, {
   3.429068296534577*^9, 3.42906829666225*^9}, 3.429068631783618*^9}],

Cell[CellGroupData[{

Cell[BoxData["Highbirthratefemale"], "Input",
 CellChangeTimes->{{3.429068284026382*^9, 3.4290682911726923`*^9}}],

Cell[BoxData["38205.29278929418`"], "Output",
 CellChangeTimes->{
  3.4290683044111233`*^9, 3.429068367980124*^9, 3.429068640453198*^9, 
   3.429069040853568*^9, {3.429088346585083*^9, 3.429088371996121*^9}, 
   3.429088603359001*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Lowbirthratefemale", "=", 
  RowBox[{
   RowBox[{"Sum", "[", 
    RowBox[{
     RowBox[{"ftot", "[", "i", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"-", "400"}], ",", "1950", ",", "25"}], "}"}]}], "]"}], "*", 
   "sexratiofemale", "*", "Numberofbirthslow", "*", "deathratefemales", "*", 
   "livetobreadingage"}]}]], "Input",
 CellChangeTimes->{{3.4290683323849792`*^9, 3.4290683423346567`*^9}, 
   3.429068629511086*^9, 3.429088591212325*^9}],

Cell[BoxData["21831.595879596676`"], "Output",
 CellChangeTimes->{3.42908860340716*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ExcessMalesHigh", "=", 
  RowBox[{"Highbirthratemale", "-", "Highbirthratefemale"}]}]], "Input",
 CellChangeTimes->{{3.429068377772393*^9, 3.429068406914118*^9}, {
  3.4290684687201567`*^9, 3.4290684690794363`*^9}}],

Cell[BoxData["7641.058557858836`"], "Output",
 CellChangeTimes->{
  3.429068409846486*^9, 3.429068640545321*^9, 3.429069040992423*^9, {
   3.429088346739855*^9, 3.429088372136366*^9}, 3.429088603440412*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ExcessMalesLow", "=", 
  RowBox[{"Lowbirthratemale", "-", "Lowbirthratefemale"}]}]], "Input",
 CellChangeTimes->{{3.4290684577132*^9, 3.429068512540059*^9}}],

Cell[BoxData["4366.319175919332`"], "Output",
 CellChangeTimes->{
  3.429068514247933*^9, 3.4290686405832167`*^9, 3.429069041043438*^9, {
   3.429088346885795*^9, 3.429088372203681*^9}, 3.429088603490398*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ExcessMalesHighCorrected", "=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"Highbirthratemale", "-", "Highbirthratefemale"}], ")"}], "*", 
   "Correctionfactor"}]}]], "Input",
 CellChangeTimes->{{3.429069143206047*^9, 3.429069198280048*^9}}],

Cell[BoxData["14756.915887850468`"], "Output",
 CellChangeTimes->{{3.42906919161281*^9, 3.429069200343747*^9}, {
   3.429088346942313*^9, 3.429088372273593*^9}, 3.429088603541271*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ExcessMalesLowCorrected", "=", 
  RowBox[{
   RowBox[{"(", 
    RowBox[{"Lowbirthratemale", "-", "Lowbirthratefemale"}], ")"}], "*", 
   "Correctionfactor"}]}]], "Input",
 CellChangeTimes->{{3.429069223568397*^9, 3.42906923874151*^9}}],

Cell[BoxData["8432.523364485976`"], "Output",
 CellChangeTimes->{
  3.4290692400294123`*^9, {3.429088347009157*^9, 3.4290883723403*^9}, 
   3.42908860360806*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{
   "People", " ", "going", " ", "to", " ", "celstial", " ", "kingdom"}], 
   "*)"}], "\[IndentingNewLine]"}]], "Input",
 CellChangeTimes->{{3.429070893239859*^9, 3.429070916197966*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"82000", "*", "0.002"}]], "Input",
 CellChangeTimes->{{3.429070922470756*^9, 3.429070941979672*^9}, {
  3.429071323473757*^9, 3.4290713265701637`*^9}, {3.4290713994765387`*^9, 
  3.429071401836995*^9}}],

Cell[BoxData["164.`"], "Output",
 CellChangeTimes->{{3.4290709333982353`*^9, 3.429070943526021*^9}, 
   3.429071327870224*^9, 3.429071404793665*^9, {3.4290883471038647`*^9, 
   3.4290883725163116`*^9}, 3.429088603703018*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"37000", "/", "2"}]], "Input",
 CellChangeTimes->{{3.4290710147417307`*^9, 3.429071018596841*^9}, {
  3.42907106119147*^9, 3.4290710624858303`*^9}}],

Cell[BoxData["18500"], "Output",
 CellChangeTimes->{
  3.429071019193819*^9, 3.429071063474491*^9, {3.42908834716241*^9, 
   3.429088372571039*^9}, 3.4290886038482103`*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{"NOTE", ":", " ", 
   RowBox[{
   "There", " ", "could", " ", "be", " ", "a", " ", "lot", " ", "of", " ", 
    "refinements", " ", "especially", " ", "by", " ", "applying", " ", "the", 
    " ", "parmeters", " ", "at", " ", "the", " ", "time", " ", "periouds", 
    " ", "as", " ", "they", " ", "change", " ", "through", " ", 
    RowBox[{"time", ".", " ", "I"}], " ", "just", " ", "dont", " ", "have", 
    " ", "the", " ", "data", " ", "so", " ", "assumped", " ", "that", " ", 
    "they", " ", "parameters", " ", "constant", " ", "through", " ", "time", 
    "*", "("}]}]}]], "Input",
 CellChangeTimes->{{3.429088501655545*^9, 3.4290885680245247`*^9}},
 EmphasizeSyntaxErrors->True]
},
WindowSize->{690, 589},
WindowMargins->{{-2, Automatic}, {Automatic, 3}},
FrontEndVersion->"6.0 for Mac OS X x86 (32-bit) (March 13, 2008)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[568, 21, 781, 16, 58, "Input"],
Cell[1352, 39, 92, 1, 27, InheritFromParent],
Cell[1447, 42, 1919, 49, 88, "Input"],
Cell[3369, 93, 94, 1, 27, InheritFromParent],
Cell[3466, 96, 205, 4, 43, "Input"],
Cell[3674, 102, 1531, 42, 73, "Input"],
Cell[5208, 146, 92, 1, 27, InheritFromParent],
Cell[5303, 149, 200, 5, 43, "Input"],
Cell[5506, 156, 1456, 41, 73, "Input"],
Cell[6965, 199, 94, 1, 27, InheritFromParent],
Cell[7062, 202, 165, 3, 27, "Input"],
Cell[7230, 207, 1460, 41, 73, "Input"],
Cell[8693, 250, 240, 3, 27, "Input"],
Cell[8936, 255, 1462, 41, 73, "Input"],
Cell[10401, 298, 94, 1, 27, InheritFromParent],
Cell[10498, 301, 1645, 46, 88, "Input"],
Cell[12146, 349, 163, 3, 27, "Input"],
Cell[12312, 354, 1472, 41, 73, "Input"],
Cell[13787, 397, 92, 1, 27, InheritFromParent],
Cell[CellGroupData[{
Cell[13904, 402, 194, 4, 27, "Input"],
Cell[14101, 408, 1298, 37, 73, "Output"]
}, Open  ]],
Cell[15414, 448, 118, 2, 27, "Input"],
Cell[CellGroupData[{
Cell[15557, 454, 476, 11, 58, "Input"],
Cell[16036, 467, 356, 12, 27, "Output"]
}, Open  ]],
Cell[16407, 482, 321, 6, 43, "Input"],
Cell[CellGroupData[{
Cell[16753, 492, 211, 4, 43, "Input"],
Cell[16967, 498, 177, 2, 27, "Output"]
}, Open  ]],
Cell[17159, 503, 1791, 48, 133, "Input"],
Cell[18953, 553, 719, 16, 43, "Input"],
Cell[CellGroupData[{
Cell[19697, 573, 109, 1, 27, "Input"],
Cell[19809, 576, 308, 4, 27, "Output"]
}, Open  ]],
Cell[20132, 583, 494, 13, 43, "Input"],
Cell[CellGroupData[{
Cell[20651, 600, 110, 1, 27, "Input"],
Cell[20764, 603, 232, 4, 27, "Output"]
}, Open  ]],
Cell[21011, 610, 544, 13, 43, "Input"],
Cell[CellGroupData[{
Cell[21580, 627, 113, 1, 27, "Input"],
Cell[21696, 630, 234, 4, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[21967, 639, 497, 12, 43, "Input"],
Cell[22467, 653, 87, 1, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22591, 659, 239, 4, 27, "Input"],
Cell[22833, 665, 206, 3, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23076, 673, 181, 3, 27, "Input"],
Cell[23260, 678, 208, 3, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23505, 686, 263, 6, 43, "Input"],
Cell[23771, 694, 183, 2, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23991, 701, 259, 6, 43, "Input"],
Cell[24253, 709, 161, 3, 27, "Output"]
}, Open  ]],
Cell[24429, 715, 238, 6, 43, "Input"],
Cell[CellGroupData[{
Cell[24692, 725, 225, 4, 27, "Input"],
Cell[24920, 731, 224, 3, 27, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[25181, 739, 171, 3, 27, "Input"],
Cell[25355, 744, 172, 3, 27, "Output"]
}, Open  ]],
Cell[25542, 750, 725, 13, 73, "Input"]
}
]
*)

(* End of internal cache information *)