fx.m

来自「简单实现了复化梯形公式的求解」· M 代码 · 共 26 行

clear
clc
a=0;
b=1;
n=100;
start=10;%最小分段数
x=sym('x','real');
f=x^2;
jf=double(int(f,x,a,b));

% for i=10:n
%     cc=ft(a,b,i)-jf;
%     plot(i,cc,'*-r')
%    hold on    
% end
c=zeros(10,1);
for i=1:n
    c(i)=ft(a,b,i+start)-jf;
end
x=1+start:n+start;
plot(x,c,'-r')
title('误差和分段之间的关系图')
xlabel('分割的段数')
ylabel('误差')