Hello, Can anybody explain to me which of these results returned by DSolve is correct and why? If I try
resA=DSolve[u''[x]-(s/D)*u[x]==-(k/D)*(u0/s)*(1-Exp[-Sqrt[s/D]*x]),u[x],x];
resB=FullSimplify[resA];
Export["res.txt",resB];
I get "Complex infinity" as a result. But if I simplify the above, by setting a=s/D and b=(k/D)*(u0/s), and try
resA=DSolve[u''[x]-a*u[x]==-b*(1-Exp[-Sqrt[a]*x]),u[x],x];
resB=FullSimplify[resA];
Export["res.txt",resB];
I get
{u[x] -> b/a + E^(Sqrt[a]*x)*C[1] - (b + 2*Sqrt[a]*b*x - 4*a*C[2])/(4*a*E^(Sqrt[a]*x))}
which is clearly a different result.
Leslaw