复化梯形公式matlab实验报告

复化梯形公式matlab实验报告

资源描述

复化积分法(复化梯形求积，复化Simpson公式，变步长求积法)MATLAB编程实验报告

一、 问题描述：

编写函数实现复化积分法。

二、 实验步骤(过程)：

(1) 复化求积法

(1) 复化梯形求积：用复化梯形求积公式求解

function [f]=Tn(a,b,n,y)

syms t;

h=(b-a)/n;

f=0;

for k=1:n+1

x(k)=a+(k-1)*h

z(k)=subs(y,t,x(k));

end

for i=2:n

f=f+z(i);

end

q=subs(y,t,a);

if y==sin(t)/t&&a==0

q=1;

end

p=subs(y,t,b);

T=h/2*(q+p+2*f);

T=vpa(T,7)

clc,clear;

syms t;

a=0;b=1;

y=sin(t)/t;

n=8;

Tn(a,b,n,y);

(2) 复化Simpson公式：用复化Simpson公式求解

function [f]=simpson(a,b,n,y)

syms t;

h=(b-a)/n;

f=0;l=0;