“…Then, the integral constants of each region can be obtained through matrix operation. In addition, the detailed derivation of l 2,3k , s 2,3k , γ 2,3k , f 2,3k , g 2,3k , l 3k,4ik , s 3k,4ik , f 3k,4ik , g 3k,4ik , l 4ik,4jk1 , s 4ik,4jk1 and the method of matrix operation can be seen in [26].…”