Add a check for non-increasing V
@@ -1138,6 +1138,9 @@ def solveConsRiskyContrib(solution_next,ShockDstn,IncomeDstn,RiskyDstn,
# Evaluate realizations of value next period conditional on shocks, assets
# and shares
+ if not vFuncAdj_next(0.02,0.6) >= vFuncAdj_next(0.02,0.05):
+ print('VFUNC NOT INCREASING!!!')
+
vAdj_next = vFuncAdj_next(mNrm_next, nNrm_next)
if AdjustPrb < 1.:
vFxd_next = vFuncFxd_next(mNrm_next, nNrm_next, Share_next)
