def moving_average(a, n=3):

    b = np.empty(len(a))
    for i in range(len(a)):
        if i < n:
            b[i] = np.sum(a[:i+1])/(i+1)
        else:
            b[i] = np.sum(a[i-n:i+1])/(n+1)
    return b