holonomic function