Improving a maximum horizontal gradient algorithm to determine geological body boundaries and fault systems based on gravity data