用MATLAB程序编写复化梯形公式和复化simpsion公式求解积分

%Compound trapezoid formula复化梯形法

function y = traint(a,b,n,f)

h = (b - a) /搜旦 n

x = linspace(a,b,n+1)

y1 = h * feval(f,x)

y1(1) = y1(1) / 2

y1(n+1) = y1(n+1) /罩漏袭 2

y = sum(y1)

%compound Simpson formula复化辛普森公式

function y = sraint(a,b,n,f)

h = (b - a) / n

x = linspace(a,b,2*n+1)

y1 = feval(f,x)

y1(2:2:2*n) = 4 * y1(2:2:2*n)

y1(3:2:2*n-1) = 2 * y1(3:2:2*n-1)

y = h / 6 * sum(y1)

>>f=inline('4./(1+x.*x)')

>>traint(-1,1,8,f)

ans =

6.2624

>物兄>sraint(-1,1,8,f)

ans =

6.2832

首先，你需梁哗提供已知的f（x）函数，再进行数值计算。

利用Matlab软件，使用复化梯形公式计算数值积分的方法橡悄行如下。

1、自定义复化梯形法函数，traint(）

function y = traint(a,b,n,func)

h = (b - a) / n

x = linspace(a,b,n+1)

y1 = h * feval(func,x)

y1(1) = y1(1) / 2

y1(n+1) = y1(n+1) / 2

y = sum(y1)

end

2、自定义已知的f（x）函数，func(）

function y = func(x)

y=f（x） %要具体的函数表达式

end

3、运扒在命令窗口中，输入

a=1.8b=2.6n=10 %回车

traint(-1,0,8,f) %回车

这个是三维绘图，构建x－y平面上的“格点”矩阵。即［X，Y］＝meshgrid（xjny）你好像差了一个变量，就是变量z，变量z是和变棚扮量X，Y有关的，也就是说它是（X，郑锋Y）的函数5173之后才能做出来三维图形链丛灶只有两个变量只能做出二维图形！

---