Hi,
using the constraint programming engine could be an option?
using CP;
int scale=10;
int D=1;
float d=0.1;
int C=30; execute
{ cp.param.timelimit=10;
} dvar
int xs2 in -1000*scale..1000*scale; dvar
int ys2 in -1000*scale..1000*scale; dexpr
float xs=xs2/scale;; dexpr
float ys=ys2/scale;
int x[1..3]=[0,10,10];
int y[1..3]=[0,10,0];
int P1=1000;
int P2=500; maximize P2*(1+d)^(-sum (i in 1..3)((xs-x[i])^2+(ys-y[i])^2)^0.5)+P2*(1+d)^(-sum (i in 1..2)((xs-x[i])^2+(ys-y[i])^2)^0.5) - C*(1+d)^(-sum (i in 1..3)((xs-x[i])^2+(ys-y[i])^2)^0.5) ; subject to
{ 0<= xs <=10; 0<= ys <=10;
}
regards
#DecisionOptimization#MathematicalProgramming-General