How to change CSpice initial velocity from an Asteroid SPK file?

[forofo]

TL;DR Summary

Change spice spk file velocity to determinate the position of an asteroid

Hello, I've made a SPK file for asteroid 7482 (1994 PC1) with Horizon. I wan't to change the initial velocity with cspice, because I want to know where it will be in a future time at the speed changed. (now I'm using Newtons calc but is slow an error increases with time). This is for calculating NEO's asteroid position desviation.

 

[forofo]

I found the answer: In this link https://github.com/fding/nbody (Newtonian mechanics and a four-order symplectic Candy-Rozmus integration), there is c++ code and phyton. I've modify the c++, so now I calc the inital positions and speeds with cspice and pass then to Newton candy with other asteroid speed. Then I run both to the future to see the asteroid distance to any planet. Just take care about units conversion, cspice frame(ECLIPJ2000) and obverser naif Solar System(0), Used de440.bsp and 2007482(1994 PC1).bsp

 

[sysprog]

You might want have APK Editor in your tool kit.