On 5/6/2009 10:00 AM Talbot, Gerry apparently wrote: > for n in xrange(1,N): > y[n] = A*x[n] + B*y[n-1] So, x is known before you start? How big is N? Also, is y.shape (N,)? Do you need all of y or only y[N]? Alan Isaac