fightlive的博客

原创 数值分析 MATLAB实现Romberg积分法（最大N=16）

使用Romberg积分法求解其他积分函数可以按注释对代码进行修改。积分方法为先使用复化梯形公式求解T序列，再使用Romberg法求解复化simpson，即S序列，再使用Romberg法求解复化cotes，即C序列，最后同理求R序列。代码支持N=1，2，4，8，16;再添加需要按模块手动添加。function Romberg_myx1=1:2;x2=1:3;x3=1:5;x4=1:9;x5=1:17;x0=1;%x0，xn为积分上下限xn=3;x1=[x0,xn];x2(