这是北航数值分析的第一道题目，要求求一501*501三角形矩阵的特征值，需要用到幂法、反幂法、LU分解解方程组等各个子程序

LU矩阵分解代码，直接输入系数矩阵，输出解，A=LU，LUx=b

研究生上有限元数值分析课时，老师布置的相关作业，此时刚接触Python，也没有多少编程经验，编程过程对我来说并不简单，抄抄改改，还好总算能够运行，达到自己的目标值。内容包含：LU分解解线性方程组，Jacobi迭代，implicit,explicit,总刚度矩阵的组装，以及简单的4单元9节点有限元分析（得到各个节点的位移和力）。