复化积分法（复化梯形求积，复化Simpson公式，变步长求积法）MATLAB编程实验报告一、问题描述：编写函数实现复化积分法二实验步骤(过程)：(一)复化求积法(1)复化梯形求积：用复化梯形求积公式求解J1聾dx0xfunction[f]二Tn(a,b,n,y)symst;h=(b-a)/n;f=0;fork=1:n+1x(k)=a+(k-1)*hz(k)二subs(y,t,x(k));endfori=2:nf=f+z(i);endq=subs(y,t,a);ify=='sin(t)/t'&&a==0q=1;endp二subs(y,t,b);T二h/2*(q+p+2*f);T=vpa(T,7)clc,clear;symst;a=0;b=1;y=sin(t)/t;n=8;Tn(a,b,n,y);(2)复化Simpson公式：用复化Simpson公式求解€2exdxifunction[f]二simpson(a,b,n,y)symst;h=(b-a)/n;f=0;l=0;fork=1:n+1x(k)=a+(k-1)*hw(k)=0.5*h+x(k)z(k)=subs(y,t,x(k));endfori=2:nf二f+z(i);endfori二上nl=l+w(i);endq二subs(y,t,a);ify=='sin(t)/t'&&a==0q=1;endp=subs(y,t,b);T=h/2*(q+p+2*f);T=vpa(T,7)clc,clear;symst;a=1;b=2;y二exp(1/t);n=5;simpson(a,b,n,y);(3)变步长求积法：以书本例4.5为例function[f]二TN(a,b,y,R0)。