function p = RKV_test(b, A, c)

ok = 1;
p = 0;

while ok,
    p = p+1;
    for l = 1:p,
        for k = 0:p-l,
            soll = 1 ./ prod( l:(l+k) );
            ok = ok & abs(b' * A^k * (c.^(l-1)) - soll) < eps;
            %ok = ok & ( b' * A^k * (c.^(l-1)) == soll);
        end
    end
end
p = p-1;