Octave / task_8

#19
a = 0.0;
b = 1.8;
h = (b - a) / 30;
x = [a : h : b];
ydop = x .^ 2 .* sqrt(4 - x .^ 2);
y = @(x)(x ^ 2 * sqrt(4 - x ^ 2));
res = y(x(1)) + y(x(31));
for i = 2 : 1 : 30 #в цикле от 2 до 30 (т. е. от 1 до 29)
if ((i - 1) / 2 != floor((i - 1) / 2))
res += 4 * y(x(i));
else
res += 2 * y(x(i));
endif
endfor
res *= h / 3;
disp("Значение интеграла по формуле Симпсона: ");
disp(res);
disp("Значение интеграла, вычисленное методом трапеций:");
disp(trapz(x, ydop));
disp("Значение интеграла, вычисленное методом квадратур:");
disp(quad(y, a, b, [0.0006, 0.0006]));