Time of the seasons

[Gil]

My reason of my request was because
I try to reproduce the Meeus' algorithm to get the June solstice for the year 1962,with correct answer in Terrestrial Time or (Terrestrial) Dynamical Time within 2-4 seconds.

I work out author's example 27.b (pages 180-181),
using the tables given in his book at pages 418-421.

I get the following results on my HP50G:


:JDE0: 2437837.38588 <Jean Meeus... 35889>
tau = (JDEO-2451545)/365230
:SL0: 176069483.423
:SL1: 628332128985.
:SL2: 47714.5774468
:SL3: -192.115254877
:SL4: -109.4139133
:SL5: -.999998731728
'(SL0+SL1*tau+SL2*tau^2+SL3*tau^3+SL4*tau^4+SL5*tau^5)/100000000'
:L[rad]: -234.048595755 :
L[d.d]: 270.0032623
:R: 1.01630134777

with these L and R values diverging slightly
from the ones given in the book on the top of page 181.