diff --git a/code/SIRR Project.ipynb b/code/SIRR Project.ipynb new file mode 100644 index 00000000..1e2ed551 --- /dev/null +++ b/code/SIRR Project.ipynb @@ -0,0 +1,553 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 199, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [ + "%matplotlib inline\n", + "\n", + "from modsim import *\n", + "\n", + "from random import random\n" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [] + }, + { + "cell_type": "code", + "execution_count": 200, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [ + "init = State(Susceptible=349, Infected=1, Removed=0, Recovered=0) \n", + " \n" + ] + }, + { + "cell_type": "code", + "execution_count": 201, + "metadata": { + "collapsed": true + }, + "outputs": [], + "source": [ + "def make_system():\n", + " \n", + " init = State(Susceptible=349, Infected=1, Removed=0, Recovered=0)\n", + " infected=init.Infected\n", + " recovered=init.Recovered\n", + " removed=init.Removed\n", + " susceptible=init.Susceptible\n", + " \n", + " init_infection_exit = State(new_susceptible=0, new_infected=0, new_recovered=0, new_removed=0)\n", + " \n", + " t0=0\n", + " t_end=98\n", + " \n", + " \n", + " system = System(t0=t0, t_end=t_end, init=init, init_infection_exit=init_infection_exit,\n", + " susceptible=susceptible, infected=infected, recovered=recovered, removed=removed)\n", + " \n", + " return system" + ] + }, + { + "cell_type": "code", + "execution_count": 202, + "metadata": {}, + "outputs": [], + "source": [ + "system = make_system()" + ] + }, + { + "cell_type": "code", + "execution_count": 203, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
| \n", + " | value | \n", + "
|---|---|
| t0 | \n", + "0 | \n", + "
| t_end | \n", + "98 | \n", + "
| init | \n", + "Susceptible 349\n", + "Infected 1\n", + "Removed ... | \n", + "
| init_infection_exit | \n", + "new_susceptible 0\n", + "new_infected 0\n", + "new_... | \n", + "
| susceptible | \n", + "349 | \n", + "
| infected | \n", + "1 | \n", + "
| recovered | \n", + "0 | \n", + "
| removed | \n", + "0 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| Susceptible | \n", + "235 | \n", + "
| Infected | \n", + "115 | \n", + "
| Removed | \n", + "0 | \n", + "
| Recovered | \n", + "0 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| Susceptible | \n", + "235 | \n", + "
| Infected | \n", + "115 | \n", + "
| Removed | \n", + "0 | \n", + "
| Recovered | \n", + "0 | \n", + "
| \n", + " | Susceptible | \n", + "Infected | \n", + "Removed | \n", + "Recovered | \n", + "
|---|---|---|---|---|
| 0 | \n", + "349.0 | \n", + "1.0 | \n", + "0.0 | \n", + "0.0 | \n", + "
| 1 | \n", + "152.0 | \n", + "165.0 | \n", + "33.0 | \n", + "0.0 | \n", + "
| 2 | \n", + "104.0 | \n", + "171.0 | \n", + "75.0 | \n", + "0.0 | \n", + "
| 3 | \n", + "71.0 | \n", + "161.0 | \n", + "118.0 | \n", + "0.0 | \n", + "
| 4 | \n", + "45.0 | \n", + "151.0 | \n", + "154.0 | \n", + "0.0 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "10 | \n", + "
| wellesley | \n", + "2 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "10 | \n", + "
| wellesley | \n", + "2 | \n", + "
| babson | \n", + "0 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "3 | \n", + "
| wellesley | \n", + "9 | \n", + "
| babson | \n", + "0 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "-1 | \n", + "
| wellesley | \n", + "13 | \n", + "
| babson | \n", + "0 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "10 | \n", + "
| wellesley | \n", + "2 | \n", + "
| babson | \n", + "0 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "9 | \n", + "
| wellesley | \n", + "3 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "10 | \n", + "
| wellesley | \n", + "2 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "8 | \n", + "
| wellesley | \n", + "4 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "9 | \n", + "
| wellesley | \n", + "3 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "9 | \n", + "
| wellesley | \n", + "3 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "10 | \n", + "
| wellesley | \n", + "2 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "10 | \n", + "
| wellesley | \n", + "2 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "11 | \n", + "
| wellesley | \n", + "1 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "9 | \n", + "
| wellesley | \n", + "3 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| olin | \n", + "10 | \n", + "
| wellesley | \n", + "2 | \n", + "
| \n", + " | United States Census Bureau (2015)[18] | \n", + "Population Reference Bureau (1973–2015)[6] | \n", + "United Nations Department of Economic and Social Affairs (2015)[7] | \n", + "Maddison (2008)[8] | \n", + "HYDE (2007)[15] | \n", + "Tanton (1994)[9] | \n", + "Biraben (1980)[10] | \n", + "McEvedy & Jones (1978)[11] | \n", + "Thomlinson (1975)[12] | \n", + "Durand (1974)[13] | \n", + "Clark (1967)[14] | \n", + "
|---|---|---|---|---|---|---|---|---|---|---|---|
| Year | \n", + "\n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " |
| 1950 | \n", + "2557628654 | \n", + "2.516000e+09 | \n", + "2525149000 | \n", + "2.544000e+09 | \n", + "2.527960e+09 | \n", + "2.400000e+09 | \n", + "2.527000e+09 | \n", + "2.500000e+09 | \n", + "2.400000e+09 | \n", + "NaN | \n", + "2.486000e+09 | \n", + "
| 1951 | \n", + "2594939877 | \n", + "NaN | \n", + "2572850917 | \n", + "2.571663e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1952 | \n", + "2636772306 | \n", + "NaN | \n", + "2619292068 | \n", + "2.617949e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1953 | \n", + "2682053389 | \n", + "NaN | \n", + "2665865392 | \n", + "2.665959e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1954 | \n", + "2730228104 | \n", + "NaN | \n", + "2713172027 | \n", + "2.716927e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| \n", + " | United States Census Bureau (2015)[18] | \n", + "Population Reference Bureau (1973–2015)[6] | \n", + "United Nations Department of Economic and Social Affairs (2015)[7] | \n", + "Maddison (2008)[8] | \n", + "HYDE (2007)[15] | \n", + "Tanton (1994)[9] | \n", + "Biraben (1980)[10] | \n", + "McEvedy & Jones (1978)[11] | \n", + "Thomlinson (1975)[12] | \n", + "Durand (1974)[13] | \n", + "Clark (1967)[14] | \n", + "
|---|---|---|---|---|---|---|---|---|---|---|---|
| Year | \n", + "\n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " |
| 2011 | \n", + "6944055583 | \n", + "6.986951e+09 | \n", + "6997998760 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2012 | \n", + "7022349283 | \n", + "7.057075e+09 | \n", + "7080072417 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2013 | \n", + "7101027895 | \n", + "7.136796e+09 | \n", + "7162119434 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2014 | \n", + "7178722893 | \n", + "7.238184e+09 | \n", + "7243784000 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2015 | \n", + "7256490011 | \n", + "7.336435e+09 | \n", + "7349472000 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| \n", + " | census | \n", + "prb | \n", + "un | \n", + "maddison | \n", + "hyde | \n", + "tanton | \n", + "biraben | \n", + "mj | \n", + "thomlinson | \n", + "durand | \n", + "clark | \n", + "
|---|---|---|---|---|---|---|---|---|---|---|---|
| Year | \n", + "\n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " |
| 1950 | \n", + "2557628654 | \n", + "2.516000e+09 | \n", + "2525149000 | \n", + "2.544000e+09 | \n", + "2.527960e+09 | \n", + "2.400000e+09 | \n", + "2.527000e+09 | \n", + "2.500000e+09 | \n", + "2.400000e+09 | \n", + "NaN | \n", + "2.486000e+09 | \n", + "
| 1951 | \n", + "2594939877 | \n", + "NaN | \n", + "2572850917 | \n", + "2.571663e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1952 | \n", + "2636772306 | \n", + "NaN | \n", + "2619292068 | \n", + "2.617949e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1953 | \n", + "2682053389 | \n", + "NaN | \n", + "2665865392 | \n", + "2.665959e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1954 | \n", + "2730228104 | \n", + "NaN | \n", + "2713172027 | \n", + "2.716927e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1955 | \n", + "2782098943 | \n", + "NaN | \n", + "2761650981 | \n", + "2.769074e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1956 | \n", + "2835299673 | \n", + "NaN | \n", + "2811572031 | \n", + "2.822502e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1957 | \n", + "2891349717 | \n", + "NaN | \n", + "2863042795 | \n", + "2.879934e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1958 | \n", + "2948137248 | \n", + "NaN | \n", + "2916030167 | \n", + "2.939254e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1959 | \n", + "3000716593 | \n", + "NaN | \n", + "2970395814 | \n", + "2.995909e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1960 | \n", + "3043001508 | \n", + "NaN | \n", + "3026002942 | \n", + "3.041507e+09 | \n", + "3.042000e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1961 | \n", + "3083966929 | \n", + "NaN | \n", + "3082830266 | \n", + "3.082161e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1962 | \n", + "3140093217 | \n", + "NaN | \n", + "3141071531 | \n", + "3.135787e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "3.036000e+09 | \n", + "
| 1963 | \n", + "3209827882 | \n", + "NaN | \n", + "3201178277 | \n", + "3.201354e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1964 | \n", + "3281201306 | \n", + "NaN | \n", + "3263738832 | \n", + "3.266477e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1965 | \n", + "3350425793 | \n", + "NaN | \n", + "3329122479 | \n", + "3.333138e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1966 | \n", + "3420677923 | \n", + "NaN | \n", + "3397475247 | \n", + "3.402224e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "3.288000e+09 | \n", + "
| 1967 | \n", + "3490333715 | \n", + "NaN | \n", + "3468521724 | \n", + "3.471464e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1968 | \n", + "3562313822 | \n", + "NaN | \n", + "3541674891 | \n", + "3.543086e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1969 | \n", + "3637159050 | \n", + "NaN | \n", + "3616108749 | \n", + "3.615743e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1970 | \n", + "3712697742 | \n", + "NaN | \n", + "3691172616 | \n", + "3.691157e+09 | \n", + "3.710000e+09 | \n", + "NaN | \n", + "3.637000e+09 | \n", + "NaN | \n", + "3.600000e+09 | \n", + "3,600,000,000– 3,700,000,000 | \n", + "3.632000e+09 | \n", + "
| 1971 | \n", + "3790326948 | \n", + "NaN | \n", + "3766754345 | \n", + "3.769818e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1972 | \n", + "3866568653 | \n", + "NaN | \n", + "3842873611 | \n", + "3.846499e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1973 | \n", + "3942096442 | \n", + "NaN | \n", + "3919182332 | \n", + "3.922793e+09 | \n", + "3.923000e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "3.860000e+09 | \n", + "
| 1974 | \n", + "4016608813 | \n", + "NaN | \n", + "3995304922 | \n", + "3.997677e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1975 | \n", + "4089083233 | \n", + "NaN | \n", + "4071020434 | \n", + "4.070671e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "3.900000e+09 | \n", + "4.000000e+09 | \n", + "NaN | \n", + "NaN | \n", + "
| 1976 | \n", + "4160185010 | \n", + "NaN | \n", + "4146135850 | \n", + "4.141445e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1977 | \n", + "4232084578 | \n", + "NaN | \n", + "4220816737 | \n", + "4.213539e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1978 | \n", + "4304105753 | \n", + "NaN | \n", + "4295664825 | \n", + "4.286317e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1979 | \n", + "4379013942 | \n", + "NaN | \n", + "4371527871 | \n", + "4.363144e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| ... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "... | \n", + "
| 1986 | \n", + "4940571232 | \n", + "NaN | \n", + "4953376710 | \n", + "4.920968e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1987 | \n", + "5027200492 | \n", + "NaN | \n", + "5045315871 | \n", + "5.006672e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1988 | \n", + "5114557167 | \n", + "NaN | \n", + "5138214688 | \n", + "5.093306e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1989 | \n", + "5201440110 | \n", + "NaN | \n", + "5230000000 | \n", + "5.180540e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1990 | \n", + "5288955934 | \n", + "NaN | \n", + "5320816667 | \n", + "5.269029e+09 | \n", + "5.308000e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1991 | \n", + "5371585922 | \n", + "NaN | \n", + "5408908724 | \n", + "5.351922e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1992 | \n", + "5456136278 | \n", + "NaN | \n", + "5494899570 | \n", + "5.435722e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1993 | \n", + "5538268316 | \n", + "NaN | \n", + "5578865109 | \n", + "5.518127e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1994 | \n", + "5618682132 | \n", + "NaN | \n", + "5661086346 | \n", + "5.599396e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1995 | \n", + "5699202985 | \n", + "5.760000e+09 | \n", + "5741822412 | \n", + "5.681575e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1996 | \n", + "5779440593 | \n", + "NaN | \n", + "5821016750 | \n", + "5.762212e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1997 | \n", + "5857972543 | \n", + "5.840000e+09 | \n", + "5898688337 | \n", + "5.842122e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1998 | \n", + "5935213248 | \n", + "NaN | \n", + "5975303657 | \n", + "5.921366e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1999 | \n", + "6012074922 | \n", + "NaN | \n", + "6051478010 | \n", + "5.999622e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2000 | \n", + "6088571383 | \n", + "6.067000e+09 | \n", + "6127700428 | \n", + "6.076558e+09 | \n", + "6.145000e+09 | \n", + "NaN | \n", + "NaN | \n", + "5.750000e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2001 | \n", + "6165219247 | \n", + "6.137000e+09 | \n", + "6204147026 | \n", + "6.154791e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2002 | \n", + "6242016348 | \n", + "6.215000e+09 | \n", + "6280853817 | \n", + "6.231704e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2003 | \n", + "6318590956 | \n", + "6.314000e+09 | \n", + "6357991749 | \n", + "6.308364e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2004 | \n", + "6395699509 | \n", + "6.396000e+09 | \n", + "6435705595 | \n", + "6.374056e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2005 | \n", + "6473044732 | \n", + "6.477000e+09 | \n", + "6514094605 | \n", + "6.462987e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2006 | \n", + "6551263534 | \n", + "6.555000e+09 | \n", + "6593227977 | \n", + "6.540214e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2007 | \n", + "6629913759 | \n", + "6.625000e+09 | \n", + "6673105937 | \n", + "6.616689e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2008 | \n", + "6709049780 | \n", + "6.705000e+09 | \n", + "6753649228 | \n", + "6.694832e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2009 | \n", + "6788214394 | \n", + "6.809972e+09 | \n", + "6834721933 | \n", + "6.764086e+09 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2010 | \n", + "6866332358 | \n", + "6.892319e+09 | \n", + "6916183482 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2011 | \n", + "6944055583 | \n", + "6.986951e+09 | \n", + "6997998760 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2012 | \n", + "7022349283 | \n", + "7.057075e+09 | \n", + "7080072417 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2013 | \n", + "7101027895 | \n", + "7.136796e+09 | \n", + "7162119434 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2014 | \n", + "7178722893 | \n", + "7.238184e+09 | \n", + "7243784000 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 2015 | \n", + "7256490011 | \n", + "7.336435e+09 | \n", + "7349472000 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
66 rows × 11 columns
\n", + "| \n", + " | value | \n", + "
|---|---|
| 1950 | \n", + "2.557629 | \n", + "
| \n", + " | Population Reference Bureau (1973–2015)[6] | \n", + "United Nations Department of Economic and Social Affairs (2015)[7] | \n", + "Maddison (2008)[8] | \n", + "HYDE (2010)[citation needed] | \n", + "Tanton (1994)[9] | \n", + "Biraben (1980)[10] | \n", + "McEvedy & Jones (1978)[11] | \n", + "Thomlinson (1975)[12] | \n", + "Durand (1974)[13] | \n", + "Clark (1967)[14] | \n", + "
|---|---|---|---|---|---|---|---|---|---|---|
| Year | \n", + "\n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " |
| -10000 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "2M[15] | \n", + "NaN | \n", + "NaN | \n", + "4.0 | \n", + "1–10M | \n", + "NaN | \n", + "NaN | \n", + "
| -9000 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "4. | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| -8000 | \n", + "5.0 | \n", + "NaN | \n", + "NaN | \n", + "5. | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "5–10M | \n", + "NaN | \n", + "
| -7000 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "8. | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| -6000 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "11. | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| \n", + " | Population Reference Bureau (1973–2015)[6] | \n", + "United Nations Department of Economic and Social Affairs (2015)[7] | \n", + "Maddison (2008)[8] | \n", + "HYDE (2010)[citation needed] | \n", + "Tanton (1994)[9] | \n", + "Biraben (1980)[10] | \n", + "McEvedy & Jones (1978)[11] | \n", + "Thomlinson (1975)[12] | \n", + "Durand (1974)[13] | \n", + "Clark (1967)[14] | \n", + "
|---|---|---|---|---|---|---|---|---|---|---|
| Year | \n", + "\n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " | \n", + " |
| 1913 | \n", + "NaN | \n", + "NaN | \n", + "1793. | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1920 | \n", + "NaN | \n", + "1860.0 | \n", + "1863. | \n", + "1912. | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "1968. | \n", + "
| 1925 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "2000.0 | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "
| 1930 | \n", + "NaN | \n", + "2070.0 | \n", + "NaN | \n", + "2092. | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "2145. | \n", + "
| 1940 | \n", + "NaN | \n", + "2300.0 | \n", + "2299. | \n", + "2307. | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "NaN | \n", + "2340. | \n", + "
| \n", + " | value | \n", + "
|---|---|
| t0 | \n", + "1950.000000 | \n", + "
| t_end | \n", + "2015.000000 | \n", + "
| p0 | \n", + "2.557629 | \n", + "
| alpha | \n", + "0.025000 | \n", + "
| beta | \n", + "-0.001800 | \n", + "
| \n", + " | United States Census Bureau (2015)[18] | \n", + "Population Reference Bureau (1973-2015)[6] | \n", + "United Nations Department of Economic and Social Affairs (2015)[7] | \n", + "
|---|---|---|---|
| Year | \n", + "\n", + " | \n", + " | \n", + " |
| 2016 | \n", + "7.334772e+09 | \n", + "NaN | \n", + "7.432663e+09 | \n", + "
| 2017 | \n", + "7.412779e+09 | \n", + "NaN | \n", + "NaN | \n", + "
| 2018 | \n", + "7.490428e+09 | \n", + "NaN | \n", + "NaN | \n", + "
| 2019 | \n", + "7.567403e+09 | \n", + "NaN | \n", + "NaN | \n", + "
| 2020 | \n", + "7.643402e+09 | \n", + "NaN | \n", + "7.758157e+09 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| S | \n", + "89 | \n", + "
| I | \n", + "1 | \n", + "
| R | \n", + "0 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| S | \n", + "0.988889 | \n", + "
| I | \n", + "0.011111 | \n", + "
| R | \n", + "0.000000 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| init | \n", + "S 0.988889\n", + "I 0.011111\n", + "R 0.000000\n", + "dtyp... | \n", + "
| t0 | \n", + "0 | \n", + "
| t_end | \n", + "98 | \n", + "
| beta | \n", + "0.333333 | \n", + "
| gamma | \n", + "0.25 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| S | \n", + "0.985226 | \n", + "
| I | \n", + "0.011996 | \n", + "
| R | \n", + "0.002778 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| S | \n", + "0.520453 | \n", + "
| I | \n", + "0.000615 | \n", + "
| R | \n", + "0.478933 | \n", + "
| \n", + " | S | \n", + "I | \n", + "R | \n", + "
|---|---|---|---|
| 0 | \n", + "0.988889 | \n", + "0.011111 | \n", + "0.000000 | \n", + "
| 1 | \n", + "0.985226 | \n", + "0.011996 | \n", + "0.002778 | \n", + "
| 2 | \n", + "0.981287 | \n", + "0.012936 | \n", + "0.005777 | \n", + "
| 3 | \n", + "0.977055 | \n", + "0.013934 | \n", + "0.009011 | \n", + "
| 4 | \n", + "0.972517 | \n", + "0.014988 | \n", + "0.012494 | \n", + "
| \n", + " | S | \n", + "I | \n", + "R | \n", + "
|---|---|---|---|
| 0 | \n", + "0.988889 | \n", + "0.011111 | \n", + "0.000000 | \n", + "
| 1 | \n", + "0.985230 | \n", + "0.011992 | \n", + "0.002778 | \n", + "
| 2 | \n", + "0.981296 | \n", + "0.012929 | \n", + "0.005776 | \n", + "
| 3 | \n", + "0.977071 | \n", + "0.013921 | \n", + "0.009008 | \n", + "
| 4 | \n", + "0.972541 | \n", + "0.014970 | \n", + "0.012488 | \n", + "
| \n", + " | 0.1 | \n", + "0.3 | \n", + "0.5 | \n", + "0.7 | \n", + "
|---|---|---|---|---|
| 0.10 | \n", + "0.084693 | \n", + "0.005444 | \n", + "0.002736 | \n", + "0.001827 | \n", + "
| 0.18 | \n", + "0.708623 | \n", + "0.015914 | \n", + "0.006118 | \n", + "0.003783 | \n", + "
| 0.26 | \n", + "0.900780 | \n", + "0.055380 | \n", + "0.011639 | \n", + "0.006427 | \n", + "
| 0.34 | \n", + "0.956888 | \n", + "0.267864 | \n", + "0.022115 | \n", + "0.010191 | \n", + "
| 0.42 | \n", + "0.977045 | \n", + "0.524563 | \n", + "0.047816 | \n", + "0.015946 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| temp | \n", + "90 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| init | \n", + "temp 90\n", + "dtype: int64 | \n", + "
| volume | \n", + "300 | \n", + "
| r | \n", + "0.01 | \n", + "
| T_env | \n", + "22 | \n", + "
| t0 | \n", + "0 | \n", + "
| t_end | \n", + "30 | \n", + "
| dt | \n", + "1 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| temp | \n", + "89.32 | \n", + "
| \n", + " | temp | \n", + "
|---|---|
| 0 | \n", + "90.000000 | \n", + "
| 1 | \n", + "89.320000 | \n", + "
| 2 | \n", + "88.646800 | \n", + "
| 3 | \n", + "87.980332 | \n", + "
| 4 | \n", + "87.320529 | \n", + "
| 5 | \n", + "86.667323 | \n", + "
| 6 | \n", + "86.020650 | \n", + "
| 7 | \n", + "85.380444 | \n", + "
| 8 | \n", + "84.746639 | \n", + "
| 9 | \n", + "84.119173 | \n", + "
| 10 | \n", + "83.497981 | \n", + "
| 11 | \n", + "82.883001 | \n", + "
| 12 | \n", + "82.274171 | \n", + "
| 13 | \n", + "81.671430 | \n", + "
| 14 | \n", + "81.074715 | \n", + "
| 15 | \n", + "80.483968 | \n", + "
| 16 | \n", + "79.899128 | \n", + "
| 17 | \n", + "79.320137 | \n", + "
| 18 | \n", + "78.746936 | \n", + "
| 19 | \n", + "78.179466 | \n", + "
| 20 | \n", + "77.617672 | \n", + "
| 21 | \n", + "77.061495 | \n", + "
| 22 | \n", + "76.510880 | \n", + "
| 23 | \n", + "75.965771 | \n", + "
| 24 | \n", + "75.426114 | \n", + "
| 25 | \n", + "74.891852 | \n", + "
| 26 | \n", + "74.362934 | \n", + "
| 27 | \n", + "73.839305 | \n", + "
| 28 | \n", + "73.320912 | \n", + "
| 29 | \n", + "72.807702 | \n", + "
| 30 | \n", + "72.299625 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| temp | \n", + "5.17 | \n", + "
| \n", + " | temp | \n", + "
|---|---|
| 0 | \n", + "5.000000 | \n", + "
| 1 | \n", + "5.170000 | \n", + "
| 2 | \n", + "5.338300 | \n", + "
| 3 | \n", + "5.504917 | \n", + "
| 4 | \n", + "5.669868 | \n", + "
| 5 | \n", + "5.833169 | \n", + "
| 6 | \n", + "5.994837 | \n", + "
| 7 | \n", + "6.154889 | \n", + "
| 8 | \n", + "6.313340 | \n", + "
| 9 | \n", + "6.470207 | \n", + "
| 10 | \n", + "6.625505 | \n", + "
| 11 | \n", + "6.779250 | \n", + "
| 12 | \n", + "6.931457 | \n", + "
| 13 | \n", + "7.082143 | \n", + "
| 14 | \n", + "7.231321 | \n", + "
| 15 | \n", + "7.379008 | \n", + "
| \n", + " | temp | \n", + "
|---|---|
| 0 | \n", + "0.000000e+00 | \n", + "
| 1 | \n", + "0.000000e+00 | \n", + "
| 2 | \n", + "0.000000e+00 | \n", + "
| 3 | \n", + "0.000000e+00 | \n", + "
| 4 | \n", + "0.000000e+00 | \n", + "
| 5 | \n", + "0.000000e+00 | \n", + "
| 6 | \n", + "0.000000e+00 | \n", + "
| 7 | \n", + "0.000000e+00 | \n", + "
| 8 | \n", + "0.000000e+00 | \n", + "
| 9 | \n", + "0.000000e+00 | \n", + "
| 10 | \n", + "0.000000e+00 | \n", + "
| 11 | \n", + "0.000000e+00 | \n", + "
| 12 | \n", + "0.000000e+00 | \n", + "
| 13 | \n", + "-1.421085e-14 | \n", + "
| 14 | \n", + "-1.421085e-14 | \n", + "
| 15 | \n", + "0.000000e+00 | \n", + "
| 16 | \n", + "-1.421085e-14 | \n", + "
| 17 | \n", + "-1.421085e-14 | \n", + "
| 18 | \n", + "-1.421085e-14 | \n", + "
| 19 | \n", + "0.000000e+00 | \n", + "
| 20 | \n", + "-1.421085e-14 | \n", + "
| 21 | \n", + "0.000000e+00 | \n", + "
| 22 | \n", + "0.000000e+00 | \n", + "
| 23 | \n", + "0.000000e+00 | \n", + "
| 24 | \n", + "0.000000e+00 | \n", + "
| 25 | \n", + "0.000000e+00 | \n", + "
| 26 | \n", + "-1.421085e-14 | \n", + "
| 27 | \n", + "0.000000e+00 | \n", + "
| 28 | \n", + "-1.421085e-14 | \n", + "
| 29 | \n", + "0.000000e+00 | \n", + "
| 30 | \n", + "0.000000e+00 | \n", + "
| \n", + " | glucose | \n", + "insulin | \n", + "
|---|---|---|
| time | \n", + "\n", + " | \n", + " |
| 0 | \n", + "92 | \n", + "11 | \n", + "
| 2 | \n", + "350 | \n", + "26 | \n", + "
| 4 | \n", + "287 | \n", + "130 | \n", + "
| 6 | \n", + "251 | \n", + "85 | \n", + "
| 8 | \n", + "240 | \n", + "51 | \n", + "
| 10 | \n", + "216 | \n", + "49 | \n", + "
| 12 | \n", + "211 | \n", + "45 | \n", + "
| 14 | \n", + "205 | \n", + "41 | \n", + "
| 16 | \n", + "196 | \n", + "35 | \n", + "
| 19 | \n", + "192 | \n", + "30 | \n", + "
| 22 | \n", + "172 | \n", + "30 | \n", + "
| 27 | \n", + "163 | \n", + "27 | \n", + "
| 32 | \n", + "142 | \n", + "30 | \n", + "
| 42 | \n", + "124 | \n", + "22 | \n", + "
| 52 | \n", + "105 | \n", + "15 | \n", + "
| 62 | \n", + "92 | \n", + "15 | \n", + "
| 72 | \n", + "84 | \n", + "11 | \n", + "
| 82 | \n", + "77 | \n", + "10 | \n", + "
| 92 | \n", + "82 | \n", + "8 | \n", + "
| 102 | \n", + "81 | \n", + "11 | \n", + "
| 122 | \n", + "82 | \n", + "7 | \n", + "
| 142 | \n", + "82 | \n", + "8 | \n", + "
| 162 | \n", + "85 | \n", + "8 | \n", + "
| 182 | \n", + "90 | \n", + "7 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| G | \n", + "287.03 | \n", + "
| X | \n", + "0.00 | \n", + "
| \n", + " | G | \n", + "X | \n", + "
|---|---|---|
| 0.0 | \n", + "290.000000 | \n", + "0.000000 | \n", + "
| 0.5 | \n", + "287.030000 | \n", + "0.000000 | \n", + "
| 1.0 | \n", + "284.104550 | \n", + "-0.000044 | \n", + "
| 1.5 | \n", + "281.229268 | \n", + "-0.000090 | \n", + "
| 2.0 | \n", + "278.403529 | \n", + "-0.000096 | \n", + "
| 2.5 | \n", + "275.620863 | \n", + "-0.000020 | \n", + "
| 3.0 | \n", + "272.869334 | \n", + "0.000179 | \n", + "
| 3.5 | \n", + "270.131906 | \n", + "0.000541 | \n", + "
| 4.0 | \n", + "267.386792 | \n", + "0.001052 | \n", + "
| 4.5 | \n", + "264.615400 | \n", + "0.001636 | \n", + "
| 5.0 | \n", + "261.809705 | \n", + "0.002223 | \n", + "
| 5.5 | \n", + "258.971601 | \n", + "0.002740 | \n", + "
| 6.0 | \n", + "256.112255 | \n", + "0.003160 | \n", + "
| 6.5 | \n", + "253.245932 | \n", + "0.003498 | \n", + "
| 7.0 | \n", + "250.384283 | \n", + "0.003770 | \n", + "
| 7.5 | \n", + "247.536483 | \n", + "0.003992 | \n", + "
| 8.0 | \n", + "244.709378 | \n", + "0.004176 | \n", + "
| 8.5 | \n", + "241.907820 | \n", + "0.004334 | \n", + "
| 9.0 | \n", + "239.134995 | \n", + "0.004478 | \n", + "
| 9.5 | \n", + "236.392528 | \n", + "0.004620 | \n", + "
| 10.0 | \n", + "233.680592 | \n", + "0.004763 | \n", + "
| 10.5 | \n", + "230.998839 | \n", + "0.004906 | \n", + "
| 11.0 | \n", + "228.347255 | \n", + "0.005044 | \n", + "
| 11.5 | \n", + "225.726142 | \n", + "0.005176 | \n", + "
| 12.0 | \n", + "223.136086 | \n", + "0.005300 | \n", + "
| 12.5 | \n", + "220.577749 | \n", + "0.005417 | \n", + "
| 13.0 | \n", + "218.051663 | \n", + "0.005528 | \n", + "
| 13.5 | \n", + "215.558239 | \n", + "0.005633 | \n", + "
| 14.0 | \n", + "213.097773 | \n", + "0.005732 | \n", + "
| 14.5 | \n", + "210.670550 | \n", + "0.005825 | \n", + "
| ... | \n", + "... | \n", + "... | \n", + "
| 167.5 | \n", + "90.548441 | \n", + "-0.000454 | \n", + "
| 168.0 | \n", + "90.590749 | \n", + "-0.000465 | \n", + "
| 168.5 | \n", + "90.632952 | \n", + "-0.000477 | \n", + "
| 169.0 | \n", + "90.675052 | \n", + "-0.000488 | \n", + "
| 169.5 | \n", + "90.717051 | \n", + "-0.000499 | \n", + "
| 170.0 | \n", + "90.758950 | \n", + "-0.000511 | \n", + "
| 170.5 | \n", + "90.800752 | \n", + "-0.000522 | \n", + "
| 171.0 | \n", + "90.842458 | \n", + "-0.000534 | \n", + "
| 171.5 | \n", + "90.884070 | \n", + "-0.000545 | \n", + "
| 172.0 | \n", + "90.925591 | \n", + "-0.000557 | \n", + "
| 172.5 | \n", + "90.967023 | \n", + "-0.000568 | \n", + "
| 173.0 | \n", + "91.008367 | \n", + "-0.000580 | \n", + "
| 173.5 | \n", + "91.049626 | \n", + "-0.000591 | \n", + "
| 174.0 | \n", + "91.090802 | \n", + "-0.000603 | \n", + "
| 174.5 | \n", + "91.131899 | \n", + "-0.000614 | \n", + "
| 175.0 | \n", + "91.172917 | \n", + "-0.000626 | \n", + "
| 175.5 | \n", + "91.213861 | \n", + "-0.000638 | \n", + "
| 176.0 | \n", + "91.254731 | \n", + "-0.000649 | \n", + "
| 176.5 | \n", + "91.295531 | \n", + "-0.000661 | \n", + "
| 177.0 | \n", + "91.336264 | \n", + "-0.000673 | \n", + "
| 177.5 | \n", + "91.376932 | \n", + "-0.000684 | \n", + "
| 178.0 | \n", + "91.417539 | \n", + "-0.000696 | \n", + "
| 178.5 | \n", + "91.458085 | \n", + "-0.000708 | \n", + "
| 179.0 | \n", + "91.498576 | \n", + "-0.000719 | \n", + "
| 179.5 | \n", + "91.539013 | \n", + "-0.000731 | \n", + "
| 180.0 | \n", + "91.579400 | \n", + "-0.000743 | \n", + "
| 180.5 | \n", + "91.619739 | \n", + "-0.000755 | \n", + "
| 181.0 | \n", + "91.660033 | \n", + "-0.000767 | \n", + "
| 181.5 | \n", + "91.700286 | \n", + "-0.000779 | \n", + "
| 182.0 | \n", + "91.740501 | \n", + "-0.000791 | \n", + "
365 rows × 2 columns
\n", + "| \n", + " | G | \n", + "X | \n", + "
|---|---|---|
| time | \n", + "\n", + " | \n", + " |
| 0 | \n", + "290.000000 | \n", + "0.000000 | \n", + "
| 2 | \n", + "278.503218 | \n", + "-0.000072 | \n", + "
| 4 | \n", + "267.407044 | \n", + "0.001329 | \n", + "
| 6 | \n", + "255.962566 | \n", + "0.003315 | \n", + "
| 8 | \n", + "244.552487 | \n", + "0.004232 | \n", + "
| 10 | \n", + "233.570125 | \n", + "0.004808 | \n", + "
| 12 | \n", + "223.075469 | \n", + "0.005331 | \n", + "
| 14 | \n", + "213.094204 | \n", + "0.005750 | \n", + "
| 16 | \n", + "203.651201 | \n", + "0.006054 | \n", + "
| 19 | \n", + "190.540814 | \n", + "0.006312 | \n", + "
| 22 | \n", + "178.671483 | \n", + "0.006494 | \n", + "
| 27 | \n", + "161.389180 | \n", + "0.006705 | \n", + "
| 32 | \n", + "146.879018 | \n", + "0.006894 | \n", + "
| 42 | \n", + "124.440473 | \n", + "0.007088 | \n", + "
| 52 | \n", + "109.230932 | \n", + "0.006417 | \n", + "
| 62 | \n", + "99.405510 | \n", + "0.005610 | \n", + "
| 72 | \n", + "93.178906 | \n", + "0.004775 | \n", + "
| 82 | \n", + "89.501328 | \n", + "0.003855 | \n", + "
| 92 | \n", + "87.562395 | \n", + "0.002954 | \n", + "
| 102 | \n", + "86.786289 | \n", + "0.002279 | \n", + "
| 122 | \n", + "86.832321 | \n", + "0.001258 | \n", + "
| 142 | \n", + "88.339444 | \n", + "0.000219 | \n", + "
| 162 | \n", + "90.119786 | \n", + "-0.000323 | \n", + "
| 182 | \n", + "91.759885 | \n", + "-0.000789 | \n", + "
| \n", + " | G | \n", + "X | \n", + "
|---|---|---|
| 2.0 | \n", + "-0.035794 | \n", + "34.378714 | \n", + "
| 4.0 | \n", + "-0.007574 | \n", + "-20.903093 | \n", + "
| 6.0 | \n", + "0.058481 | \n", + "-4.666342 | \n", + "
| 8.0 | \n", + "0.064154 | \n", + "-1.328219 | \n", + "
| 10.0 | \n", + "0.047295 | \n", + "-0.934767 | \n", + "
| 12.0 | \n", + "0.027173 | \n", + "-0.582564 | \n", + "
| 14.0 | \n", + "0.001675 | \n", + "-0.312059 | \n", + "
| 16.0 | \n", + "-0.028931 | \n", + "-0.016229 | \n", + "
| 19.0 | \n", + "-0.082709 | \n", + "0.204179 | \n", + "
| 22.0 | \n", + "-0.136334 | \n", + "0.202042 | \n", + "
| 27.0 | \n", + "-0.218693 | \n", + "0.299434 | \n", + "
| 32.0 | \n", + "-0.286757 | \n", + "0.172747 | \n", + "
| 42.0 | \n", + "-0.384822 | \n", + "0.409781 | \n", + "
| 52.0 | \n", + "-0.455501 | \n", + "0.567695 | \n", + "
| 62.0 | \n", + "-0.467895 | \n", + "0.466529 | \n", + "
| 72.0 | \n", + "-0.440646 | \n", + "0.559205 | \n", + "
| 82.0 | \n", + "-0.394390 | \n", + "0.516638 | \n", + "
| 92.0 | \n", + "-0.339047 | \n", + "0.565098 | \n", + "
| 102.0 | \n", + "-0.275395 | \n", + "0.164876 | \n", + "
| 122.0 | \n", + "-0.170525 | \n", + "0.545027 | \n", + "
| 142.0 | \n", + "-0.102344 | \n", + "-0.885106 | \n", + "
| 162.0 | \n", + "-0.048847 | \n", + "1.056507 | \n", + "
| 182.0 | \n", + "-0.021124 | \n", + "0.266167 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| S_G | \n", + "0.027514 | \n", + "
| S_I | \n", + "0.000849 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| init | \n", + "I 360\n", + "dtype: int64 | \n", + "
| I0 | \n", + "10 | \n", + "
| gamma | \n", + "0.004 | \n", + "
| k | \n", + "0.25 | \n", + "
| G_T | \n", + "80 | \n", + "
| G | \n", + "<scipy.interpolate.interpolate.interp1d object... | \n", + "
| ts | \n", + "Int64Index([ 0, 2, 4, 6, 8, 10, 12,... | \n", + "
| \n", + " | I | \n", + "
|---|---|
| time | \n", + "\n", + " |
| 0 | \n", + "360.000000 | \n", + "
| 2 | \n", + "219.652719 | \n", + "
| 4 | \n", + "137.739127 | \n", + "
| 6 | \n", + "89.505898 | \n", + "
| 8 | \n", + "61.636133 | \n", + "
| 10 | \n", + "45.765330 | \n", + "
| 12 | \n", + "37.053121 | \n", + "
| 14 | \n", + "32.989732 | \n", + "
| 16 | \n", + "31.406589 | \n", + "
| 19 | \n", + "31.811071 | \n", + "
| 22 | \n", + "32.581632 | \n", + "
| 27 | \n", + "34.014287 | \n", + "
| 32 | \n", + "33.758216 | \n", + "
| 42 | \n", + "30.932902 | \n", + "
| 52 | \n", + "24.554216 | \n", + "
| 62 | \n", + "15.767804 | \n", + "
| 72 | \n", + "7.782053 | \n", + "
| 82 | \n", + "-0.382598 | \n", + "
| 92 | \n", + "0.598298 | \n", + "
| 102 | \n", + "1.943431 | \n", + "
| 122 | \n", + "3.415756 | \n", + "
| 142 | \n", + "4.413573 | \n", + "
| 162 | \n", + "11.170251 | \n", + "
| 182 | \n", + "25.702699 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| y | \n", + "381 meter | \n", + "
| v | \n", + "0.0 meter / second | \n", + "
| \n", + " | y | \n", + "v | \n", + "
|---|---|---|
| 0.000 | \n", + "381.000000 | \n", + "0.0000 | \n", + "
| 0.191 | \n", + "380.821243 | \n", + "-1.8718 | \n", + "
| 0.382 | \n", + "380.284972 | \n", + "-3.7436 | \n", + "
| 0.573 | \n", + "379.391188 | \n", + "-5.6154 | \n", + "
| 0.764 | \n", + "378.139890 | \n", + "-7.4872 | \n", + "
| \n", + " | y | \n", + "v | \n", + "
|---|---|---|
| 18.336 | \n", + "-1266.423590 | \n", + "-179.6928 | \n", + "
| 18.527 | \n", + "-1300.923672 | \n", + "-181.5646 | \n", + "
| 18.718 | \n", + "-1335.781268 | \n", + "-183.4364 | \n", + "
| 18.909 | \n", + "-1370.996377 | \n", + "-185.3082 | \n", + "
| 19.100 | \n", + "-1406.569000 | \n", + "-187.1800 | \n", + "
| \n", + " | y | \n", + "v | \n", + "
|---|---|---|
| 0.0 | \n", + "381.0 | \n", + "8.0 | \n", + "
| 1.0 | \n", + "384.1 | \n", + "-1.8 | \n", + "
| 2.0 | \n", + "377.4 | \n", + "-11.6 | \n", + "
| 3.0 | \n", + "360.9 | \n", + "-21.4 | \n", + "
| 4.0 | \n", + "334.6 | \n", + "-31.2 | \n", + "
| 5.0 | \n", + "298.5 | \n", + "-41.0 | \n", + "
| 6.0 | \n", + "252.6 | \n", + "-50.8 | \n", + "
| 7.0 | \n", + "196.9 | \n", + "-60.6 | \n", + "
| 8.0 | \n", + "131.4 | \n", + "-70.4 | \n", + "
| 9.0 | \n", + "56.1 | \n", + "-80.2 | \n", + "
| 10.0 | \n", + "-29.0 | \n", + "-90.0 | \n", + "
| \n", + " | value | \n", + "
|---|---|
| init | \n", + "y 381 meter\n", + "v 0.0 meter / secon... | \n", + "
| g | \n", + "9.8 meter / second ** 2 | \n", + "
| mass | \n", + "0.0025 kilogram | \n", + "
| rho | \n", + "1.2 kilogram / meter ** 3 | \n", + "
| C_d | \n", + "0.4445009981135434 dimensionless | \n", + "
| area | \n", + "0.0002835287369864788 meter ** 2 | \n", + "
| ts | \n", + "[0.0 second, 0.3 second, 0.6 second, 0.8999999... | \n", + "
| \n", + " | value | \n", + "
|---|---|
| init | \n", + "x 0 meter\n", + "y ... | \n", + "
| g | \n", + "9.8 meter / second ** 2 | \n", + "
| mass | \n", + "0.145 kilogram | \n", + "
| area | \n", + "0.004185386812745002 meter ** 2 | \n", + "
| rho | \n", + "1.2 kilogram / meter ** 3 | \n", + "
| C_d | \n", + "0.3 | \n", + "
| ts | \n", + "[0.0 second, 0.051 second, 0.102 second, 0.153... | \n", + "
| \n", + " | x | \n", + "y | \n", + "vx | \n", + "vy | \n", + "
|---|---|---|---|---|
| 0.000 | \n", + "0.000000 | \n", + "1.000000 | \n", + "10.352762 | \n", + "38.637033 | \n", + "
| 0.051 | \n", + "0.525223 | \n", + "2.947460 | \n", + "10.244833 | \n", + "37.737037 | \n", + "
| 0.102 | \n", + "1.045033 | \n", + "4.849486 | \n", + "10.140414 | \n", + "36.855149 | \n", + "
| 0.153 | \n", + "1.559603 | \n", + "6.706980 | \n", + "10.039371 | \n", + "35.990594 | \n", + "
| 0.204 | \n", + "2.069104 | \n", + "8.520808 | \n", + "9.941577 | \n", + "35.142638 | \n", + "
| \n", + " | x | \n", + "y | \n", + "vx | \n", + "vy | \n", + "
|---|---|---|---|---|
| 4.896 | \n", + "39.033848 | \n", + "42.325221 | \n", + "6.726804 | \n", + "-15.869290 | \n", + "
| 4.947 | \n", + "39.376127 | \n", + "41.505021 | \n", + "6.695814 | \n", + "-16.294823 | \n", + "
| 4.998 | \n", + "39.716812 | \n", + "40.663211 | \n", + "6.664295 | \n", + "-16.716737 | \n", + "
| 5.049 | \n", + "40.055876 | \n", + "39.799976 | \n", + "6.632259 | \n", + "-17.134973 | \n", + "
| 5.100 | \n", + "40.393294 | \n", + "38.915507 | \n", + "6.599718 | \n", + "-17.549471 | \n", + "