polynomial equality constraints at certain points

[Eike F]

Hi!



I have a polynomial P in 4 variables over which I would like to optimise. One constraint is, that P(v)==1 for a certain vector v = [v1;v2;v3;v4], e.g. v=[1;1;1;1].

Is there a way to state this constraint directly, or more generally: Is there a command to evaluate a polynomial at certain points?
If not, can you help me to state it in a different manner? If I understand it correctly, I could state it like:

coefficients(P-1- [c1,c2,c3,c4]*[x(1)-1;x(2)-1;x(3)-1;x(4)-1;],x) ==0

where c1 - c4 are additional parameterised polynomials. Is that right?

Best,

Eike

[Johan Löfberg]

replace(P,v,[1 1 1 1]) == 1

would be the easiest I think

[Eike F]

Thank you very much!