A recursive formulation for the tesseral disturbing function in equinoctial variables