例2-63.m

这是一个MATLAB7.0基础与提高例题的所有源码

%用trapz在区间-1<x<3上计算y=humps(x)下面的面积
>> x=-1:0.12:3;                     %粗略计算积分产生x的序列
>> y=humps(x);
>> area=trapz(x,y)                   %计算积分
>> x=-1:0.02:3;                      %以更好的精度计算
>> y=humps(x);
>> area=trapz(x,y)
%用quad和quadl在区间-1<x<3上计算y=humps(x)下面的面积
>> format long
>> area=quad('humps',-1,3)
>> area=quadl('humps',-1,3)