function res=number_pts(r,start,inc)
res = 1
for j = start:inc:(length( r ) - 1) 
res = res * r( j );
end