Rev | Line | |
---|
[37] | 1 | function gp |
---|
| 2 | |
---|
| 3 | sdpvar x y |
---|
| 4 | css=set(x>1)+set(y>1)+set(x/y<4)+set(y<8); |
---|
| 5 | css=css+set(x^2/y==1.5); |
---|
| 6 | obj = x+y/x; |
---|
| 7 | sol = solvesdp(css,obj) |
---|
| 8 | |
---|
| 9 | mbg_asserttolequal(sol.problem,0); |
---|
| 10 | mbg_asserttolequal(double(obj),2.04124145231932,1e-5); |
---|
| 11 | mbg_asserttolequal(double([x y]), [ 1.22474487139159 1.00000000000000], 1e-5); |
---|
| 12 | |
---|
| 13 | sdpvar x y z |
---|
| 14 | css=set(x>1)+set(y>1)+set(x/y<4)+set(y<8) + set(x*z == 10) %+ set(1<z<16); |
---|
| 15 | css=css+set(x^2/y==1.5); |
---|
| 16 | obj = x+y/x; |
---|
| 17 | sol = solvesdp(css,obj) |
---|
| 18 | |
---|
| 19 | mbg_asserttolequal(sol.problem,0); |
---|
| 20 | mbg_asserttolequal(double(obj),2.04124145231932,1e-5); |
---|
| 21 | mbg_asserttolequal(double([x y z]), [ 1.22474487139159 1.00000000000000 8.16496580927726], 1e-5); |
---|
| 22 | |
---|
Note: See
TracBrowser
for help on using the repository browser.