In[]:=
funx=-1/10+5Log[a1x+a1y-14]-2+95/100(-110(1-ix/40)+5Log[a2x+a2y-14])funy=-1/10+5Log[a1x+a1y-14]-2+95/100(-110(1-ix/40)+100Log[a2x+a2y-14])dWx/da1x=-2/10a1x+5/(a1x+a1y-14)0dWy/da1y=-2/10a1y+5/(a1x+a1y-14)0dWx/da2x=95/100(-2/10a2x(1-ix/40)+5/(a2x+a2y-14))0dWy/da2y=95/100(-2/10a2y(1-ix/40)+5/(a2x+a2y-14))0dWx/dix=-+(95/100)(1/400)0dW/diy=-0Reduce95100(-2/10a2x(1-ix/40)+5/(a2x+a2y-14))0,95100(-2/10a2y(1-ix/40)+5/(a2x+a2y-14))0,+(95/100)(1/400)0,a2x≥0,a2y≥0,0<=ix≤24,0<=iy≤24,{a2x,a2y,ix,iy}
2
a1x
1/2
ix
2
a2x
2
a1y
1/2
iy
2
a2y
-1/2
ix
2
a2x
-1/2
iy
-1/2
ix
2
a2x
Out[]=
--2-1-+5Log[-14+a2x+a2y]
2
a1x
10
ix
+5Log[-14+a1x+a1y]+19
20
1
10
2
a2x
ix
40
Out[]=
--2-1-+100Log[-14+a2x+a2y]
2
a1y
10
iy
+5Log[-14+a1x+a1y]+19
20
1
10
2
a2y
ix
40
dWx
da1x
Out[]=
-+0
a1x
5
5
-14+a1x+a1y
dWy
da1y
Out[]=
-+0
a1y
5
5
-14+a1x+a1y
dWx
da2x
Out[]=
19
20
5
-14+a2x+a2y
1
5
ix
40
dWy
da2y
Out[]=
19
20
5
-14+a2x+a2y
1
5
ix
40
dWx
dix
Out[]=
19
2
a2x
8000
1
ix
dW
diy
Out[]=
-0
1
iy
Out[]=
False
