function V_0=initial_a(d,k)
    V_0=zeros(k,2)
    for j=1:k
    V_0(j,1)=d^(j-1)/(d^k-1)
    end
end 
