changed the parameters for the subspace simulation, the and the results is much better. However, the cubatured method are still not good.