复合梯形方法-matlab实现

数值分析数值积分-复合梯形方法的matlab实现

%复合梯形公式
a = 1;b = 3;
m = 5000;


h = （b-a）/m;
x = a:h:b;
y = （100./（x.^2））.*（sin（10./x））;
Im = h/2*（y（1）+y（m+1）+2*sum（y（2:m）））;

%变步长梯形公式
n = 1;
h1 = b-a;
x1 = [ab];
y1 = （100./（x1.^2））.*（sin（10./x1））;
s1 = h1/2*（y1（1）+y1（2））;
for i = 1:10000
    n = 2*n;
    h2 = h1/2;
    h1 = h2;
    x1 = a:h2:b;
    y1 = （100./（x1.^2））.*（sin（10./x1））;
    s2 = s1/2 + h2*sum（y1（2:2:n））;
    if abs（s2-s1）< 0.5*10^（-6）
        break;
    end
    s1 = s2;
end
s1
n
i
    
    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         537  2011-11-14 01:05  E6_韩明岗\O_R.m
     文件         509  2011-11-10 11:34  E6_韩明岗\fuhe_tixing.m
     目录           0  2011-11-14 12:31  E6_韩明岗\