Mercurial > hg > octave-nkf

view examples/oregonator.m @ 3112:92394a9e4784

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1997-11-27 02:02:27 by jwe]
author jwe
date Thu, 27 Nov 1997 02:02:32 +0000
parents 8c7955a8d49f
children 7c96e85c76db
line wrap: on
line source

## The `oregonator'.

function dx = oregonator (x, t)

  dx = zeros (3, 1);

  dx(1) = 77.27*(x(2) - x(1)*x(2) + x(1) - 8.375e-06*x(1)^2);
  dx(2) = (x(3) - x(1)*x(2) - x(2)) / 77.27;
  dx(3) = 0.161*(x(1) - x(3));

end