Integration by parts: The algorithm to calculate β-functions in 4 loops