powered by
dX_t = (1.5X_t-0.4X_tY_t)dt+sqrt(0.05X_t)dW_t dY_t = (-1.5Y_t+0.4X_tY_t-0.2Y_t^2)dt+sqrt(0.1X_t)dB_t
where dW_t and dB_t are standard Brownian motions, t is time and X_0 = 5, Y_0 = 5.
dW_t
dB_t
t
X_0 = 5
Y_0 = 5
data("SDEsim3")
time
time[i+1]-time[i] = 1/4
mx
my