Dear all, we want to find the maximum of the function G over [0,1]x[0,1]. The code below gives 0 which is clearly wrong.
Using instead sage: minimize_constrained(-G, [[0, 1], [0, 1]],[1/2, 1/2]) seems to provide the right answer... var('x,y') c1(x,y)=x c2(x,y)=y c3(x,y)=1-x c4(x,y)=1-y G(x,y) = -(56*x^2*y^2 - 4*x^2*y - 4*x*y^2 - 4*x^2 - x*y - 4*y^2 - x - y - 1)*(x - 1)*(y - 1) M = minimize_constrained(-G(x,y), [c1(x,y),c2(x,y),c3(x,y),c4(x,y)],[0.5,0.5]) print(G(M[0],M[1])) -- You received this message because you are subscribed to the Google Groups "sage-support" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-support+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sage-support/cd93a61b-2547-4365-bd0b-9394272e98c5n%40googlegroups.com.