求教高手一个编程的问题！！

star198311

以下是一个复合simpson求积分的函数。用函数求解是f=inline('sin(x)<FONT color=#000000>. </FONT>/x','x')。但这个函数不能求解函数sin(x)/x在[0 1]间的积分。因为它在0点的值用feval（）算的结果是inf。这个问题怎么解决？function [Isic]=simpsonc(a,b,M,f,varargin)<BR>H=(b-a)/M;<BR>x=linspace(a,b,M+1)<BR>fpm=feval(f,x,varargin{:});<BR>fpm(2:end-1)=2*fpm(2:end-1);<BR>Isic=H*sum(fpm)/6;<BR>x=linspace(a+H/2,b-H/2,M);<BR>fpm=feval(f,x,varargin{:});<BR>Isic=Isic+2*H*sum(fpm)/3;