Fixed several bug:

* Formula for the velocity
* Sun's grav par introduced to compute Jupiter orbital velocity
* Correction for baricentric coordinate
Check to see if the perijove is inside Jupiter
Added several data sources and improved naming in the loop
little fixes
Complete rewrite because of wrong assumptions
clear all -> clearvars (so that breackpoints do not disappear)
Little comment fix