f(u) = u^2/(1+u^2)
param alpha=1, beta=1, gamma=1, tau=1
x' = -x + alpha*f(x) - beta*f(y) + beta
y' = (-y + gamma*f(x))/tau
@ MAXSTOR=1000000
done