The following were calculated with various incarnations of my conjunction program.
The latest runs are being tested on a Mac in 2015.
The SEPARATION column is geocentric for the conjunction in ecliptic longitude. The distance of the planetary centers is shown. The SIN RM column helps to convert this to the closest approach distance. Then the parallax shift is considered for the best viewing location on Earth to yield a maximum overlap. If this overlap has a positive value an occultation is flagged.
I'm presently assessing the accuracy of this run. For example, it's possible the positions of Mercury (perhaps more than most) can drift a few arcseconds over the 2000 year integration. It's interesting that the upcoming MERCURY NEPTUNE occultation in 2067 is a rare event with a recurrence only around every 2000 years. I'm experimenting with some even longer runs and I'm starting to find some Jupiter-Saturn occultations in the distant past with the most recent around 8000 years ago.
These are geocentric conjunctions (run on a UNIX system in 2003).
These are geocentric conjunctions in celestial longitude (VAX version run around 1994).
These are three-planet conjunctions when the planets fit in the smallest circle.
Additional information on these calculations is found in "Mutual Occultations of Planets: 1557-2230", Sky and Telescope, March 1979.
There was a remarkable coincidence that the printed version of the Sky and Telescope article appeared on the same page as an article about Galileo and his lunar drawings. My article also suggested examination of historical observations to see if Uranus or Neptune could have been mistaken for a moon of Jupiter and thus being seen before it was discovered. Stillman Drake followed up on that notion using Galileo's observations from 1613. It turns out that Neptune actually shows up as a moving "star" on Galileo's drawings so Galileo himself unwittingly observed Neptune right around the time of the invention of the telescope. This was published in Scientific American in 1980. My role in this was recognized in receiving the R.R. Newton award by the International Journal of the Scientific History in 2009.
Some additional occultations were found by Larry Bogan on his website.