Thanks. Found it in my archives. (SAE paper 650507). In VB in looks like this:
Function spindtAFR(O2percent, COpercent, CO2percent, HCppmC, HCratio)
O2 = O2percent
CO2 = CO2percent
CO = COpercent
HC = HCppmC
q = O2 / CO2
r = CO / CO2
Fb = (CO2 + CO) / (CO2 + CO + (0.0001 * HC))
Fc = 12.01 / (12.01 + 1.008 * HCratio)
spindtAFR = Fb * ((11.492 * Fc * ((1 + 0.5 * r + q) / _
(1 + r))) + ((120 * (1 - Fc)) / (3.5 + r)))
End Function