m0_50772268的博客

原创 自适应求积算法 MatLab版

%自适应simpson算法 clc clear close all; syms x f(x)=100/(x^2)*sin(10/x);%目标函数 a=1;b=3;%初始区间 TOL=1E-4;%误差容限 [s,fcount]=acsimpson1(f,a,b,TOL); t = 1:0.001:3; plot(t,f(t)) hold on disp("自适应求积结果:") result = vpa(s); fprintf("result= %f",result); % 画出f(x)图像 functio.