quadl_humps.m

matlab从入门到精通第6章书后源码对初学者很有帮助

% quadl_humps.m
% 积分求解humps函数的积分解
% 通过函数cumtrapz累计计算humps函数的积分数值序列
x=linspace(-1,2,100);
y=humps(x);
z=cumtrapz(x,y);
[AX,H1,H2]=plotyy(x,y,x,z,'plot');
set(get(AX(1),'Ylabel'),'String','humps(x)')
set(get(AX(2),'Ylabel'),'String','integral of humps(x)')
grid on;xlabel('x');
title('cumulative integral of humps(x)');
% 通过quad和quadl函数求积分
humps_hdl=@humps;
humps1=quad(humps_hdl,-1,2)
humps2=quadl(humps_hdl,-1,2)