diff options
author | Georgey | 2017-07-05 11:43:47 +0530 |
---|---|---|
committer | Georgey | 2017-07-05 11:43:47 +0530 |
commit | ecdd7ba6762a15bdb867be900e7128201be2e690 (patch) | |
tree | dc36ba4879b83eb45d042d76c0650ab7ae59144b /demos/intfminbnd.dem.sce | |
parent | 66089674c189f557b401f2ad0cf6b35354caadfa (diff) | |
download | symphony-ecdd7ba6762a15bdb867be900e7128201be2e690.tar.gz symphony-ecdd7ba6762a15bdb867be900e7128201be2e690.tar.bz2 symphony-ecdd7ba6762a15bdb867be900e7128201be2e690.zip |
Added demo files
Diffstat (limited to 'demos/intfminbnd.dem.sce')
-rw-r--r-- | demos/intfminbnd.dem.sce | 55 |
1 files changed, 55 insertions, 0 deletions
diff --git a/demos/intfminbnd.dem.sce b/demos/intfminbnd.dem.sce new file mode 100644 index 0000000..adbc9fa --- /dev/null +++ b/demos/intfminbnd.dem.sce @@ -0,0 +1,55 @@ +mode(1) +// +// Demo of intfminbnd.sci +// + +//Find x in R^6 such that it minimizes: +//f(x)= sin(x1) + sin(x2) + sin(x3) + sin(x4) + sin(x5) + sin(x6) +//-2 <= x1,x2,x3,x4,x5,x6 <= 2 +//Objective function to be minimised +function y=f(x) +y=0 +for i =1:6 +y=y+sin(x(i)); +end +endfunction +//Variable bounds +x1 = [-2, -2, -2, -2, -2, -2]; +x2 = [2, 2, 2, 2, 2, 2]; +intcon = [2 3 4] +//Options +options=list("MaxIter",[1500],"CpuTime", [100]) +[x,fval] =intfminbnd(f ,intcon, x1, x2, options) +// Press ENTER to continue +halt() // Press return to continue + +//Find x in R such that it minimizes: +//f(x)= 1/x^2 +//0 <= x <= 1000 +//Objective function to be minimised +function y=f(x) +y=1/x^2; +endfunction +//Variable bounds +x1 = [0]; +x2 = [1000]; +intcon = [1]; +[x,fval,exitflag,output,lambda] =intfminbnd(f,intcon , x1, x2) +// Press ENTER to continue +halt() // Press return to continue + +//The below problem is an unbounded problem: +//Find x in R^2 such that it minimizes: +//f(x)= -[(x1-1)^2 + (x2-1)^2] +//-inf <= x1,x2 <= inf +//Objective function to be minimised +function y=f(x) +y=-((x(1)-1)^2+(x(2)-1)^2); +endfunction +//Variable bounds +x1 = [-%inf , -%inf]; +x2 = [ %inf , %inf]; +//Options +options=list("MaxIter",[1500],"CpuTime", [100]) +[x,fval,exitflag,output,lambda] =intfminbnd(f,intcon, x1, x2, options) +//========= E N D === O F === D E M O =========// |