- #1
KevinL
- 37
- 0
I have no programming experience and trying to get mathematica to do what I need it to do is frustrating. I have the following functions that I need to iterate. For notational purposes, k[t+1] is the value of K in the next period. w is a parameter.
k[t+1] = -x[t] - y[t] + w
x[t+1] = k[t]*sqrt(y[t]/(-y[t] + w)) - y[t]
y[t+1] = k[t]*sqrt(x[t]/(-x[t] + w)) - x[t]
As I care more about the behavior of x and y, this system can be reduced to:
x[t+2] = (-x[t]-y[t]+w)*sqrt(y[t+1]/(-y[t+1] + w)) - y[t+1]
y[t+2] = (-x[t]-y[t]+w)*sqrt(x[t+1]/(-x[t+1] + w)) - x[t+1]
I have already found for what w the system becomes unstable by taking the Jacobian of the 3-system and seeing where it exceeds unity in absolute value. Now I just want to iterate the system and get some pretty graphs.
k[t+1] = -x[t] - y[t] + w
x[t+1] = k[t]*sqrt(y[t]/(-y[t] + w)) - y[t]
y[t+1] = k[t]*sqrt(x[t]/(-x[t] + w)) - x[t]
As I care more about the behavior of x and y, this system can be reduced to:
x[t+2] = (-x[t]-y[t]+w)*sqrt(y[t+1]/(-y[t+1] + w)) - y[t+1]
y[t+2] = (-x[t]-y[t]+w)*sqrt(x[t+1]/(-x[t+1] + w)) - x[t+1]
I have already found for what w the system becomes unstable by taking the Jacobian of the 3-system and seeing where it exceeds unity in absolute value. Now I just want to iterate the system and get some pretty graphs.