二维有限元解法的程序实现