对于离散的数值积分，只能用梯形积分法来求解。而dblquad则需要具体函数表达式。求解命令函数，trapz（）求解格式，z=trapz(x,y) %注：y与x是同维数的向量。例：x=-1:0.5:1;y=[0.367879 0.778801 1 0.778801 0.367879];z=trapz(x,y)
这是不行的。但可用trapz这个函数实现数值积分。
