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);