例6-28  quad8函數使用自適應法(Simpson法)來求解積分。
函數 
先編寫M文件如下。
function f=fun2(x)
f=x.^2./exp(-x);
end
在命令窗口中輸入如下命令，并按Enter鍵確認。
>> quad8('fun2',0,3)
ans =
   98.4277
>> a=quad8('fun2',1,3,[1e-10,1e-11])
a =
  97.7094
>> vpa(a,10)
ans =
97.70940279
>> quad8('fun2',1,3,1e-6,1)
      18     1.0000000000    1.00000000e+000    97.7094029629
ans =
   97.7094
>>