exm040811_1.m

同上

syms x;IS=int('exp(-x*x)','x',0,1)	%求解析积分
vpa(IS)									%求所得解析积分的32位精度近似值  

fun=inline('exp(-x.*x)','x');		%注意：数组乘符号.*的采用是必须的。
Isim=quad(fun,0,1),IL=quadl(fun,0,1)  

Ig=gauss10(fun,0,1)  

xx=0:0.1:1.5;ff=exp(-xx.^2);		%产生被积函数的“表格”数据
pp=spline(xx,ff);					%由“表格”数据构成样条函数
int_pp=fnint(pp);					%求样条积分
Ssp=ppval(int_pp,[0,1])*[-1;1]	%据样条函数计算[0,1]区间的定