b y l r = c D r ∗ ( 1 + A G T E C r ) t − 1 ∗ K A G r C D A L F r , s = 1 ∗ L A B S r , s = 1 ( 1 − C D A L F r , s = 1 ) ∗ S A T K r {\displaystyle byl_{r}=cD_{r}*{(1+AGTEC_{r})}_{t-1}*KAG_{r}^{CDALF_{r,s=1}}*LABS_{r,s=1}^{{(1-CDALF}_{r,s=1)}}*SATK_{r}} <math>CD_{r}=YL_{r,t=1}/{KAG_{r,t=1}^CDALF_{r,s=1}*LABS_{r,s=1,t=1}^{{(1-CDALF_{r,s=1})}