for k = 1:1024
polys10(5,k) = p10(k)*polys10(4,k) - r10(3)*polys10(4,k);
end

 r10(2) = integre2(polys10(3,:),p10,p10)/integre2(polys(1,:),p10,p10);