This model is easier to build if polar coordinates are used. The same approach is used in the pendulum model and can be used to investigate Kepler's law.
|
|
Vensim Version
|
STELLA Version
|
Ang_Position(t) = Ang_Position(t - dt) + (Change_Position) * dt INIT Ang_Position = 0 Change_Position = Ang_velocity Ang_velocity = SQRT(Force/(mass*radius)) Force = 10 mass = 2 radius = 2 x_position = radius*COS(Ang_Position) y_position = radius*SIN(Ang_Position)
Time Specs
Standard: Range: 0 - 12, dT = 0.25, Integration Method = Euler's
