function[vec] = P(n)
if n == 0
    vec = 0;
    return
else
    vec = [(1/3*P(n-1)),(2/3 + 1/3*P(n-1))];
end
end


