求lingo高手帮我看下我的lingo程序到底哪里出错啦

廖文强

是2002年的国赛彩票问题，用lingo8计算了5个小时都没有出结果，用lingo9以上的版本就总是报错。。。跪求高手解答！！跪谢

Model:
Sets:
Ways/1..7/:p,x,ux;
ab/1..6/:a,b;
Price/1..3/:r;
Endsets
Max=@sum(ways:p*ux);
data:
       a,b=10,233,4,54,3,17,4,20,2,10,2,10;
enddata

p(1)=@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/@exp(@lgm(n+1));
p(2)=m*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/@exp(@lgm(n+1));
p(3)=m*(n-m-1)*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/@exp(@lgm(n+1));
p(4)=0.5*m*(m-1)*(n-m-1)*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/@exp(@lgm(n+1));
p(5)=0.25*m*(m-1)*(n-m-1)*(n-m-2)*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/@exp(@lgm(n+1));
p(6)=m*(m-1)*(m-2)*(n-m-1)*(n-m-2)*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/(12*@exp(@lgm(n+1)));
p(7)=m*(m-1)*(m-2)*(n-m-1)*(n-m-2)*(n-m-3)*@exp(@lgm(m+1))*@exp(@lgm(n-m+1))/(36*@exp(@lgm(n+1)));

@for(ways(i)|i#le#3:x(i)=(1-p(4)*x(4)+p(5)*x(5)+p(6)*x(6)+p(7)*x(7))*r(i)/p(i));
@for(ways(i):ux(i)=1-@exp(-(x(i)/ 630589)^2));
r(1)+r(2)+r(3)=1;
@bnd(0.5,r(1),0.8);
@bnd(600000,x(1),5000000);
@for(ways(i)|i#le#6:a(i)<(x(i)/x(i+1)));
@for(ways(i)|i#le#6: b(i)>(x(i)/x(i+1)));
@for(ways(i)|i#le#6:p(i)<p(i+1));
@bnd(5,m,7);
@bnd(29,n,60);
@for(ways(i)bnd(0,p(i),1));
@for(ways(i):x(i)>0);
@for(price(j):r(j)>0);
@gin(m);@gin(n);
end

廖文强

这个是题目

廖文强

这个是概率P