The paper discusses new cubature formulas for the Riesz potential and the fractional Laplacian (−∆) α/2 , 0 < α < 2, in the framework of the method approximate approximations. This approach, combined with separated representations, makes the method successful also in high dimensions. We prove error estimates and report on numerical results illustrating that our formulas are accurate and provide the predicted convergence rate 2, 4, 6, 8 up to dimension 10 4 .