Question: how to fix bug on DirectSearch package?

i'm doing an optimaztion for my final project. when i execute the coding, it says: 

Warning, initial point [m1 = .900000000000000, m2 = .900000000000000, m3 = .900000000000000, sigma1 = .900000000000000, sigma2 = .900000000000000, sigma3 = .900000000000000] does not satisfy the inequality constraints; trying to find a feasible initial point
Error, (in DirectSearch:-Search) cannot find feasible initial point; last infeasible trial point is [m1 = HFloat(0.9), m2 = HFloat(1.9), m3 = HFloat(1.9), sigma1 = HFloat(1.9), sigma2 = HFloat(0.9), sigma3 = HFloat(0.9)].

here my full coding:


restart;
with(linalg);
with(Optimization);
with(DirectSearch);
[BoundedObjective, CompromiseProgramming, DataFit, 

  ExponentialWeightedSum, GlobalOptima, GlobalSearch, Minimax, 

  ModifiedTchebycheff, Search, SolveEquations, WeightedProduct, 

  WeightedSum]
q := Array(0 .. .111, [0.586e-2, 0.67475e-3, 0.52476e-3, 0.419785e-3, 0.354814e-3, 0.324859e-3, 0.319948e-3, 0.31e-3, 0.295013e-3, 0.28e-3, 0.259974e-3, 0.25496e-3, 0.26e-3, 0.270027e-3, 0.280026e-3, 0.284987e-3, 0.274934e-3, 0.279893e-3, 0.294824e-3, 0.324765e-3, 0.374672e-3, 0.439555e-3, 0.509439e-3, 0.56934e-3, 0.60923e-3, 0.634201e-3, 0.634233e-3, 0.624319e-3, 0.61442e-3, 0.624491e-3, 0.6495e-3, 0.684465e-3, 0.714452e-3, 0.729463e-3, 0.749451e-3, 0.789387e-3, 0.864347e-3, 0.964331e-3, 0.1064261e-2, 0.119413e-2, 0.1333876e-2, 0.150354e-2, 0.1683293e-2, 0.1883023e-2, 0.2102591e-2, 0.2356929e-2, 0.2656062e-2, 0.3004981e-2, 0.3403637e-2, 0.3841874e-2, 0.4349482e-2, 0.4931339e-2, 0.5587622e-2, 0.6293742e-2, 0.7044467e-2, 0.780967e-2, 0.8563865e-2, 0.9291539e-2, 0.9987678e-2, 0.10677379e-1, 0.1139506e-1, 0.12195833e-1, 0.13108819e-1, 0.14262384e-1, 0.15557439e-1, 0.17002799e-1, 0.18574592e-1, 0.20263971e-1, 0.22103523e-1, 0.24043337e-1, 0.26165465e-1, 0.28548531e-1, 0.31489568e-1, 0.3458584e-1, 0.37891201e-1, 0.41511365e-1, 0.45444649e-1, 0.49867283e-1, 0.54877269e-1, 0.60478547e-1, 0.66532745e-1, 0.7397857e-1, 0.81900977e-1, 0.90034141e-1, 0.98921462e-1, .108746024, .11839162, .130752255, .14368791, .157201822, .169770195, .177233464, .188823516, .204733363, .224639756, .249617627, .273409998, .298674064, .304378882, .319592743, .342821691, .367512872, .395257421, .426270115, .460795318, .498484563, .537417055, .577542345, .618093683, .660441764, .703856817, 1]);
variables = [m1, m2, m3, sigma1, sigma2, sigma3];
        variables = [m1, m2, m3, sigma1, sigma2, sigma3]
psi1 := 0.11480601e-1; psi2 := 0.8890123e-2; psi3 := .979629276;
                          0.011480601
                          0.008890123
                          0.979629276
s1(x):=exp(-(x/(m1))^((m1/(sigma1)))):  s1(x+1):=exp(-((x+1)/(m1))^((m1/(sigma1)))):  s2(x):=1-exp(-(x/(m2))^((-m2/(sigma2)))):  s2(x+1):=1-exp(-((x+1)/(m2))^((-m2/(sigma2)))):  s3(x):=exp(exp(-m3/(sigma3))-exp((x-m3)/(sigma3))): s3(x+1):=exp(exp(-m3/(sigma3))-exp((x+1-m3)/(sigma3)))  :  s(x):=(psi1*s1(x)+psi2*s2(x)+psi3*s3(x)): s(x+1):=(psi1*s1(x+1)+psi2*s2(x+1)+psi3*s3(x+1)):   qtopi(x):=1-((s(x+1))/(s(x))):
fungsikerugian := add((1-qtopi(x)/q[x])^2, x = 0 .. .110);
for x from 0 to 111 do c1[x] := 0 <= s1(x+1); c2[x] := 1 >= s1(x+1); c3[x] := 0 <= s2(x+1); c4[x] := 1 >= s2(x+1); c5[x] := 0 <= s3(x+1); c6[x] := 1 >= s3(x+1); c7[x] := 0 <= s1(x); c8[x] := 1 >= s1(x); c9[x] := 0 <= s2(x); c10[x] := 1 >= s2(x); c11[x] := 0 <= s3(x); c12[x] := 1 >= s3(x); c13[x] := 0 <= s(x); c14[x] := 1 >= s(x); c15[x] := 0 <= s(x+1); c16[x] := 1 >= s(x+1); c17[x] := 0 <= qtopi(x); c18[x] := 1 >= qtopi(x) end do;

constr := {m2 > sigma2, m3 > sigma3, seq(c1[x], x = 0 .. .110), seq(c10[x], x = 0 .. .111), seq(c11[x], x = 0 .. .111), seq(c12[x], x = 0 .. .111), seq(c13[x], x = 0 .. .111), seq(c14[x], x = 0 .. .111), seq(c15[x], x = 0 .. .110), seq(c16[x], x = 0 .. .110), seq(c17[x], x = 0 .. .110), seq(c18[x], x = 0 .. .110), seq(c2[x], x = 0 .. .110), seq(c3[x], x = 0 .. .110), seq(c4[x], x = 0 .. .110), seq(c5[x], x = 0 .. .110), seq(c6[x], x = 0 .. .110), seq(c7[x], x = 0 .. .111), seq(c8[x], x = 0 .. .111), seq(c9[x], x = 0 .. .111), m1 < sigma1, 0 <= m1 and m1 < 17, 17 <= m2 and m2 < 33, 33 <= m3 and m3 < 111};
solusi := DirectSearch:-Search(fungsikerugian, constr, assume = positive, evaluationlimit = 5555);
Warning, initial point [m1 = .900000000000000, m2 = .900000000000000, m3 = .900000000000000, sigma1 = .900000000000000, sigma2 = .900000000000000, sigma3 = .900000000000000] does not satisfy the inequality constraints; trying to find a feasible initial point
Error, (in DirectSearch:-Search) cannot find feasible initial point; last infeasible trial point is [m1 = HFloat(0.9), m2 = HFloat(1.9), m3 = HFloat(1.9), sigma1 = HFloat(1.9), sigma2 = HFloat(0.9), sigma3 = HFloat(0.9)]

thanks before

Please Wait...