% ENM 8.4
% Denna kod skrevs på övningen 2017-09-28 i SF1545

clear

% Given funktion
p = @(a,b,c) (a+b+c)/2;
T = @(a,b,c) sqrt(p(a,b,c)*(p(a,b,c)-a)*(p(a,b,c)-b)*(p(a,b,c)-c));

a = 12.5; b = 10.25; c = 2.75;
Area = T(a,b,c) % ostörda värdet

% Störningsräkning
% Stör en parameter i taget
Ea = 0.005; Eb = Ea; Ec = Ea;
part1 = abs(T(a+Ea,b,c) - Area)
part2 = abs(T(a,b+Eb,c) - Area)
part3 = abs(T(a,b,c+Ec) - Area)
ET = part1 + part2 + part3

% (Här kan vi också se vilken parameter som bidrar mest till
% osäkerheten. I detta fall tydligen c, men skillnaden är inte så
% stor mellan bidragen.)
