k=.5;i=6;n=2;s=0;
forj=1:i-n,
h=k;b1=j;c=n;d=n;
cj=dima(h,b1,c,d);
h=k;b1=i-j;c=b1;d=n;
cst=dima(h,b1,c,d);
s=s+cj*cst;
end
cid=s*n/(i-n)
h=k;b1=i;c=i;d=n;
cis=dima(h,b1,c,d)
dif=cis-cid