read DOUBLE_SUM_SINGLE_INTEGRAL : P:=1: F:=rf(al+1,n)*rf(be+1,n)/n!/rf(al+be+1,n)/t^(n+1)* (1+t)^(-al-be-1)*rf(al+be+1,2*(m+n1))/n1!/m!/rf(al+1,n1)/rf(be+1,m)* (t/(4*(1+t)^2))^(m+n1)* (1+x)^m*(1+y)^m*(1-x)^n1*(1-y)^n1; ORDER:=2: resh:=[al,be,x,y,n]: S1:=1/(t+1): S2:=(n1+al)*n1: S3:=m*(m+be): disorcon:=continuous: cer:=findope(P,F,x,t,n1,m,ORDER,resh,D_x,S1,S2,S3,disorcon); quit: #-------------end typical input file--------------------