Shosei Yu
Faculty of Economics, University of Tokyo
From Miranda and Fackler, Applied Computational Economics and Finance, 2002, Section 7.6.5
Julia translation of the Python version
using QuantEcon using Plots plotlyjs()