You should rewrite the differential equation as a canonical state-space model, as dx/dt = f(x). And then, build a ANFIS model for the last row of the f(x), using function approximation capability of ANFIS structure.
For this purpose, you can use the following codes: