I have the following task to do:

I developed the following code:
function I = simpson(f,a,b,n)
%I = SIMPSON(f,a,b,n) przybliża wartość całki funkcji f
%na przedziale [a,b] stosując metodę Simpsona z n podprzedziałami

x = linspace(a, b, n+1); % Określ listę wektor punktów przy użyciu LINSPACE
Delta = (b – a) / n; % Określ stałą Delta

% Zakładamy wstępną wartość całki jako 0, a następnie obliczamy sumę
I = 0;
for i = 1:n % Ustaw zakres dla i w pętli
dI = f(a(i)) + 4 * f(a(i) + a(i+1) / 2) + f(a(i+1)); % Określ składnik sumy dla i-tego podprzedziału
% Zwiększa sumę całkowitą o nową wartość
I = I + dI;
end
% Mnoży sumę przez Delta
I = Delta*I;
end
Auxiliary function:
% Test 1
a = 0; b = 1; n = 10; f = @(x)1-x+x.^3;
I = simpson(f,a,b,n)
fprintf(‘Test1: Błąd trapezów: %e | Błąd Simpsona: %e’,abs(trapz((b-a)/n,f(linspace(a,b,n+1)))-0.75),abs(I-0.75))

% Test 2
a = 0; b = pi/2; n = 50; f = @(x)sin(x);
I = simpson(f,a,b,n)
fprintf(‘Test2: Błąd trapezów: %e | Błąd Simpsona: %e’,abs(trapz((b-a)/n,f(linspace(a,b,n+1)))-1),abs(I-1))
Error message:

Where is the problem?I have the following task to do:

I developed the following code:
function I = simpson(f,a,b,n)
%I = SIMPSON(f,a,b,n) przybliża wartość całki funkcji f
%na przedziale [a,b] stosując metodę Simpsona z n podprzedziałami

x = linspace(a, b, n+1); % Określ listę wektor punktów przy użyciu LINSPACE
Delta = (b – a) / n; % Określ stałą Delta

% Zakładamy wstępną wartość całki jako 0, a następnie obliczamy sumę
I = 0;
for i = 1:n % Ustaw zakres dla i w pętli
dI = f(a(i)) + 4 * f(a(i) + a(i+1) / 2) + f(a(i+1)); % Określ składnik sumy dla i-tego podprzedziału
% Zwiększa sumę całkowitą o nową wartość
I = I + dI;
end
% Mnoży sumę przez Delta
I = Delta*I;
end
Auxiliary function:
% Test 1
a = 0; b = 1; n = 10; f = @(x)1-x+x.^3;
I = simpson(f,a,b,n)
fprintf(‘Test1: Błąd trapezów: %e | Błąd Simpsona: %e’,abs(trapz((b-a)/n,f(linspace(a,b,n+1)))-0.75),abs(I-0.75))

% Test 2
a = 0; b = pi/2; n = 50; f = @(x)sin(x);
I = simpson(f,a,b,n)
fprintf(‘Test2: Błąd trapezów: %e | Błąd Simpsona: %e’,abs(trapz((b-a)/n,f(linspace(a,b,n+1)))-1),abs(I-1))
Error message:

Where is the problem? I have the following task to do:

I developed the following code:
function I = simpson(f,a,b,n)
%I = SIMPSON(f,a,b,n) przybliża wartość całki funkcji f
%na przedziale [a,b] stosując metodę Simpsona z n podprzedziałami

x = linspace(a, b, n+1); % Określ listę wektor punktów przy użyciu LINSPACE
Delta = (b – a) / n; % Określ stałą Delta

% Zakładamy wstępną wartość całki jako 0, a następnie obliczamy sumę
I = 0;
for i = 1:n % Ustaw zakres dla i w pętli
dI = f(a(i)) + 4 * f(a(i) + a(i+1) / 2) + f(a(i+1)); % Określ składnik sumy dla i-tego podprzedziału
% Zwiększa sumę całkowitą o nową wartość
I = I + dI;
end
% Mnoży sumę przez Delta
I = Delta*I;
end
Auxiliary function:
% Test 1
a = 0; b = 1; n = 10; f = @(x)1-x+x.^3;
I = simpson(f,a,b,n)
fprintf(‘Test1: Błąd trapezów: %e | Błąd Simpsona: %e’,abs(trapz((b-a)/n,f(linspace(a,b,n+1)))-0.75),abs(I-0.75))

% Test 2
a = 0; b = pi/2; n = 50; f = @(x)sin(x);
I = simpson(f,a,b,n)
fprintf(‘Test2: Błąd trapezów: %e | Błąd Simpsona: %e’,abs(trapz((b-a)/n,f(linspace(a,b,n+1)))-1),abs(I-1))
Error message:

Where is the problem? #simpson integral MATLAB Answers — New Questions

​