代码如下
h=Symbol('h')
i=Symbol('i')
solved_value=solve([0.396*(800-676+20*h-100*h-2*i)**2+0.594*(600-676+20*h+100*h-2*i)**2+0.004*(400-100*h+199*i-676-20*h)**2+0.006*(200+100*h+199*i-676-20*h)**2==3600,0.396*(800-676-20*h-100*h-2*i)**2+0.594*(600-676-20*h+100*h-2*i)**2+0.004*(400-100*h+199*i-676-20*h)**2+0.006*(200+100*h+199*i-676-20*h)**2==3600],[h,i])
print(solved_value)