How can I switch from Simpson code to composite simpson code?

1 次查看（过去 30 天）

显示更早的评论

Hazel Can 2022-12-25

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1884277-how-can-i-switch-from-simpson-code-to-composite-simpson-code

评论： Torsten 2022-12-25

在 MATLAB Online 中打开

Composite simpson rule :

Simpson rule's code :

% MATLAB code for syms function that creates a variable
% dynamically and automatically assigns
% to a MATLAB variable with the same name
syms x
  
% Lower Limit
a = 4;
  
% Upper Limit
b = 5.2;
  
% Number of Segments
n = 6;
  
% Declare the function
f1 = log(x);
  
% inline creates a function of string containing in f1
f = inline(f1);
  
% h is the segment size
h = (b - a)/n;
  
% X stores the summation of first and last segment
X = f(a)+f(b);
  
% variables Odd and Even to store 
% summation of odd and even
% terms respectively
Odd = 0;
Even = 0;
for i = 1:2:n-1
    xi=a+(i*h);
    Odd=Odd+f(xi);
end
for i = 2:2:n-2
    xi=a+(i*h);
    Even=Even+f(xi);
end
  
% Formula to calculate numerical integration 
% using Simpsons 1/3 Rule
I = (h/3)*(X+4*Odd+2*Even);
  
disp('The approximation of above integral is: ');
disp(I);