This code and its documentation is a result of solving the programming example for the lecture Numerical Methods in Fluid Mechanics and Heat Transfer (LV-Nr.: 321.023 VO, WS 2012/13) @ Graz university of Technology.