梯形积分及辛普生积分

看看是不是这样的：

function zint
clear all;clc;
a=0;b=1;d=0.05;
x=a:d:b;y=func(x);
format long
y_trapz=trapz(y)*d %梯形积分
y_simps=quad(@func,a,b) %辛普生积分
y_accurt=eval(int('exp(-t)',a,b))% 精确值
%误差
y_trapz_error=y_trapz-y_accurt
y_simps_error=y_simps-y_accurt
%图形
plot(x,exp(-x),'o-')

function y=func(x)
y=exp(-x);

结果：
y_trapz = 0.632252245124816
y_simps = 0.632120560205377
y_accrt = 0.632120558828558

y_trapz_error = 1.316862962582599e-004
y_simps_error = 1.376819436238463e-009

大按时代发粉我下 按时打发上；打发上；地方按时地方按时地方；老师傅拉；四六风；啊；色浪费；挖人；‘大师傅山豆根 士大夫敢死队风格

计算方法的书里有写，上学学过忘了