correct bug when SEP_FXN==0
add test function

this gives results from before changes to allow different densities,
Poisson's ratio and Young's modulus happened;
and also before the kernel matrix/forcing vectors were moved to a
separate function