[fanxing39]

上面的方法是
1.先产生x，y的网格矩阵
2.将（x，y）每个点的 f（x，y）利用quad （数值积分）进行求解
3. surf （P） 画出曲面图

程序里用了2个for 循环，这个是我不希望看到的，但是又没有找到更好的方法，这里修改下的话，运算速度才比较快！