i=6;j=2;
q=[190 490 346 82 19 12;24 71 45 13 5 1];r=[214 561 391 95 24 13];
qt=[1139 159];n=1298;
r1=(1/n).*r;l=ones(size(r1));
%r2=l-r1;
q1=(1/n).*qt;
h=q1(1)*n.*r1;
g=q1(2)*n.*r1;
a=[h;g];d=q-a;
d2=d.^2;
f=d2./a;s=sum(f);
xx=sum(s)