mbudisic

import numpy as np

def xyindependent(x,b,c,dt,randomfun):
    for idx in range(2,len(x)):
        x[idx] =(2-b*dt)*x[idx-1] + (b*dt-c*dt*dt-1)*x[idx-2] + dt*dt*randomfun()
    return x

import matplotlib.pyplot as plt

x0 = np.zeros(100,)

mbudisic