import numpy as np
import matplotlib.pyplot as plt
from matplotlib import animation
dt = 0.01
tfinal = 1
x0 = 0
sqrtdt = np.sqrt(dt)
n = int(tfinal/dt)
xtraj = np.zeros(n+1, float)
trange = np.linspace(start=0,stop=tfinal ,num=n+1)
xtraj[0] = x0
for i in range(n):
xtraj[i+1] = xtraj[i] + np.random.normal()
x = trange
y = xtraj
fig, ax = plt.subplots(1, 1, figsize = (6, 6))
def animate(i):
ax.cla()
ax.plot(x[:i], y[:i])
ax.set_xlim([x0, tfinal])
ax.set_ylim([1.1*np.min(y), 1.1*np.max(y)])
anim = animation.FuncAnimation(fig, animate, frames = len(x) + 1, interval = 1, blit = False)
plt.show()