复化梯形积分公式matlab程序

function[y]=ft(x)y=exp(-x*x);%此处放置要积分的函数function[s]=fhtxinga=0.2;b=1.5;n=8;%a为下限，b为上限，n为要划分的段数h=(b-a)/n;s=0;fork=1:nXk=a+h*(k-1);Xkk=a+h*(k);s=s+ft(Xk)+ft(Xkk);ends=h/2*s;