Have you posted the formulas that you use for that somewhere? I'd be interested to see how the process works.
I did, but I don't remember which thread it was...
But the idea is simple: I can easily calculate the speed of any object at the periapsis and at the apoapsis of its trajectory. That also applies to any point of a circular orbit of course.
For example, in the case of the transfer Earth to Saturn, I consider a transfer orbit, which periapsis is at Earth level, and apoapsis at Saturn level. I can calculate the ship's speed at Earth level, and the Earth's speed, which gives me the speed difference between Earth and the ship (I consider the velocities are aligned). I also have to take into account the speed loss due to Earth's gravity, but the formula for that is simple. This finally gives me the speed I should reach at Earth's level. By considering I am already at 1655 m/s (the speed required to be in LEO), that gives me the delta-V to provide from LEO.
When I want to calculate the delta-V to provide for the insertion burn on Saturn (let's assume I aim for Low Saturn Orbit), that's the same logic. I can calculate the speed difference between Saturn and the ship, and the speed gained from Saturn's gravity. That gives me the speed at LSO level. I can calculate the speed required to be in LSO, so my arrival speed minus LSO speed is the delta-V I'll have to provide.
I already have a file to make all those calculations, so it's relatively fast. The most tedious part is extracting the necessary planet data from the files for the calculations.