... cov2(1,1)=sum(pw2x.*[(clustertot(1 ,:) -q(2,1)).* (clustertot(1 ,:) -q(2,1))])/sum(pw2x); cov2(1,2)=sum(pw2x.*[(clustertot(1 ,:) -q(2,1)).* (clustertot(2 ,:) -q(2,2))])/sum(pw2x); cov2(2,1)=sum(pw2x.*[(clustertot(2 ,:) -q(2,2)).* ... (clustertot(1 ,:) -q(2,1))])/sum(pw2x); cov2(2,2)=sum(pw2x.*[(clustertot(2 ,:) -q(2,2)).* (clustertot(2 ,:) -q(2,2))])/sum(pw2x); cov3(1,1)=sum(pw3x.*[(clustertot(1 ,:) -q(3,1)).* (clustertot(1 ,:) -q(3,1))])/sum(pw3x); ... cov3(1,2)=sum(pw3x.*[(clustertot(1 ,:) -q(3,1)).* (clustertot(2 ,:) -q(3,2))])/sum(pw3x); cov3(2,1)=sum(pw3x.*[(clustertot(2 ,:) -q(3,2)).* (clustertot(1 ,:) -q(3,1))])/sum(pw3x); cov3(2,2)=sum(pw3x.*[(clustertot(2 ,:) -q(3,2)).*...