;k=.5; g=50;a=(1+k*(g-1))/g;i=150;n=3;
h=k;b1=i-1;c1=b1;d=n;

e=1;s=dima(h,b1,c1,d,e)/(1-k);
for j=2:i-n;
h=k;b1=j;c=1;d=1;e=0;
cc(j)=dima(h,b1,c,d,e);
h=k;b1=i-j;c=b1;d=n;e=1;
cst(j)=dima(h,b1,c,d,e);
s=s+cc(j)*cst(j);
end
cid=s*n/(i-n)
h=k;b1=i;c=i;d=n;e=1;
cis=dima(h,b1,c,d,e)
dif=cis-cid
v=(a^(k/(1-k))*(1-a));pin=cis*(a/(1-a))^n*v^i