Try implementing value function iteration using Mata (Stata's matrix programming language), based on some problem/notes you have (or can find). For example, the following notes from Jerome Adda, use very basic Matlab and comment on the general algorithm: