celestial objects and orbital paths

I need a physicist to create an algorithm, that takes 2 celestial objects (one is stationary for simplicity) as parameters with mass, global position, velocity and acceleration vectors, and return a path (of max 20-50 points for drawing) of the second celestial body around the stationary (huge difference in mass!). For simplicity one object is stationary and of course the SOI is a predefined radius! Looking for future cooperation too! I need a lot of performance, so the algorithm must be fast! (should work based on mathematical formulas, no loop iterations, I have tried it!!) It's just a simple algorithm for the one who knows what he is doing! So no big bids!!! Thanks in advance!