A0=[1 0 0 ; 2/5 2/5 1/5 ; 2/5 1/5 2/5];
A1=[2/5 2/5 1/5; 0 1 0; 1/5 2/5 2/5];
A2=[2/5 1/5 2/5; 1/5 2/5 2/5 ; 0 0 1];
 
 
D0=[1 2/9 2/9 ; 0 2/9 -1/9; 0 -1/9 2/9];
D1=[2/9 0 -1/9; 2/9 1 2/9; -1/9 0 2/9];
D2=[2/9 -1/9 0 ; -1/9 2/9 0; 2/9 2/9 1];
 
q0=[1/2 1]; q1=[0 0]; q2=[1 0];
 
E0=1/125.*[47 -3 -3; 14 9 -6; 14 -6 9];
E1=1/125.*[9 14 -6; -3 47 -3; -6 14 9];
E2=1/125.*[9 -6 14; -6 9 14; -3 -3 47];