Added http://www.plasmaphysics.org.uk/programs/coll3d_cpp.htm
Added coll3d.c
coll3d doesn't seem to work the way my interact function works
Initial speed can be both positive or negative