慕工大数值分析第五课第一题

clc;
clear;

a = 0;
b = 4;

%% 5.1.1 Mittelpunktregel

I_M = (b - a) .* (f((b - a)./ 2));

%% 5.1.2 Trapezregel

I_T = (b - a) .* (f(a) + f(b)) ./ 2;

%% Function

function y = f(x)
y = power(x, 5) ./ power((1+x), 5);
end

%% 1. Gauß-Quadratur mit einer Stützstelle

I_G_1 = (b - a) ./2 .* sum(gw(1) .* f(((b - a) ./ 2) .* gx(1) + (a + b) ./ 2));

%% 2. Gauß-Quadratur mit zwei Stützstellen

I_G_2 = (b - a) ./2 .* sum(gw(2) .* f(((b - a) ./ 2) .* gx(2) + (a + b) ./ 2));

%% 3. Gauß-Quadratur mit drei Stützstellen

I_G_3 = (b - a) ./2 .* sum(gw(3) .* f(((b - a) ./ 2) .* gx(3) + (a + b) ./ 2));

%% Vergleichender Fehler mit der exakten Integration I = 0.556543771162832.

I = 0.556543771162832;
Fehler_1 = abs(I - I_G_1);
Fehler_2 = abs(I - I_G_2);
Fehler_3 = abs(I - I_G_3);

%% Function

function y = f(x)
y = (x ./ (1 + x)) .^ 5;
end

function gaussw = gw(n)

if n == 1

gaussw = 2;

elseif n == 2

gaussw = [1, 1];

elseif n == 3

gaussw = [5/9, 8/9, 5/9];

else

disp('false n')

end

function gaussx = gx(n)

if n == 1

gaussx = 0;

elseif n == 2

gaussx = [-1/sqrt(3), 1/sqrt(3)];

elseif n == 3

gaussx = [-sqrt(3/5), 0, sqrt(3/5)];

else

disp('false n')

end