**Appendix A**

function *z* = hanshu(*x*) syms *A*1 *b*1 *A*2 *b*2 *y*1 *y*2 *x*1 *x*2 *r c*; syms *x*1\_1 *x*1\_2 *x*1\_3 *x*1\_4 *x*1\_5 *x*1\_6 *x*1\_7 *x*1\_8 *x*1\_9 *x*1\_10 *x*1\_11; syms *x*2\_1 *x*2\_2 *x*2\_3 *x*2\_4 *x*2\_5 *x*2\_6 *x*2\_7 *x*2\_8 *x*2\_9 *x*2\_10 *x*2\_11; *A*1 = *x*(1) *b*1 = *x*(2) *A*2 = *x*(3) *b*2 = *x*(4) *y*1 = [3.76 14.27 19.61 30.25 33.92 42.75 54.32 62.31 70.06 78.77 88.84]; '*t*<sup>10</sup> of anthracite *y*2 = [3.63 8.14 11.33 12.55 15.68 21.52 26.69 30.95 34.84 37.56 42.72]; '*t*<sup>10</sup> of bituminous coal *x*1 = −log(1 − *y*1./*A*1)/*b*1; 'specific energy of anthracite *x*2 = −log(1 − *y*2./*A*2)/*b*2; 'specific energy of bituminous coal *r* = 0.25. × *x*1 + 0.75. × *x*2; *c* = [0.21 0.41 0.61 0.81 1.00 1.37 1.71 2.19 2.66 3.04 3.74]; *z* = abs(abs(*c*(1,1) − *r*(1,1))./max(*c*(1,1),*r*(1,1)) + abs(*c*(1,2) − (1,2))./max(*c*(1,2),*r*(1,2)) +

abs( (1,3) − *r*(1,3))./max(*c*(1,3),*r*(1,3)) + abs(*c*(1,4) − *r*(1,4))./max(*c*(1,4),*r*(1,4)) + abs(*c*(1,5) − *r*(1,5))./max(*c*(1,5),*r*(1,5)) + abs(*c*(1,6) − *r*(1,6))./max(*c*(1,6),*r*(1,6)) + abs(*c*(1,7) − *r*(1,7))./max(*c*(1,7),*r*(1,7)) + abs(*c*(1,8) − *r*(1,8))./max(*c*(1,8),*r*(1,8)) + abs(*c*(1,9) − *r*(1,9))./max(*c*(1,9),*r*(1,9)) + abs(*c*(1,10) − *r*(1,10))./max(*c*(1,10),*r*(1,10)))
