LotkaVolterra
As originally coded by Ignasi Alvarez Garriga (I think...)
; Lotka Volterra model for predator-prey relations dynamics
;
; dx/dt = Ax-Bxy
; dy/dt= -Cy+Dxy
;
; parameters : A, B, C, D
;
; Ignasi Alvarez, 1998
; Rev. 2005
izoom = p4
iamp = p5
ipanL = sqrt(p6)
ipanR = sqrt(1-p6)
iparammod = p7
ka = <paramA>
kb = <paramB>
kc = <paramC>
kd = <paramD>
kdt = <dt>
kdt2 = <dt>/2
kgainx = <gainX>
kgainy = <gainY>
ky init 1
kx init 1
kxo init 0
kyo init 0
kynew init 0
kxnew init 0
kynewo init 0
kxnewo init 0
kcont init 0
ka = ka+iparammod
kb = kb+iparammod
kc = kc+iparammod
kd = kd+iparammod
loop:
;sistema d'equacions
kxnew=ka*kx-kb*kx*ky
kynew=(-1)*kc*ky+kd*kx*ky
kxo=kx
kyo=ky
kxnewo=kxnew
kynewo=kynew
kx=kxo+kxnewo*kdt
ky=kyo+kynewo*kdt
;sistema d'equacions
kxnew=ka*kx-kb*kx*ky
kynew=(-1)*kc*ky+kd*kx*ky
kx=kxo+kdt2*(kxnew+kxnewo)
ky=kyo+kdt2*(kynew+kynewo)
kcont=kcont+1
if kcont!=izoom kgoto loop
kcont=0
asig2 = kx*kgainx+ky*kgainy
asig1 dcblock asig2*iamp
; Optionally output to a Display object:
;outvalue "kdisplax", 0.25+kx*kx/25
;outvalue "kdisplay", 0.25+ky/25
blueMixerOut asig1*ipanL,asig1*ipanR
paramA
61
34
true
0.01
10.0
1.1089
80
true
79
135
paramB
154
34
true
0.01
10.0
0.55945
80
true
172
135
paramC
245
34
true
0.01
10.0
1.154854
80
true
263
135
paramD
334
34
true
0.01
10.0
3.675331
80
true
352
135
dt
61
164
true
2.0E-4
0.36
0.193269
80
true
93
266
gainX
152
164
true
-2000.0
2000.0
1164.8
80
true
173
264
gainY
242
164
true
-2000.0
2000.0
521.60004
80
true
263
264