\(\begin{align} &\Delta E = \sqrt{ \left(\dfrac{\Delta L'}{K_LS_L}\right)^2+ \left(\dfrac{\Delta C'}{K_CS_C}\right)^2+ \left(\dfrac{\Delta H'}{K_HS_H}\right)^2+ R_T\left(\dfrac{\Delta C'}{K_CS_C}\right) \left(\dfrac{\Delta H'}{K_HS_H}\right)}\\ &\bar{L}' = 0.5(L_1+L_2)\\ &C_1 = \sqrt{a_1^2+b_1^2}\\ &C_2 = \sqrt{a_2^2+b_2^2}\\ &\bar{C} = 0.5(C_1+C_2)\\ &G = 0.5\left(1-\sqrt{\dfrac{\bar{C}^7}{\bar{C}^7+25^7}}\right)\\ &a_1' = a_1(1+G)\\ &a_2' = a_2(1+G)\\ &C_1' = \sqrt{a_1'^2+b_1^2}\\ &C_2' = \sqrt{a_2'^2+b_2^2}\\ &\bar{C}' = 0.5\left(C_1'+C_2'\right)\\ &h_1' = \begin{cases} \text{arctan}\left(\dfrac{b_1}{a_1'}\right)&\text{if }\text{arctan}\left(\dfrac{b_1}{a_1'}\right) \ge 0\\ \text{arctan}\left(\dfrac{b_1}{a_1'}\right)+360&\text{else} \end{cases}\\ &h_2' = \begin{cases} \text{arctan}\left(\dfrac{b_2}{a_2'}\right)&\text{if }\text{arctan}\left(\dfrac{b_2}{a_2'}\right) \ge 0\\ \text{arctan}\left(\dfrac{b_2}{a_2'}\right)+360&\text{else} \end{cases}\\ &\bar{H}' = \begin{cases} 0.5(h_1'+h_2'+360)&\text{if }|h_1'-h_2'| \gt 180\\ 0.5(h_1'+h_2')&\text{else} \end{cases}\\ &T = 1-0.17\cos\left(\bar{H}'-30\right)+0.24\cos\left(2\bar{H}'\right)+0.32\cos\left(3\bar{H}'+6\right)-0.20\cos\left(4\bar{H}'-63\right)\\ &\Delta h' = \begin{cases} h_2'-h_1'&\text{if }|h_2'-h_1'| \le 180\\ h_2'-h_1'+360&\text{else if }h_2' \le h_1'\\ h_2'-h_1'-360&\text{else}\\ \end{cases}\\ &\Delta L = L_2-L_1\\ &\Delta C = C_2-C_1\\ &\Delta H = 2\sqrt{C_1'C_2'}\sin(0.5\Delta h')\\ &S_L = 1+\dfrac{0.015\left(\bar{L}'-50\right)^2}{\sqrt{20+\left(\bar{L}'-50\right)^2}}\\ &S_C = 1+0.045\bar{C}'\\ &S_H = 1+0.015\bar{C}'T\\ &\Delta\theta = 30\exp\left(-\left(\dfrac{\bar{H}'-275}{25}\right)^2\right)\\ &R_C = 2\sqrt{\dfrac{\bar{C}'^7}{\bar{C}'^7+25^7}}\\ &R_T = -R_C\sin(2\Delta\theta)\\ &K_L = 1\text{ default}\\ &K_C = 1\text{ default}\\ &K_H = 1\text{ default} \end{align}\)