例7-67  quad函數使用自適應法(Simpson法)來求解積分。
函數 
先編制M文件如下。
function f=fun1(x)
f=x.^2./(1+sin(x)+x.^2);
end
將該M文件以fun1.m為函數名保存。繼續在命令窗口中輸入如下命令，并按Enter鍵確認。
>> quad('fun1',0,1)
ans =
    0.2301
>>