MATLAB习题答案_2007(张志涌著)北京航空航天大学出版社

d=0.1;tt=0:d:10;t=tt+(tt==0)*eps;y=sin(t)./t;s=d*trapz(y)%计算出积分值ss=d*(cumtrapz(y))%计算累计积分并积分曲线plot(t,y,t,ss,'r'),holdon%插值法计算y(4.5)yi=interp1(t,ss,4.5)plot(4.5,yi,'r+')%符号解法symsxty1y2y1i,y1=sin(t)./t,y1i=int(y1,t,0,x),y2=subs(y1i,x,4.5)holdon,plot(4.5,y2,'*m')