Assumed the following example of double pendulum by Thison (Answered by Jens):
Imagine that in the given animation(i.e. http://i.stack.imgur.com/p4TgL.gif ), the second pendulum (Second Linkage) can not rotate 360 degree. In other words, I want to somehow put a constraint on this system with these given equations and initial conditions and limit their motions. For instance, Second pendulum can have a rotation like 0<phi2[t]<180
.
I wonder how it is possible to apply this constraint in NDSolve[]
or any possible solution for this problem.
Thank you in advance.