Hello everybody.

I have a function:

f(x,y)=GAMMA(y, -ln(x))/GAMMA(y)

seq(sum(f(x, y), y = 0 .. 1), x = 0 .. 5)

and I got a error message:

**Error, (in ln) numeric exception: division by zero ??**

This is normal behavior in** seq **function or** Bug**?

but when I'm first calculate the sum sol := sum(f(x, y), y = 0 .. 1) -> x,

and evalf([seq(sol, x = 0 .. 5)]) ->[0., 1., 2., 3., 4., 5.] works fine.

Seq-division_by_zero.mw

Mariusz Iwaniuk