Hello, Mr. Rogers,
Using DSolve[eq13, q, {x, s}, GeneratedParameters -> B1] will produce a pure function and the following calculation will be based on function.
But using DSolve[eq13, q[x,s], {x, s}, GeneratedParameters -> B1] will produce an expression and the following calculation will be based on substitution.
Is it correct?
So if I redefine q[x,s] and it can work likewise.
In[45]:= q[x_,
s_] := -(E^(-((Sqrt[s] x)/Sqrt[Dp]) -
Sqrt[s/Dp] (lh +
x)) (-(E^(
Sqrt[s/Dp] (2 lh + x)) (-Sqrt[s] +
Sqrt[Dp] (Sqrt[s/Dp] - \[Alpha]1)) +
E^(Sqrt[s/Dp] (lh + 2 x)) (Sqrt[s] + Sqrt[Dp] \[Alpha]1) +
E^(lh Sqrt[s/Dp] + (2 Sqrt[s] x)/Sqrt[
Dp]) (Sqrt[s] + Sqrt[Dp] \[Alpha]1) +
E^((2 lh Sqrt[s])/Sqrt[Dp] +
Sqrt[s/Dp] x) (Sqrt[s] -
Sqrt[Dp] (Sqrt[s/Dp] + \[Alpha]1))) Sqrt[
Dp s (s + \[Theta]c)] (-Dc s + Dp (s + \[Theta]c)) (-A2 +
pinf (s + \[Theta]c)) -
A1 Sqrt[Dp s] (s + \[Theta]c) (Dp^(
3/2) (E^(Sqrt[s/Dp] (2 lh + x)) -
E^((2 lh Sqrt[s])/Sqrt[Dp] + Sqrt[s/Dp] x)) Sqrt[s/
Dp] (s + \[Theta]c) +
Dc (E^(lh Sqrt[s/Dp] + (2 Sqrt[s] x)/Sqrt[
Dp] + (lh - x) Sqrt[(s + \[Theta]c)/Dc]) (Sqrt[s] +
Sqrt[Dp] \[Alpha]1) (-s +
Dp Sqrt[s/Dp] Sqrt[(s + \[Theta]c)/Dc]) -
E^(lh Sqrt[s/Dp] +
2 Sqrt[s/Dp] x + (lh - x) Sqrt[(s + \[Theta]c)/
Dc]) (Sqrt[s] + Sqrt[Dp] \[Alpha]1) (s +
Dp Sqrt[s/Dp] Sqrt[(s + \[Theta]c)/Dc]) +
E^((2 lh Sqrt[s])/Sqrt[Dp] +
Sqrt[s/Dp] x) (-s^(3/2) +
Sqrt[Dp] s (Sqrt[s/Dp] + \[Alpha]1) +
Dp Sqrt[s] Sqrt[s/Dp] Sqrt[(s + \[Theta]c)/Dc] -
Dp^(3/2) Sqrt[s/Dp] \[Alpha]1 Sqrt[(s + \[Theta]c)/
Dc]) + E^(
Sqrt[s/Dp] (2 lh + x)) (s^(3/2) +
Sqrt[Dp] s (-Sqrt[(s/Dp)] + \[Alpha]1) +
Dp Sqrt[s] Sqrt[s/Dp] Sqrt[(s + \[Theta]c)/Dc] +
Dp^(3/2) Sqrt[s/Dp] \[Alpha]1 Sqrt[(s + \[Theta]c)/
Dc]))) -
2 Dp E^(Sqrt[s/Dp] (lh + x)) s Sqrt[s/
Dp] (E^((2 lh Sqrt[s])/Sqrt[
Dp]) (Sqrt[s] - Sqrt[Dp] \[Alpha]1) +
E^((2 Sqrt[s] x)/Sqrt[
Dp]) (Sqrt[s] + Sqrt[Dp] \[Alpha]1)) (s + \[Theta]c)^(
3/2) (-Dc s + Dp (s + \[Theta]c)) B2[1][s]))/(2 Dp^2 (s/Dp)^(
3/2) (Sqrt[s] + Sqrt[Dp] \[Alpha]1) (s + \[Theta]c)^(
3/2) (-Dc s + Dp (s + \[Theta]c)))
In[46]:= eq34 =
FullSimplify[
eq14, (lh - x) < 0 && s > 0 && s + \[Theta]c > 0 && Dp > 0 &&