0
$\begingroup$

I am trying to plot slope field for $\ddot{a}(t)$ for this differential equation but I don't know how to

$\frac{d}{dt}\ddot{a}(t)=2\frac{\dot{a}(t)}{a(t)}\ddot{a}(t)$

I want to show that already accelerating solutions (with $\ddot{a}$ small and positive) are driven to accelerate even more, and already decelerating ones (with $\ddot{a}$ small and negative) are driven to further deceleration. So essentially I want to find out what $\dddot{a}(t)$ would look like with respect to time.

$\endgroup$
5
  • $\begingroup$ 4612 and 43478 might help. $\endgroup$
    – C. E.
    Jun 1, 2014 at 18:53
  • $\begingroup$ I tried using DSolve, NDSolve and VectorPlot but I'm not sure how to begin. I know how to plot first or second order differential equation but this one is different. $\endgroup$
    – user14701
    Jun 1, 2014 at 18:54
  • $\begingroup$ sol = ParametricNDSolve[{f'''[t] == 2 f'[t] f''[t]/f[t], f[1] == 1, f'[1] == 1, f''[1] == a}, f, {t, 1, 2}, a]; Plot[ Evaluate[Table[f'''[a][t] /. sol, {a, -2, 2, 1}]], {t, 1, 2}, PlotRange -> All] $\endgroup$ Jun 1, 2014 at 19:12
  • $\begingroup$ Are you sure that ther is a third derivative on the left? I think you have mistyped something... Are this a jerk parameter, related to scale factor? $\endgroup$
    – m0nhawk
    Jun 1, 2014 at 19:23
  • $\begingroup$ Yeap, there is a derivative on the left, this describes the rate of change of the second derivative of a. $\endgroup$
    – user14701
    Jun 1, 2014 at 20:33

1 Answer 1

2
$\begingroup$

A slope field is a two-dimensional image that represents the value of f(t,y) for values values of t and y for a first-order differential equation y'(t)=f(t,y(t)).

You have a third-order DE. You cannot plot a slope field.

You may be thinking of the phase portrait, which is a similar representation of the slope of an autonomous second order DE y''=f(y,y').

However, your equation is third order. The best you can do is to draw the three-dimensional phase-space for a third-order DE by converting the system into a system of first-order DEs:

$x = a(t)$

$y = a'(t)$

$z = a''(t)$

Your system becomes:

$x' = y$

$y' = z$

$z' = 2yz/x$

You can plot a unit vector in this direction in space using:

VectorPlot3D[
 {y, z, 2 y z/x}/Norm[{y, z, 2 y z/x}],
 {x, -3, 3}, {y, -3, 3}, {z, -3, 3},
 VectorStyle -> "Segment"
]

You can adjust the viewing window accordingly, I'm guessing [-3,3] is probably not what you want. You can change VectorStyle to something else, but the arrowheads often obscure the image. Options like this:

VectorStyle -> "Arrow3D", VectorScale -> {0.2, Scaled[0.1]}
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.