q3=(3a3-3a1a2+a13)/3,