Computes the behaviour of a catenary,
as an inextensible cable, via the Nelder-Mead algorithm. For details,
see the plate by Newton-Raphson.
The method solves the original 4 ODE's.
The graph shows: (a) y (catenary)
vs. x;
or (b) θ (green) and q (red) vs. s.
† The coordinates of B must lie in a unit circle
and the tension, q, is given as related to the cable weight
(due to the non-dimensional formulation).
Other proposed data: (.35, -.35; -85, .1; etc.).
(Problem suggested and managed by Prof. A. Pinto da Costa,
IST.) |
