Here is what the mathcad output looks like.
What it does is
1) take the inital values of qo and po and plugs them into p and q for both equations pj+1 and qj+1
2) it then returns a new value p1 and q1 after the first iteration
3) it then plugs in p1 and q1 into p and q again and returns p2 and q2
4) it continues to do this for the number of iterations set in j
5) I can then easily graph my results with respect to the number of iterations.
To me it seems Mathematica has a lot of trouble bouncing values back and forth since everything is so "top down"
Attachment: mathcaditerate.tiff, URL: ,