mx1=[0          0                       0                           vrs(4)*vrs(19);
     vrs(1)     vrs(2)*(1-vrs(5))       0                           0;
     0          vrs(2)*vrs(5)           vrs(3)*(1-vrs(6))           0;
     0          0                       vrs(3)*vrs(6)               vrs(4)];                 
    
 mx2=[0          0                       0                           vrs(10)*vrs(19);
     vrs(7)     vrs(8)*(1-vrs(11))       0                           0;
     0          vrs(8)*vrs(11)           vrs(9)*(1-vrs(12))           0;
     0          0                       vrs(9)*vrs(12)               vrs(10)];
 
 mx3=[0          0                       0                           vrs(16)*vrs(19);
     vrs(13)     vrs(14)*(1-vrs(17))       0                           0;
     0          vrs(14)*vrs(17)           vrs(15)*(1-vrs(18))           0;
     0          0                       vrs(15)*vrs(18)               vrs(16)];
 
 mx=zeros(dims,dims);
 mx(1:4,1:4) = mx1;
 mx(5:8,5:8) = mx2; 
 mx(9:12,9:12) = mx3;