This code solves the Aiyagari Model based on lecture code of Christian Bayer from the University of Bonn. To solve the model, changes were made in the main file, as well as in the Young function such that only one stationary distribution is found in the eigenvalue decomposition.