s := 0;
  for j from 1 to 10 do
     s := s + i*y*t^j/(y + j*t)^j
  end_for;
  normal(s);
  purge(s);