Project 4 in FYS3150/4150 Computational Physics, fall of 2013