Locating periodic solutions (dynamic systems)

Hey all,

I'm trying to program a MATLAB function that can locate periodic solutions of a dynamic system. My work is based on the *shooting method* (sorry if this is not the english term).

Consider the system , where and . The flow of the system is , and is a Poincare section / application / map (again, I'm sorry, I don't know the english word). Now, consider the following system:

Then, it is not hard to see that solving would give a *x* on a periodical solution of period *t* (obtained by solving the same system, of course). We can use Newton's method applied to the system to solve it numerically. Of course, we don't have DH (H's Jacobian matrix) explicitly since and are also not explicitly known. Note that we have:

Ok, finally we get to my problem! :) I can get around calculating 's Jacobian matrix as well as , but I'm having trouble figuring out how to get a good approximation of . Does anyone have an idea?

Thank you all for reading all this!