function c = hyp(a, b)
  c = sqrt(a*a + b*b);