function res = pk(k,v )
if k == 1
    res = v(1);
    return
elseif k==2
    res = v(1)*v(2)-1;
    return
else 
    res = v(k)*pk(k-1,v) - pk(k-2,v);
end
end

