ClearAll[Q1,Q2,Q3,Q4];​​mndmd=500;​​stdmd=150;​​pri=100;​​salpri=10;​​srtcst=20;​​hol=0.1;​​bcrodr=1;​​M=1000;​​​​c1=42;​​c2=44;​​c3=40;​​c4=46;​​​​yldmx1=0.96;yldmn1=0.92;​​yldmx2=0.94;yldmn2=0.88;​​yldmx3=0.90;yldmn3=0.88;​​yldmx4=0.97;yldmn4=0.94;​​​​timmn1=0;timstd1=3;aa1=5;​​timmn2=0;timstd2=3;aa2=-2;​​timmn3=0;timstd3=5;aa3=2;​​timmn4=0;timstd4=5;aa4=-1;​​​​​​dmnfun=
1
stdmd
2π
-1
2
2
x-mndmd
stdmd

;​​​​yldfun1=
1
yldmx1-yldmn1
;​​yldfun2=
1
yldmx2-yldmn2
;​​yldfun3=
1
yldmx3-yldmn3
;​​yldfun4=
1
yldmx4-yldmn4
;​​​​tinfun1=
2
timstd1
2π
-1
2
2
l1-timmn1
timstd1

aa1
l1-timmn1
timstd1
∫
-∞
1
2π
-
2
t1
2

t1;​​tinfun2=
2
timstd2
2π
-1
2
2
l2-timmn2
timstd2

aa2
l2-timmn2
timstd2
∫
-∞
1
2π
-
2
t2
2

t2;​​tinfun3=
2
timstd3
2π
-1
2
2
l3-timmn3
timstd3

aa3
l3-timmn3
timstd3
∫
-∞
1
2π
-
2
t3
2

t3;​​tinfun4=
2
timstd4
2π
-1
2
2
l4-timmn4
timstd4

aa4
l4-timmn4
timstd4
∫
-∞
1
2π
-
2
t4
2

t4;​​​​objfun[Q1_,Q2_,Q3_,Q4_,y1_,y2_,y3_,y4_]:=N​​
yldmx4
∫
yldmn4

yldmx3
∫
yldmn3

yldmx2
∫
yldmn2

yldmx1
∫
yldmn1
​​
(y1*k1*Q1+y2*k2*Q2+y3*k3*Q3+y4*k4*Q4)
∫
0
((pri*x+salpri*((y1*k1*Q1+y2*k2*Q2+y3*k3*Q3+y4*k4*Q4)-x))dmnfun)x​​​​+
∞
∫
(y1*k1*Q1+y2*k2*Q2+y3*k3*Q3+y4*k4*Q4)
((pri*(y1*k1*Q1+y2*k2*Q2+y3*k3*Q3+y4*k4*Q4)-srtcst*(x-(y1*k1*Q1+y2*k2*Q2+y3*k3*Q3+y4*k4*Q4)))dmnfun)x*yldfun1k1*yldfun2k2*yldfun3k3*yldfun4k4+​​​​
yldmx1
∫
yldmn1
hol*
0
∫
-∞
(y1*k1*Q1*l1*tinfun1)l1​​-bcrodr*
∞
∫
0
(y1*k1*Q1*l1*tinfun1)l1​​-c1*y1*k1*Q1yldfun1k1+​​​​
yldmx2
∫
yldmn2
hol*
0
∫
-∞
(y2*k2*Q2*l2*tinfun2)l2​​-bcrodr*
∞
∫
0
(y2*k2*Q2*l2*tinfun2)l2​​-c2*y2*k2*Q2yldfun2k2+​​​​
yldmx3
∫
yldmn3
hol*
0
∫
-∞
(y3*k3*Q3*l3*tinfun3)l3​​-bcrodr*
∞
∫
0
(y3*k3*Q3*l3*tinfun3)l3​​-c3*y3*k3*Q3yldfun3k3+​​​​
yldmx4
∫
yldmn4
hol*
0
∫
-∞
(y4*k4*Q4*l4*tinfun4)l4​​-bcrodr*
∞
∫
0
(y4*k4*Q4*l4*tinfun4)l4​​-c4*y4*k4*Q4yldfun4k4​​​​​​objfun[150,150,150,150,1,1,1,1](*testoffunction.....working*)​​​​NMaximize[{objfun[Q1,Q2,Q3,Q4,y1,y2,y3,y4],​​{Q1≥0.000001,Q2≥0.000001,Q3≥0.000001,Q4≥0.000001,Q1≤y1*M,Q2≤y2*M,Q3≤y3*M,Q4≤y4*M,0≤y1≤1,0≤y2≤1,0≤y3≤1,0≤y4≤1,y1+y2+y3+y44,Element[{y1,y2,y3,y4},Integers]}},{Q1,Q2,Q3,Q4,y1,y2,y3,y4},Method{"RandomSearch","SearchPoints"100,MaxIterations10000000000000000,Tolerance0.00000001}](*Thiscalculationisgeneratingtheerror*)​​
Out[]=
21635.3
​