- #1
BeeKay
- 16
- 0
Homework Statement
Hello,
I have a question about using Eulers Method to approximate a solution to a differential equation. The problem lists forces that would be applied on an object and influences its velocity and therefore its position. I believe I am doing the Euler method correct to approximate what the "new velocity" would be, but based on somewhat wonky results I want to confirm if I am making an error using the Euler Method ( or in the Fortran code, but I feel like that is fine).
Homework Equations
[itex] m \frac{d^ {2}y}{dt^2} + b \frac{dy} {dt} + ky = F_{o} sin(νt), y(0) = y_{o}, \frac{dy} {dt} = ˙y_{o}. [/itex] when t=0
[itex] \frac{dy} {dt} = v, y(0) = y_{o}, [/itex]
[itex] \frac{dv} {dt} = \frac{− k} {m} y − \frac {b} {m} v + \frac {F_{o}} {m} sin(νt), v(0) = ˙y_{o} [/itex]
[itex] \frac {y_{n+1} − y_{n}} {∆t} = v_{n}, [/itex]
[itex] \frac{v_{n+1} − v_{n}} {∆t} = \frac {− k }{m} y_{n} − \frac {b} {m} v_{n} + \frac {Fo} {m} sin(νt_{n}), [/itex]
[itex] t_{n+1} = t_{n + ∆t} [/itex]
The Attempt at a Solution
The way I interpreted it was to first find [itex] \frac{v_{n+1} − v_{n}} {∆t} [/itex] by find it at the initial values. If I use 0.02 as the timestep, I solved by substituting in the initial given y value and v value. Some other parameters were given, but they would not be changing. I multiplied this entire solution by the time step to find how much the velocity changed in that time. I then added this change to the original velocity. From there, I am a little confused as to how it applies to the position estimation. If [itex] \frac {y_{n+1} − y_{n}} {∆t} = v_{n}, [/itex], then shouldn't I just take the solution to [itex] \frac{v_{n+1} − v_{n}} {∆t} = \frac {− k }{m} y_{n} − \frac {b} {m} v_{n} + \frac {Fo} {m} sin(νt_{n}), [/itex] multiply it by the timestep, and then add it to the previous y value? Thanks in advance and sorry if I messed up LaTeX. If you need more information about the problem or my process, let me know.