下面给出了根据Simpson的1-3（h / 3）方法找到以下积分的C代码。适当地填写代码上的空白。我想在下面用Matlab解决这个问题，但是我没有做。这是一个简单的问题，但我做不到。如果有人能帮助我，我会很高兴。
[C代码版本2
我在Matlab中尝试了以下代码块：
% Ask for user input
% Lower bound (a)
a = input('What is your lower bound (a)?')
% Upper bound (b)
b = input('What is your upper bound (b)?')
% Subintervals
N = input('How many subintervals (N)?')
% Defining function
f = @(x,e) (e*x+sin(x))
% Finding h
h=(b-a)/N;
% Finding the values of x for each interval
x=linspace(a,b,N);
% Calculating the integral
for i = 1:N-1
I(i)= (h/3)*(f(x(i))+(4*f((x(i)+x(i+1))/2))+f(x(i+1)));
end
answer1 = sum(I)
disp(I)
% adding f(b) to sum
val2=ff(length(xx));
sum=val1+val2+sum;% set sum
% set result
result=sum*h/3;
请注意，MATLAB不会将符号e用作尼泊尔编号（欧拉编号）。为了产生MATLAB欧拉数，你可以使用指数函数exp(x)，e = exp(1)，
因此，首先，请更正您的函数定义：
F = @(x) exp(1).^x + sin(x) % Always try to use Upper-Case letters for your variable/function name
然后，您可以使用以下代码段使用Simpson的1/3计算积分：
a = 0; b = 3; N = 1e4;
F = @(x) exp(1).^x + sin(x);
h = ((b-a)/2)/N;
x = linspace(a,b,N);
I = 0;
for i = 1:N-1
I = I + h/3*(F(x(i)) + 4*F((x(i)+x(i+1))/2) + F(x(i+1)));
end
disp(I)
% To compare your result:
Itz = trapz(x, F(x))
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句