def __init__(self, A=2.87, k0=3.5, delta=0.08, s=0.1, n=0.015, alpha=0.36, t0=1956, tmax=2060):
    self._A = A
    self._k0 = k0
    self._k = k0
    self._delta = delta
    self._s = s
    self._n = n
    self._alpha = alpha
    self._t0 = t0
    self._tmax = tmax
    self._t = range(t0, (tmax + 1))
    self._y = np.zeros(len(self._t))
    self._y[0] = (self._A * (self._k0 ** self._alpha))
    self._time_passed = 0
def __init__(self, A=2.87, k0=3.5, delta=0.08, s=0.1, n=0.015, alpha=0.36, t0=1956, tmax=2060):
    self._A = A
    self._k0 = k0
    self._k = k0
    self._delta = delta
    self._s = s
    self._n = n
    self._alpha = alpha
    self._t0 = t0
    self._tmax = tmax
    self._t = range(t0, (tmax + 1))
    self._y = np.zeros(len(self._t))
    self._y[0] = (self._A * (self._k0 ** self._alpha))
    self._time_passed = 0