以下程序求解主应力大小及主方向，以矩阵形式输入点的应力分量,可以以上三角形式输入-The following program to solve the principal stress and the principal direction, the stress components of the matrix form of the input points, can be more than triangular form input