为什么C122错误？
y112 = 0.92 + 0.0013 * (T - 270)
y122 = 0.92 + 0.0013 * (T - 270)
C111 = c0H0 + c1H1 * T + c2H0 * T ^ 2 + Hch * (c0H1 + c1H1 * T + c2H1 * T ^ 2) + (Hch ^ 2) * (c0H2 + c1H2 * T + c2H2 * T ^ 2)
C222 = c0222 + c1222 * T + c2222 * T ^ 2
C333 = c0333 + c1333 * T + c2333 * T ^ 2
C444 = c0444 + c1444 * T + c2444 * T ^ 2
C115 = c0115 + c1115 * T + c2115 * T ^ 2
C223 = c0223 + c1223 * T + c2223 * T ^ 2
C233 = c0233 + c1233 * T + c2233 * T ^ 2
C112 = y112 * (C111 * C111 * C222) ^ (1 / 3)
C113 = y113 * (C111 * C111 * C333) ^ (1 / 3)
C114 = y114 * (C111 * C111 * C444) ^ (1 / 3)

C122 = y122 * (C111 * C222 * C222) ^ (1 / 3)

负数的三分之一次方算不了吗

C112这么改就可以了，C113，C114也类似

C112 = y112 * sgn(C111 * C111 * C222)*abs(C111 * C111 * C222) ^ (1 / 3)