a1=(1+a,1,1,1)^T,a2=(2,2+a,2,2)^T,a3=(3,3,3+a,3)^T,a1=(4,4,4