I am trying to solve PDE containing integral. Its form is showed below:

D[Pa[t],t] == Integrate[-k[x]*Pa[t],{x,0,1}]
k[x_] := 2*x

with initial conditions:

Pa[0] == 100

And NDSolve command gives fine result:

sol1 = NDSolve[{D[Pa[t], t] == Integrate[-k[x]*Pa[t], {x, 0, 1}], Pa[0] == 100}, Pa, {t, 0, 5}]

However, is it possible to get rid of integral and solve the equation as PDE:

D[Pa[t,x],t,x] == -k[x]*Pa[t,x]

What kind of Initial conditions are needed in this case?


1 Answer 1


The first part of your question is a differential equation

D[Pa[t],t] == Integrate[-k[x],{x,0,1}]*Pa[t]

which doesn't depend on x directly and can be easily solved analytical

Pa[t]== Pa0 Exp[-C t] 

The last part of your question concerning the pde might be solved with the separation of variables


Your pde evaluates to

X'[x]T'[t] == -k[x]*X[x]T[t] 

which can be separated to the form

T'[t]/T[t] == -((k[x]*X[x])/X'[x]) ==  -const

These two ode's can be solved analytically. const must be choosen to fullfill your boundary conditions X[0] and X[1] ...


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.