Line | |
---|
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.