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 /...