function v = f508(t, y)
  global a
  v = a*y;
  