A Is it possible to solve this differential equation?

Donghyun Kim
Messages
1
Reaction score
0
Hello. I'm solving the second order nonlinear ODE, and I'm not sure that it's possible or not.

Please help.

A_0, B_0, C_0, a, b, c, e, p is all known constant.

dA/dt = -(a+b*C)*B
dB/dt = -p*c*A
dC/dt = -(1-p)*e*A

then I want to get a solution like A(t) = function of (A_0, B_0, C_0, a, b, c, e, p)
Can it be possible?

please help :(.
 
Physics news on Phys.org
Find C in terms of A and B from the first equation and then put it in the third equation. This gives you an equation involving only A and B. Then you can write one of them in terms of the other and put it in the second equation. Now you have an ordinary differential equation involving only one unknown function.
 
Donghyun Kim said:
Hello. I'm solving the second order nonlinear ODE, and I'm not sure that it's possible or not.

Please help.

A_0, B_0, C_0, a, b, c, e, p is all known constant.

dA/dt = -(a+b*C)*B
dB/dt = -p*c*A
dC/dt = -(1-p)*e*A

then I want to get a solution like A(t) = function of (A_0, B_0, C_0, a, b, c, e, p)
Can it be possible?

please help :(.

If pc = 0 then B is constant, and the system for A and C is then linear, so in principle you can solve it analytically. Similarly if e(1 -p) = 0 then C is constant, and the system for A and B is then linear, so again in principle you can solve it analytically.

Aside from those cases, the system is non-linear, so probably cannot be solved analytically and you may have to rely on numerical methods. The quantity X = e(1-p)B - pcC is conserved so the trajectory is confined to a plane of constant X. Thus you can eliminate C and consider the second-order system <br /> \dot A = -\left(a + \frac{b}{pc}(e(1-p)B - X)\right)B \\<br /> \dot B = -pcA.

EDIT: We have \ddot B = -pc\dot A = (pca + b(e(1-p)B - X)B which is a function of B; therefore one can multiply both sides by \dot B to obtain \frac{d}{dt}\left( \tfrac12 (\dot B)^2 - \tfrac13 be(1-p)B^3 + \tfrac12 bXB^2 - pcaB \right) = 0 which you can then integrate easily to obtain an equation for \dot B in terms of B (be careful with the choice of sign when taking the root). Unfortunately you probably can't then solve analytically for B, but if you could you would then have <br /> A(t) = A_0 - \int_0^t \left(a + \frac{b}{pc}(e(1-p)B(s) - X)\right)B(s)\,dt but you probably can't do that integral analytically even if you knew B.
 
Last edited:
Thread 'Direction Fields and Isoclines'
I sketched the isoclines for $$ m=-1,0,1,2 $$. Since both $$ \frac{dy}{dx} $$ and $$ D_{y} \frac{dy}{dx} $$ are continuous on the square region R defined by $$ -4\leq x \leq 4, -4 \leq y \leq 4 $$ the existence and uniqueness theorem guarantees that if we pick a point in the interior that lies on an isocline there will be a unique differentiable function (solution) passing through that point. I understand that a solution exists but I unsure how to actually sketch it. For example, consider a...
Back
Top