Mathematica Stack Exchange is a question and answer site for users of Mathematica. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I am now trying to solve second - order linear partial differential equation in that interested eq. has been separated into variables to simplify the procedure for Mathematica. Here is my equation;

Differential Equation

or separated form,

Separated Form

Also, the last equation is coded in Mathematica as;

    L[t_, z_] := T[t] Z[z]
    Eq = (D[T[t], t]/(A*T[t])) == R*(D[D[Z[z], z], z]/Z[z]) - B*(D[Z[z], z]/Z[z])
    L[t, z] /. Flatten@{DSolve[Eq[[1]] == F, T, t, GeneratedParameters -> CT], 
    DSolve[Eq[[2]] == F, Z, z, GeneratedParameters -> CS]}

Then, the solution is given as;

The Solution

However at that point, I should plug IC and BCs into related eq. to both get rid of constants and proceed with my calculations. These are my IC and BCs, respectively;




How can I solve this problem? Is there any way to insert these IC and BCs into the code above? or Does there exist any method, command etc. to come up with solution without integration constants?

Note: Any of the values included in BCs is neither a function of z nor that of t.

share|improve this question
Enter your IC and BCs at the same time as your eqn, i.e. as a list of eqns. There are many examples in the documentation for DSolve. But of course you need to transform your IC and BCs in the same way that you transformed the eqn. Try doing that and then let us know if you have problems. – Mike Honeychurch Jul 16 '14 at 22:04

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.