3.3b: New routine, plots intensity for various tau for fixed temperature as function of u.