Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- lam = 0.5; om = Pi; tau = 1;
- NDSolve[{x'[t] == lam x[t] + om y[t] - k[t] (x[t] - x[t - tau]),
- y'[t] == -om x[t] + lam y[t] - k[t] (y[t] - y[t - tau]),
- k'[t] == gam ((x[t] - x[t - tau]) (x[t] - 2 x[t - tau] +
- x[t - 2 tau]) + (y[t] - y[t - tau]) (y[t] - 2 y[t - tau] +
- y[t - 2 tau])), x[t /; t <= 0] == 0,
- y[t /; t <= 0] == 0, k[t /; t <= 2 tau] == 0}, {x, y,
- k}, {t, - tau, 40}]
Add Comment
Please, Sign In to add comment