i=6;j=2;
q=[150 399 261 56 13 6;64 162 130 39 11 7];r=[214 561 391 95 24 13];
qt=[885 413];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.*r2;
a=[h;g];d=q-a;
d2=d.^2;
f=d2./a;s=sum(f);
xx=sum(s)