Simulation of general relativistic corrections in long term numerical integrations of planetary orbits