Guest User

Untitled

a guest
Jun 13th, 2018
67
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.38 KB | None | 0 0
  1. lam = 0.5; om = Pi; tau = 1;
  2.  
  3. NDSolve[{x'[t] == lam x[t] + om y[t] - k[t] (x[t] - x[t - tau]),
  4. y'[t] == -om x[t] + lam y[t] - k[t] (y[t] - y[t - tau]),
  5. k'[t] == gam ((x[t] - x[t - tau]) (x[t] - 2 x[t - tau] +
  6. x[t - 2 tau]) + (y[t] - y[t - tau]) (y[t] - 2 y[t - tau] +
  7. y[t - 2 tau])), x[t /; t <= 0] == 0,
  8. y[t /; t <= 0] == 0, k[t /; t <= 2 tau] == 0}, {x, y,
  9. k}, {t, - tau, 40}]
Add Comment
Please, Sign In to add comment