function kg = trout_vector ( y )
  kg = 25.1 * ( 1.0 - exp ( - 1.19 * y ) ) .^ 3;
  return
end
