test function added

old code from before kernel matrix and forcing vectors were moved to
separate function and before Poisson's ratio and density were allowed
to be different for both sides of the crack