A n-body problem solver with implicit euler scheme