Solving ordinary differential equation

[thavamaran]

Hi guys! Attached is the differential equation that I want to solve both numerically and analytically, numerically done.

But analytically what method could I use? There is so many methods in differential equations. Please advice on this. Thank you very much

 

[HallsofIvy]

Your equation is
\frac{dZ}{dt}+ \Gamma\frac{I(t)}{qV}= Z(t)\left(\frac{-1+ \beta}{\tau_n}\right)


Are \Gamma, q, V, \beta, and \tau_n all constants? If so that is a "linear non-homogeneous first order differential equation" and there are a number of ways of solving it. For simplicity, I am going to write A for the constant
\frac{-1+ \beta}{\tau_n}
and B for
\frac{\Gamma}{qV}

Now the equation is
\frac{dZ}{dt}+ BI(t)= AZ
or, equivalently,
\frac{dZ}{dt}- AZ= -BI(t)

The simplest method is to find an "integrating factor". That is, find a function, u(t), so that the left side of that equation, multiplied by u(t), is simply
\frac{d(uZ)}{dt}.

Of course, by the chain rule,
\frac{d(uZ)}{dt}= u\frac{dZ}{dt}+ \frac{du}{dt}Z
while just multiplying by u would make the left side
u\frac{dZ}{dt}- AuZ

Comparing those, it is easy to see that we need to have
\frac{du}{dt}= -Au[/itex]<br /> which is the same as<br /> \frac{du}{u}= -Adt[/itex]&lt;br /&gt; &lt;br /&gt; Integrating both sides gives ln|u|= -At+ C or &lt;br /&gt; U(t)= e^{-At+ C}[/itex]&amp;lt;br /&amp;gt; &amp;lt;br /&amp;gt; Since we only need one function, take C= 0 so that u(t)= e^{-At}.&amp;lt;br /&amp;gt; &amp;lt;br /&amp;gt; That is, multiplying the entire equation by e^{-At} we get&amp;lt;br /&amp;gt; e^{-At}\frac{dZ}{dt}- Ae^{-at}Z= \frac{d(e^{-at}Z)}{dt}= -Be^{-At}I(t)[/itex]&amp;amp;lt;br /&amp;amp;gt; and so&amp;amp;lt;br /&amp;amp;gt; e^{-At}Z(t)= -B\int^t e^{-Ax}I(x)dx+ C&amp;amp;lt;br /&amp;amp;gt; or&amp;amp;lt;br /&amp;amp;gt; Z(t)= -Be^{At}\int^t e^{-Ax}I(x)dx+ Ce^{At}&amp;amp;lt;br /&amp;amp;gt; where I have changed the &amp;amp;amp;quot;dummy&amp;amp;amp;quot; variable in the integral to &amp;amp;amp;quot;x&amp;amp;amp;quot; so as not to confuse it with the actual variable t. In particular, the two exponentials, one inside the integral and one outside, do not cancel. The upper limit on the integral is to indicate that the variable becomes &amp;amp;amp;quot;t&amp;amp;amp;quot; after the integration and the &amp;amp;amp;quot;C&amp;amp;amp;quot; is the constant of integration. &amp;amp;lt;br /&amp;amp;gt; &amp;amp;lt;br /&amp;amp;gt; Of course, exactly what that integral is depends upon the function I(t) and, of course, you want to replace A and B with the original expressions.

 

[thavamaran]

Hi HallsofIvy, thank you very much and you are simply great to offer such help. But I still have trouble with solving the part as attached.

Im having trouble solving e^{-At}\frac{dZ}{dt}- Ae^{-at}Z= \frac{d(e^{-at}Z)}{dt}= -Be^{-At}I(t) where I am having trouble with integration by parts cause of the variables, its confusing me too much, can you explain to me a bit further on this, I am extremely sorry, but i have been trying to solve this since yesterday, deriving and reaching again to different types of false answers.

By the way, the equation is \frac{dZ}{dt}=\Gamma\frac{I(t)}{qV}+ Z(t)\left(\frac{-1+ \beta}{\tau_n}\right) and for integrating factor it should be:

\frac{dZ}{dt}-Z(t)\left(\frac{-1+ \beta}{\tau_n}\right)=\Gamma\frac{I(t)}{qV} .

thank you a millions times!