function v = rigaP(i)
q = i/(i+1);
v = q.^linspace(10,20,11);
v = v/sum(v);
