Possibly they are entirely working, but they are hard to debug because doing the numerical derivatives correctly is very tricky, especially shifting the mole fractions while keeping molar density constant. A new method might be needed for testing use only.