The main version of VSOP87 defines elliptic variables for the planets as a function of time.
The VSOP data provides the following elliptic variables:
A the semi-major axis in Astronomic Units (AU where 1AU = 149597870.7km)
L the mean longitude
k = e cos π where e is the eccentricity and π is the longitude of the perihelion
h = e sin π
q = sin i/2 cos Ω where i is inclination and Ω is the longitude of the ascending node
p = sin i/2 sin Ω
The calculations are based on the time in Julian centuries after the J2000 epoch 2000-01-01 12:00:00.
T = (JDE - 2451545)/365250
Meysam Mahooti (2020). Sun Position (VSOP87) (https://www.mathworks.com/matlabcentral/fileexchange/77248-sun-position-vsop87), MATLAB Central File Exchange. Retrieved .