zdt1.mw

restart; 1 

ZDT 1 

 

`:=`(f, proc (x) options operator, arrow; `+`(1, `-`(sqrt(x))) end proc); -1; 'f(x)' = f(x); 1; `:=`(x[min], 0); -1; `:=`(x[max], 1); -1; 'x[min]' = x[min], 'x[max]' = x[max] 

 

f(x) = `+`(1, `-`(`*`(`^`(x, `/`(1, 2)))))
x[min] = 0, x[max] = 1 (1)
 

 

`:=`(delta, proc (x) options operator, arrow; `/`(`*`(sqrt(`+`(`-`(`/`(`*`(diff(f(x), x)), `*`(`+`(1, `*`(`^`(diff(f(x), x), 2))))))))), `*`(int(sqrt(`+`(`-`(diff(f(x), x)))), x = x[min] .. x[max]))) ... 

delta(x) = `+`(`/`(`*`(`/`(3, 2), `*`(`^`(x, `/`(1, 4)))), `*`(`^`(`+`(`*`(4, `*`(x)), 1), `/`(1, 2))))) (2)
 

`:=`(H1, proc (x) options operator, arrow; `+`(x, `/`(`*`(`+`(f(x), `-`(f(x[min])))), `*`(diff(f(x), x)))) end proc); -1 

`:=`(R1, Optimization:-Maximize(H1(x), x = x[min] .. x[max])); -1; 'R[1]' = R1 

R[1] = [3., [x = 1.]] (3)
 

`:=`(fd, unapply(diff(f(x), x), x)); -1; limit(fd(x), x = x[min]) 

undefined (4)