function [q,p]=dima(h,b1,c1,d,e)
a=1-h;z=b1/a;w=z:-1:z-b1+d*e+1;
t=b1-d*e:-1:1;
s=w./t;
q=(d/c1)*prod(s);
p=0;