A pendulum of length swings through arbitrary angles (not necessarily small) and…
PHYSICS · Subjective · Fermi University Physics · Mathematical Methods for Physics · Computational Physics · Simulating classical dynamical systems
Computational Physics, Ch Stepping the equations of motion: four integrators measured
Problem
A pendulum of length L swings through arbitrary angles (not necessarily small) and experiences a drag torque proportional to its angular velocity, so its equation of motion is
θ¨=−ω02sinθ−γθ˙,ω02=g/L,
where θ is the angle from the vertical, γ is a constant damping rate, and g=9.8ms−2.
(a) Write this as a first-order system: state the two-component state vector y and the function f(y) such that y˙=f(y).
(b) Take L=0.80m and γ=0.15s−1. Evaluate f numerically at the state θ=0.30rad, θ˙=0rads−1.