How could i use 3 random variables instead of one with shape=3 with pymc3

You can stack your variables with pm.math.stack([a, b, c]). Note however that unless you need to use different prior distributions, it’s better to use a single distribution with shape=n, as that results in more efficient code.