Function to give probability of observing x given age a and the viral kinetics curve

p_a(x, a, pars, viral_loads, sd_mod)