四节点等参元matlab程序,采用高斯积分3*3运算,这个是主程序