To guide the optimization towards a plausible response function, we initialize \(R(x)\)  to \(x^{0.45}\)  and add the following additional constraints based on [1]}:
\(R(0) & = 0 \\R(1) & = 1 \\R^{\prime \prime }(x) & = 0\) 
