added test cases
added util.h
added stdafx.cpp
eliminated  several errors: index of f(1,0) and such, and also inverse of matrix A