module_uparallel modified (and tested) for large x, now use Gaussians

The modification has been tested against Neufeld solution for the slab. I also compared the version with (Gaussians for x>10) and without the modif (use tables up to x=50) and they are similar.