输入函数变量:
>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,10,'mjx') 按回车键显示结果如下: f=0.7856 n =10
(2)100等份(n=100) 输入函数变量:
>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,100,'mjx') 按回车键显示结果: f =0.7854 n =100
(3)1000等份(n=1000) 输入函数变量:
>>f=inline('1/(1+x*x)');
>>[f,n]=tyszjf(f,0,1,1000,'mjx') 按回车键显示结果如下: f =0.7854 n =1000
3.2.2梯形公式调用实例(计算不同等分份数的积分值)
(一)10等份 输入函数变量:
>>f=inline('1/(1+x*x)');
7
>>[f,n]=tyszjf(f,0,1,10,'tx') 按回车键显示结果如下: f =0.7850 n =10
(二)100等份 输入函数变量:
>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,100,'tx') 按回车键显示结果如下: f =0.7854 n =100
(三)1000等份 输入函数变量:
>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,1000,'tx') 按回车键显示结果如下: f =0.7854 n=1000
3.2.3辛普森公式调用实例
(一)10等份 输入函数变量:
>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,10,'sps') 按回车键显示结果如下:
8
f =0.7854 n =10
(二)100等份 输入函数变量:
>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,100,'sps') 按回车键显示结果如下: f =0.7854 n =100
(三)1000等份 输入函数变量:
>>f=inline('1/(1+x*x)');
>>[f,n]=tyszjf(f,0,1,1000,'sps') 按回车键显示结果如下: f =0.7854 n =1000
(四)精度要求不同时
(1)wc=0.5*10^(-7) 输入函数变量:
>>f=inline('1/(1+x*x)');
>>[f,n]=tyszjf(f,0,1,n,'sps',0.5*10^(-7)) 按回车键显示结果如下: f =0.7854 n =100
(2)wc=0.5*10^(-8)
9
输入函数变量:
>>f=inline('1/(1+x*x)');
>>[f,n]=tyszjf(f,0,1,n,'sps',0.5*10^(-8)) 按回车键显示结果如下: f =0.7854 n =100
3.2.4牛顿-柯茨公式调用实例
(一)10等份 输入函数变量:
>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,10,'nc') 按回车键显示结果如下: f =0.7854 n =10
(二)100等份 输入函数变量:
>>f=inline('1/(1+x*x)'); >>[f,n]=tyszjf(f,0,1,100,'nc') 按回车键显示结果如下: f =0.7854 n =100
(三)1000等份 输入函数变量:
>>f=inline('1/(1+x*x)');
10
>>[f,n]=tyszjf(f,0,1,1000,'nc') 按回车键显示结果如下: f =0.7854 n =1000
3.2.5自适应积分调用实例
(一)wc=0.5*10^(-6) 输入函数变量:
>>f=inline('1/(1+x*x)');
>>[f,n]=tyszjf(f,0,1,n,'zsy',0.5*10^(-6)) 按回车键显示结果如下: f=0.7854 n=8
(二)wc=0.5*10^(-7) 输入函数变量:
>>f=inline('1/(1+x*x)');
>>[f,n]=tyszjf(f,0,1,n,'zsy',0.5*10^(-7)) 按回车键显示结果如下: f =0.7854 n =8
(三)wc=0.5*10^(-5) 输入函数变量:
>>f=inline('1/(1+x*x)');
>>[f,n]=tyszjf(f,0,1,n,'zsy',0.5*10^(-5)) 按回车键显示结果如下:
11