Project 3 of computational physics