faster way of computing the Roland_Sylvain integrals by pre-computing many common terms