julien

from HARK.ConsumptionSaving.ConsIndShockModel import IndShockConsumerType
import HARK.ConsumptionSaving.ConsumerParameters as param
import copy

init_lifecycle = copy.copy(param.init_idiosyncratic_shocks)
init_lifecycle['PermGroFac'] = [1.01, 1.00]
init_lifecycle['PermShkStd'] = [0.1, 0.0]
init_lifecycle['TranShkStd'] = [0.3,0.0]
init_lifecycle['LivPrb']     = [0.99,0.9]
init_lifecycle['T_cycle']    = 2
init_lifecycle['T_retire']   = 0
init_lifecycle['T_age'] = 2 # Make sure that old people die at terminal age and don't turn into newborns!

pcct = IndShockConsumerType(**init_lifecycle)

pcct.solve()
pcct.T_sim = 30
pcct.AgentCount = 1
pcct.initializeSim()
pcct.simulate()
pcct.track_vars=['mNrmNow', 'aNrmNow', 'cNrmNow', 't_age', 'TranShkNow']
pcct.initializeSim()
pcct.simulate()