(2pi G)/log(b/a) for a=coreradius, b=shellradius