Support evaluation based stuff in optimizer

[johanlofberg]

Something like this would be nice

sdpvar x y
P = optimizer([x <= sin(y)*cos(y)],x^2,sdpsettings('solver','+quadprog'),y,x)
P{pi}