%该程序计算平面四节点四边形等参元的刚度矩阵%输入弹性模量E、泊松比NU和厚度h%输入四个节点1、2、3、4的坐标x1,x2,x3,x3,y1,y2,y3,y4%输入平面问题性质参数ID（1为平面应力，2为平面应变）%输出单元刚度矩阵