Question

在MATLAB用辛普森法则计算x^2+3x在区间[1,2]上面的积分值

感谢大佬

Answer 1

用辛普森法则计算，可以得到 x≈ 6.8333。实现代码
function simpson_main()
clc,close all
f=inline('x.^2+3*x','x');
simpson(1,2,10,f)
辛普森公式代码内容：
h = (b - a) / n;
x = linspace(a,b,2*n+1);
y1 = feval(f,x);
y1(2:2:2*n) = 4 * y1(2:2:2*n);
y1(3:2:2*n-1) = 2 * y1(3:2:2*n-1);
y = h / 6 * sum(y1);