Issue vectorizing a recursive function that is used in iterative scheme to calculate Numpy array

Question

I have the following recursive function,

def subspaceiterate(A,V,v,j):
    if j == 0:
        return v
    else:
        v_jm1 = V[:,j-1]
        v_jm1 = np.r