initial add

4 files changed, 115 insertions, 0 deletions

diff --git a/demos/exam.sce b/demos/exam.sce
new file mode 100644
index 0000000..e054eff
--- /dev/null
+++ b/demos/exam.sce
@@ -0,0 +1,5 @@
+function y = f(x) 
+  y = 3*x(1)^2 + 2*x(1)*x(2) + x(2)^2 - 4*x(1) + 5*x(2) ;
+endfunction
+
+[xval, fval, status, gradient, hessian] = intfminbnd(f,[1], [1 2],[2 6])
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 =========//
diff --git a/demos/intfminunc.dem.sce b/demos/intfminunc.dem.sce
new file mode 100644
index 0000000..97cbb2d
--- /dev/null
+++ b/demos/intfminunc.dem.sce
@@ -0,0 +1,49 @@
+mode(1)
+//
+// Demo of intfminunc.sci
+//
+
+//Find x in R^2 such that it minimizes the Rosenbrock function
+//f = 100*(x2 - x1^2)^2 + (1-x1)^2
+//Objective function to be minimised
+function y= f(x)
+y= 100*(x(2) - x(1)^2)^2 + (1-x(1))^2;
+endfunction
+//Starting point
+x0=[-1,2];
+intcon = [2]
+//Options
+options=list("MaxIter", [1500], "CpuTime", [500]);
+//Calling
+[xopt,fopt,exitflag,gradient,hessian]=intfminunc(f,x0,intcon,options)
+// Press ENTER to continue
+halt()   // Press return to continue
+ 
+//Find x in R^2 such that the below function is minimum
+//f = x1^2 + x2^2
+//Objective function to be minimised
+function y= f(x)
+y= x(1)^2 + x(2)^2;
+endfunction
+//Starting point
+x0=[2,1];
+intcon = [1];
+[xopt,fopt]=intfminunc(f,x0,intcon)
+// Press ENTER to continue
+halt()   // Press return to continue
+ 
+//The below problem is an unbounded problem:
+//Find x in R^2 such that the below function is minimum
+//f = - x1^2 - x2^2
+//Objective function to be minimised
+function [y,g,h] = f(x)
+y = -x(1)^2 - x(2)^2;
+g = [-2*x(1),-2*x(2)];
+h = [-2,0;0,-2];
+endfunction
+//Starting point
+x0=[2,1];
+intcon = [1]
+options = list("gradobj","ON","hessian","on");
+[xopt,fopt,exitflag,gradient,hessian]=intfminunc(f,x0,intcon,options)
+//========= E N D === O F === D E M O =========//
diff --git a/demos/intqpipopt.dem.sce b/demos/intqpipopt.dem.sce
new file mode 100644
index 0000000..7c63476
--- /dev/null
+++ b/demos/intqpipopt.dem.sce
@@ -0,0 +1,6 @@
+mode(1)
+//
+// Demo of intqpipopt.sci
+//
+
+//========= E N D === O F === D E M O =========//