Numerical and experimental verification were carried out, and the results showed that method (2-2), i.e. the method with a matrix order of 3 and matrix modification, exhibited the highest accuracy. However, numerical verification suggested that some error will occur when method (2-2) is applied for a target material having low hardening exponent nu and very low offset of strain ε0u.