Hi all,

How to calculate this integral:

for k>0,m>0

Int(exp(-(1/2)*v/k)*v^3*exp((1/2)*v/m)*Ei(1, -(1000*I)*v+(1/2)*v/m), v = 0 .. infinity)

I'm  tried to take advantage of with(IntegrationTools) but I failed

and and I got a strange result ,like this:

Integral.mw

Hello.

I have a Pde solution in from of the sum.

pde := diff(u(x, t), t) = diff(u(x, t), x$2)

symbolic := pdsolve([pde, u(x, 0) = 1, u(0, t) = 0, u(1, t) = 0])

symbolic := u(x, t) = Sum(-(2*((-1)^_Z9-1))*sin(_Z9*Pi*x)*exp(-Pi^2*_Z9^2*t)/(Pi*_Z9), _Z9 = 1 .. infinity)

I tried a subs or eval command dosen't work.

Thanks.

pdex1.mw
 

Download pdex1.mw

Example 6 from Maple Help.

restart:

with(Optimization):

LPSolve(2*x+5*y, {3*x-y = 1, x-y <= 5}, assume = {integer, nonnegative})

Kernel crash the same with Maple 2015.

Does anyone can confirm?

Mariusz Iwaniuk

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

Hello guys,

I was just playing around with differential equations, when I noticed that symbolic solution is  different from the numerical.What is the reason for this strange behavior?

ODE := (diff(y(x), x))*(ln(y(x))+x) = 1

sol := dsolve({ODE, y(1) = 1}, y(x))

a := plot(op(2, sol), x = .75 .. 2, color = "Red");
sol2 := dsolve([ODE, y(1) = 1], numeric, range = .75 .. 2);

with(plots);
b := odeplot(sol2, .75 .. 2, thickness = 4);
display({a, b});

Strange_issue.mw

Mariusz Iwaniuk