Iterative solution difficulty

*To*: mathgroup at smc.vnet.net*Subject*: [mg99942] Iterative solution difficulty*From*: James <jborrell at umd.edu>*Date*: Tue, 19 May 2009 07:17:00 -0400 (EDT)

I am trying to solve a system of ODEs using NDSolve until a desired ending condition is reached, then the value at the final time x[tf-] is "linearly" related to the new initial conditions, x[tf+] and the ODEs are solved until the ending condition again and so on. The following code is the numerical solution to the equation of motion of a simple pendulum. Physically, the pendulum is released from horizontal and then stopped at the vertical and returned to horizontal and solved again. I expect the output to look periodic, similar to a capacitor discharging. ICx1 = Pi/2; ICxprime1 = 0; output = NDSolve[{x''[t] == -g/L Sin[x[t]], x[0] == ICx1, x'[0] == ICxprime1} /. {g ->9.8, L -> 1}, x, {t, 0, 10},Method -> {EventLocator, "Event" -> x[t] > 0,EventAction :> Throw[t1 = t; ics1 = x[t];,StopIntegration"]}] x1 = x[t1] + Pi/2 /. output[[1]] x1p = x'[t1] + 5 /. output[[1]]