function w = heron(a, b, c)
p = (a + b + c) / 2;
w = sqrt(p*(p-a)*(p-b)*(p-c));
end

