clc;
clear all;

c = 0:0.025:2;
b = 0:0.025:2;
d = 0:0.025:2;
a = c.*(b).^2 + d.^3;


samples_c = normrnd(2, 0.5, 200, 1);
std_c = std(samples_c);
samples_b = normrnd(5, 0.1, 200, 1);
std_b = std(samples_b);
samples_d = normrnd(9, 0.7, 200, 1);
std_d = std(samples_d);

absolut_a = abs(sqrt(((2 * 5^2) * sqrt(((std_c./samples_c).^2) + ((2*5./samples_b).^2)).^2 )+(3*9*std_d)^2));
meam_absoolut_a = mean(absolut_a)
vidnosna_a = absolut_a./meam_absoolut_a

subplot(3,1,1)
plot(samples_c);hold on;
plot(samples_b);hold on;
plot(samples_d);hold on;
title("Вхідні значення")
legend('Значення c', "Значення b", "Значення d")

subplot(3,1,2)
plot(vidnosna_a);hold on;
plot(absolut_a)
title("Абсолютна та відносна похибки")
legend('Відносна похибка', "Абсолютна похибка")

subplot(3,1,3)
plot(a)
title("Залежність вихідної похибки від вхідних похибок")