function area = traparea(a, b, h)
  area = 0.5 * (a + b) * h;
end