Vary hardcoded integration params and plot output

[nbatta]

Reducing the max integration on y2d_tilde from 25Theta_c to 5Theta_c affects the variance at the 1e-4 level, we can speed things up in we choose 5*Theta_c without lose of accuracy.