Simple pendulum Runge-Kutta-Fehlberg solver without error analysis

Simple pendulum problem solver with error analysis

This webpage uses the Runge-Kutta-Fehlberg fourth-order method with fifth-order error checking (RKF45) to approximate the solution to the problem of the simple pendulum: \[ \dfrac{d^2 \theta}{dt^2} = -\dfrac{g}{l} \cos{\theta} \] where \(g\) is the acceleration due to gravity in metres per second squared, \(l\) is the length of the pendulum in metres, \(\theta\) is the angle from the positive \(x\) axis (in radians) and \(t\) is the time in seconds. Below you can specify the various parameters for the problem we will solve. From this equation, the period \(T\) of the problem is approximated when the conditions for periodicity are satisfied, namely using the equation: \[ \begin{aligned} T &= 2 \left|\int_{\theta_\mathrm{min}}^{\theta_\mathrm{max}} \dfrac{d\theta}{\sqrt{\dot{\theta}_0^2 + \dfrac{2g}{l} (\sin{\theta_0} - \sin{\theta}})}\right| \end{aligned} \] where \(\theta_\mathrm{min}\) and \(\theta_\mathrm{max}\) are the two closest values for which \(\dot{\theta} = 0\). \(T\) is calculated using Chebyshev-Gauss quadrature (Simpson's rule could not be used as there are unremovable singularities at the endpoints which makes Simpson's rule markedly less accurate).

Parameter	Value	Explanation
\(g\):		Acceleration due to gravity in \(\mathrm{m} \cdot \mathrm{s}^{-2}\).
\(l\):		Length of the pendulum rod in metres (m).
\(t_0\):		Starting time for the simulation in seconds (s).
\(t_f\):		End time for the simulation in seconds. The default \(t_f - t_0\) value is equal to four times the period \(T\) of the problem.
\(\theta_0\):		Angle from the positive \(x\)-axis in radians at time \(t_0\).
\(\dot{\theta}_0\):		First derivative of \(\theta\) with respect to \(t\) at \(t_0\) (in \(\mathrm{radians}\cdot \mathrm{s}^{-1}\)).
\(N\):		Number of Chebyshev-Gauss quadrature nodes used to approximate the period of the problem. \(N \gt 1e7\) usually freezes up the page.
\(\epsilon\):		Error tolerance in both \(\theta\) and \(\dot{\theta}\). \(\epsilon \lt \)2e-12 often freezes the webpage.
dtInitial:		Initial guess for step size.

Parameter

Value

Explanation

\(g\):

Acceleration due to gravity in \(\mathrm{m} \cdot \mathrm{s}^{-2}\).

\(l\):

Length of the pendulum rod in metres (m).

\(t_0\):

Starting time for the simulation in seconds (s).

\(t_f\):

End time for the simulation in seconds. The default \(t_f - t_0\) value is equal to four times the period \(T\) of the problem.

\(\theta_0\):

Angle from the positive \(x\)-axis in radians at time \(t_0\).

\(\dot{\theta}_0\):

First derivative of \(\theta\) with respect to \(t\) at \(t_0\) (in \(\mathrm{radians}\cdot \mathrm{s}^{-1}\)).

\(N\):

Number of Chebyshev-Gauss quadrature nodes used to approximate the period of the problem. \(N \gt 1e7\) usually freezes up the page.

\(\epsilon\):

Error tolerance in both \(\theta\) and \(\dot{\theta}\). \(\epsilon \lt \)2e-12 often freezes the webpage.

dtInitial:

Initial guess for step size.

Information	Data	Notes
\(T\):		Period of the problem, auto-calculated.
\(\theta_\mathrm{min}\)		The value of \(\theta_\mathrm{min}\) used in the above period equation.
\(\theta_\mathrm{max}\)		The value of \(\theta_\mathrm{max}\) used in the above period equation.
RKF45 step number		Number of steps used in our Runge-Kutta-Fehlberg approximation of the solution.

Information

Data

Notes

\(T\):

Period of the problem, auto-calculated.

\(\theta_\mathrm{min}\)

The value of \(\theta_\mathrm{min}\) used in the above period equation.

\(\theta_\mathrm{max}\)

The value of \(\theta_\mathrm{max}\) used in the above period equation.

RKF45 step number

Number of steps used in our Runge-Kutta-Fehlberg approximation of the solution.