search for: eval_g_ineq

...return( ) } #constraint function eval_g0 <- function(x) { return( ) } # objective function eval_f0 <- function(x) { romberg(f, 0.5, 0.5001)} ARL1 <- nloptr( x0=c(0.653), eval_f=eval_f0, lb = c(0), ub = c(6), eval_g_ineq = eval_g0, opts = list("algorithm"="NLOPT_LN_COBYLA", "maxeval"=1000), ) print( ) Thanks [[alternative HTML version deleted]]

...u1) z1<-t(y1) ARLV1<-K1%*%one1 ARLV21<-t(ARLV1) return( abs((ARLV21%*%y1)-500)-1) } eval_f0 <- function(x) romberg(function(s) f(x, s), 0, 6) of <- nloptr( x0=c(0.653,0.09), eval_f=eval_f0, lb = c(0,0), ub = c(6,1), eval_g_ineq = eval_g0, opts = list("algorithm"="NLOPT_LN_COBYLA", "maxeval"=1000), ) print(of) Regards, Aya [[alternative HTML version deleted]]

...tion(x) abs(f(x,0)-500)-1 #objective function eval_f0 <- function(x) romberg(function(y) ((1+y*y)*f(x,y)/6),0,6, tol = .Machine$double.eps^(2/3))$value obj <- nloptr( x0=c(0.653,0.09), eval_f=eval_f0, lb = c(0,0), ub = c(6,1), eval_g_ineq = eval_g0, opts = list("algorithm"="NLOPT_LN_COBYLA", "maxeval"=600)) obj(ARL1) Thanks, Emara [[alternative HTML version deleted]]