@ecterrab 

Hi,

Thank you for your answer.

actually, i have recently read a paper in this way but i didn't understand the method and algorithm. The authors represented an algorithm for computing constraint structure via Grobner bases using Maple in this paper: arxiv:math/9909113v2

But the authors didn't gave Maple code and i didn't understand the algorithm too.

would you please help me in this way?

@Mariusz Iwaniuk 

Thank you very much. it helps a lot :)

@Mariusz Iwaniuk 

Thank You very much, my friend.

It is very helpful and I appreciate for this. But it gives us an answer in integral form. Is it possible to convert the answer to an Algebraic term or at least a power series?