Matlab 中的 Simpson's 1-3 (h / 3) 方法

问题描述

下面给出了根据辛普森 1-3 (h / 3) 方法求出下列积分的 C 代码。适当填写代码上的空白。 我想在 Matlab 中解决下面的这个问题，但我没有这样做。这是一个简单的问题，但我做不到。如果有人帮助我，我会很高兴。

C code version

[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;

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）