10 files changed, 725 insertions, 0 deletions
diff --git a/tests/unit_tests/cbcintlinprog.dia.ref b/tests/unit_tests/cbcintlinprog.dia.ref
new file mode 100644
index 0000000..7ebbe4e
--- /dev/null
+++ b/tests/unit_tests/cbcintlinprog.dia.ref
@@ -0,0 +1,71 @@
+// Copyright (C) 2017 - IIT Bombay - FOSSEE
+//
+// Author: Georgey John
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+// <-- JVM NOT MANDATORY -->
+// <-- ENGLISH IMPOSED -->
+
+
+//
+// assert_close --
+//   Returns 1 if the two real matrices computed and expected are close,
+//   i.e. if the relative distance between computed and expected is lesser than epsilon.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+function flag = assert_close ( computed, expected, epsilon )
+  if expected==0.0 then
+    shift = norm(computed-expected);
+  else
+    shift = norm(computed-expected)/norm(expected);
+  end
+//  if shift < epsilon then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+    flag = assert_checktrue ( shift < epsilon );
+endfunction
+//
+// assert_equal --
+//   Returns 1 if the two real matrices computed and expected are equal.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+//function flag = assert_equal ( computed , expected )
+//  if computed==expected then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+//endfunction
+
+
+//Integer Linear program with all constraint types
+c=[-100,-150]'
+intcon=[1,2]
+A=[8,4;15,30]
+b=[40,200]
+Aeq=[1,1]
+beq=[7]
+lb=[0,0]
+ub=[]
+[xopt,fopt,exitflag,output]=cbcintlinprog(c, intcon, A, b, Aeq, beq, lb, ub)
+
+assert_close ( xopt , [ 1 6 ]' , 0.0005 );
+assert_close ( fopt , [ -1000 ]' , 0.0005 );
+assert_checkequal( exitflag , 0 );
+printf("Test Successful");
+
diff --git a/tests/unit_tests/cbcintlinprog.tst b/tests/unit_tests/cbcintlinprog.tst
new file mode 100644
index 0000000..7ebbe4e
--- /dev/null
+++ b/tests/unit_tests/cbcintlinprog.tst
@@ -0,0 +1,71 @@
+// Copyright (C) 2017 - IIT Bombay - FOSSEE
+//
+// Author: Georgey John
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+// <-- JVM NOT MANDATORY -->
+// <-- ENGLISH IMPOSED -->
+
+
+//
+// assert_close --
+//   Returns 1 if the two real matrices computed and expected are close,
+//   i.e. if the relative distance between computed and expected is lesser than epsilon.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+function flag = assert_close ( computed, expected, epsilon )
+  if expected==0.0 then
+    shift = norm(computed-expected);
+  else
+    shift = norm(computed-expected)/norm(expected);
+  end
+//  if shift < epsilon then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+    flag = assert_checktrue ( shift < epsilon );
+endfunction
+//
+// assert_equal --
+//   Returns 1 if the two real matrices computed and expected are equal.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+//function flag = assert_equal ( computed , expected )
+//  if computed==expected then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+//endfunction
+
+
+//Integer Linear program with all constraint types
+c=[-100,-150]'
+intcon=[1,2]
+A=[8,4;15,30]
+b=[40,200]
+Aeq=[1,1]
+beq=[7]
+lb=[0,0]
+ub=[]
+[xopt,fopt,exitflag,output]=cbcintlinprog(c, intcon, A, b, Aeq, beq, lb, ub)
+
+assert_close ( xopt , [ 1 6 ]' , 0.0005 );
+assert_close ( fopt , [ -1000 ]' , 0.0005 );
+assert_checkequal( exitflag , 0 );
+printf("Test Successful");
+
diff --git a/tests/unit_tests/intfminbnd.dia.ref b/tests/unit_tests/intfminbnd.dia.ref
new file mode 100644
index 0000000..50df69c
--- /dev/null
+++ b/tests/unit_tests/intfminbnd.dia.ref
@@ -0,0 +1,73 @@
+// Copyright (C) 2017 - IIT Bombay - FOSSEE
+//
+// Author: Georgey John
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+// <-- JVM NOT MANDATORY -->
+// <-- ENGLISH IMPOSED -->
+
+
+//
+// assert_close --
+//   Returns 1 if the two real matrices computed and expected are close,
+//   i.e. if the relative distance between computed and expected is lesser than epsilon.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+function flag = assert_close ( computed, expected, epsilon )
+  if expected==0.0 then
+    shift = norm(computed-expected);
+  else
+    shift = norm(computed-expected)/norm(expected);
+  end
+//  if shift < epsilon then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+    flag = assert_checktrue ( shift < epsilon );
+endfunction
+//
+// assert_equal --
+//   Returns 1 if the two real matrices computed and expected are equal.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+//function flag = assert_equal ( computed , expected )
+//  if computed==expected then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+//endfunction
+
+
+//Bound constrained optimization
+function y=f(x)
+y= x(1)^2 + x(2)^2 + 19*x(1) +5*x(2);
+endfunction
+
+x1=[4.5,7];
+
+x2=[10,9];
+
+intcon=[1 2];
+
+//Calling fminunc function
+[xopt,fopt,exitflag,gradient,hessian]=intfminbnd(f,intcon,x1,x2)
+
+assert_close ( xopt , [ 5 7 ]' , 0.0005 );
+assert_close ( fopt , [ 204 ]' , 0.0005 );
+assert_checkequal( double(exitflag) , 0 );
+printf("Test Successful");
+\ No newline at end of file
diff --git a/tests/unit_tests/intfminbnd.tst b/tests/unit_tests/intfminbnd.tst
new file mode 100644
index 0000000..50df69c
--- /dev/null
+++ b/tests/unit_tests/intfminbnd.tst
@@ -0,0 +1,73 @@
+// Copyright (C) 2017 - IIT Bombay - FOSSEE
+//
+// Author: Georgey John
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+// <-- JVM NOT MANDATORY -->
+// <-- ENGLISH IMPOSED -->
+
+
+//
+// assert_close --
+//   Returns 1 if the two real matrices computed and expected are close,
+//   i.e. if the relative distance between computed and expected is lesser than epsilon.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+function flag = assert_close ( computed, expected, epsilon )
+  if expected==0.0 then
+    shift = norm(computed-expected);
+  else
+    shift = norm(computed-expected)/norm(expected);
+  end
+//  if shift < epsilon then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+    flag = assert_checktrue ( shift < epsilon );
+endfunction
+//
+// assert_equal --
+//   Returns 1 if the two real matrices computed and expected are equal.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+//function flag = assert_equal ( computed , expected )
+//  if computed==expected then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+//endfunction
+
+
+//Bound constrained optimization
+function y=f(x)
+y= x(1)^2 + x(2)^2 + 19*x(1) +5*x(2);
+endfunction
+
+x1=[4.5,7];
+
+x2=[10,9];
+
+intcon=[1 2];
+
+//Calling fminunc function
+[xopt,fopt,exitflag,gradient,hessian]=intfminbnd(f,intcon,x1,x2)
+
+assert_close ( xopt , [ 5 7 ]' , 0.0005 );
+assert_close ( fopt , [ 204 ]' , 0.0005 );
+assert_checkequal( double(exitflag) , 0 );
+printf("Test Successful");
+\ No newline at end of file
diff --git a/tests/unit_tests/intfmincon.dia.ref b/tests/unit_tests/intfmincon.dia.ref
new file mode 100644
index 0000000..b2f266d
--- /dev/null
+++ b/tests/unit_tests/intfmincon.dia.ref
@@ -0,0 +1,76 @@
+// Copyright (C) 2017 - IIT Bombay - FOSSEE
+//
+// Author: Georgey John
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+// <-- JVM NOT MANDATORY -->
+// <-- ENGLISH IMPOSED -->
+
+
+//
+// assert_close --
+//   Returns 1 if the two real matrices computed and expected are close,
+//   i.e. if the relative distance between computed and expected is lesser than epsilon.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+function flag = assert_close ( computed, expected, epsilon )
+  if expected==0.0 then
+    shift = norm(computed-expected);
+  else
+    shift = norm(computed-expected)/norm(expected);
+  end
+//  if shift < epsilon then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+    flag = assert_checktrue ( shift < epsilon );
+endfunction
+//
+// assert_equal --
+//   Returns 1 if the two real matrices computed and expected are equal.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+//function flag = assert_equal ( computed , expected )
+//  if computed==expected then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+//endfunction
+
+
+//Constrained optimization
+function [y,dy]=f(x)
+y=-x(1)^2-x(2)^2-x(2)/3;
+endfunction
+
+x0=[0 , 0];
+intcon = [1]
+A=[1,1 ; 1,1/4 ; 1,-1 ; -1/4,-1 ; -1,-1 ; -1,1];
+b=[2;1;2;1;-1;2];
+Aeq=[1,1];
+beq=[2];
+lb=[];
+ub=[];
+nlc=[];
+
+[xopt,fopt,exitflag,grad,hessian] =intfmincon(f, x0,intcon,A,b,Aeq,beq,lb,ub,nlc)
+
+assert_close ( xopt , [ 0 2 ]' , 0.0005 );
+assert_close ( fopt , [ -4.6666 ]' , 0.0005 );
+assert_checkequal( double(exitflag) , 0 );
+printf("Test Successful");
+\ No newline at end of file
diff --git a/tests/unit_tests/intfmincon.tst b/tests/unit_tests/intfmincon.tst
new file mode 100644
index 0000000..144ff2b
--- /dev/null
+++ b/tests/unit_tests/intfmincon.tst
@@ -0,0 +1,76 @@
+// Copyright (C) 2017 - IIT Bombay - FOSSEE
+//
+// Author: Georgey John
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+// <-- JVM NOT MANDATORY -->
+// <-- ENGLISH IMPOSED -->
+
+
+//
+// assert_close --
+//   Returns 1 if the two real matrices computed and expected are close,
+//   i.e. if the relative distance between computed and expected is lesser than epsilon.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+function flag = assert_close ( computed, expected, epsilon )
+  if expected==0.0 then
+    shift = norm(computed-expected);
+  else
+    shift = norm(computed-expected)/norm(expected);
+  end
+//  if shift < epsilon then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+    flag = assert_checktrue ( shift < epsilon );
+endfunction
+//
+// assert_equal --
+//   Returns 1 if the two real matrices computed and expected are equal.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+//function flag = assert_equal ( computed , expected )
+//  if computed==expected then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+//endfunction
+
+
+//Constrained optimization
+function [y]=f(x)
+y=-x(1)^2-x(2)^2-x(2)/3;
+endfunction
+
+x0=[0 , 0];
+intcon = [1]
+A=[1,1 ; 1,1/4 ; 1,-1 ; -1/4,-1 ; -1,-1 ; -1,1];
+b=[2;1;2;1;-1;2];
+Aeq=[1,1];
+beq=[2];
+lb=[];
+ub=[];
+nlc=[];
+
+[xopt,fopt,exitflag,grad,hessian] =intfmincon(f, x0,intcon,A,b,Aeq,beq,lb,ub,nlc)
+
+assert_close ( xopt , [ 0 2 ]' , 0.0005 );
+assert_close ( fopt , [ -4.6666 ]' , 0.0005 );
+assert_checkequal( double(exitflag) , 0 );
+printf("Test Successful");
+\ No newline at end of file
diff --git a/tests/unit_tests/intfminunc.dia.ref b/tests/unit_tests/intfminunc.dia.ref
new file mode 100644
index 0000000..aa63e91
--- /dev/null
+++ b/tests/unit_tests/intfminunc.dia.ref
@@ -0,0 +1,70 @@
+// Copyright (C) 2017 - IIT Bombay - FOSSEE
+//
+// Author: Georgey John
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+// <-- JVM NOT MANDATORY -->
+// <-- ENGLISH IMPOSED -->
+
+
+//
+// assert_close --
+//   Returns 1 if the two real matrices computed and expected are close,
+//   i.e. if the relative distance between computed and expected is lesser than epsilon.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+function flag = assert_close ( computed, expected, epsilon )
+  if expected==0.0 then
+    shift = norm(computed-expected);
+  else
+    shift = norm(computed-expected)/norm(expected);
+  end
+//  if shift < epsilon then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+    flag = assert_checktrue ( shift < epsilon );
+endfunction
+//
+// assert_equal --
+//   Returns 1 if the two real matrices computed and expected are equal.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+//function flag = assert_equal ( computed , expected )
+//  if computed==expected then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+//endfunction
+
+
+//Unconstrained optimization
+function y=f(x)
+y= 4*x(1)^2 + 5*x(2)^2 + 6*x(1);
+endfunction
+
+x0=[5,5];
+intcon = [1,2]
+//Calling fminunc function
+[xopt,fopt,exitflag,gradient,hessian]=intfminunc(f,x0,intcon)
+
+assert_close ( xopt , [ -1 0 ]' , 0.0005 );
+assert_close ( fopt , [ -2 ]' , 0.0005 );
+assert_checkequal( double(exitflag) , 0 );
+printf("Test Successful");
+
diff --git a/tests/unit_tests/intfminunc.tst b/tests/unit_tests/intfminunc.tst
new file mode 100644
index 0000000..42171c1
--- /dev/null
+++ b/tests/unit_tests/intfminunc.tst
@@ -0,0 +1,70 @@
+// Copyright (C) 2017 - IIT Bombay - FOSSEE
+//
+// Author: Georgey John
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+// <-- JVM NOT MANDATORY -->
+// <-- ENGLISH IMPOSED -->
+
+
+//
+// assert_close --
+//   Returns 1 if the two real matrices computed and expected are close,
+//   i.e. if the relative distance between computed and expected is lesser than epsilon.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+function flag = assert_close ( computed, expected, epsilon )
+  if expected==0.0 then
+    shift = norm(computed-expected);
+  else
+    shift = norm(computed-expected)/norm(expected);
+  end
+//  if shift < epsilon then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+    flag = assert_checktrue ( shift < epsilon );
+endfunction
+//
+// assert_equal --
+//   Returns 1 if the two real matrices computed and expected are equal.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+//function flag = assert_equal ( computed , expected )
+//  if computed==expected then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+//endfunction
+
+
+//Unconstrained optimization
+function y=f(x)
+y= 4*x(1)^2 + 5*x(2)^2 + 6*x(1);
+endfunction
+
+x0=[5,5];
+intcon = [1,2]
+//Calling fminunc function
+[xopt,fopt,exitflag,gradient,hessian]=intfminunc(f,x0,intcon)
+
+assert_close ( xopt , [ -1 0 ]' , 0.0005 );
+assert_close ( fopt , [ -2 ]' , 0.0005 );
+assert_checkequal( exitflag , 0 );
+printf("Test Successful");
+
diff --git a/tests/unit_tests/intqpipopt.dia.ref b/tests/unit_tests/intqpipopt.dia.ref
new file mode 100644
index 0000000..a3b0153
--- /dev/null
+++ b/tests/unit_tests/intqpipopt.dia.ref
@@ -0,0 +1,73 @@
+// Copyright (C) 2017 - IIT Bombay - FOSSEE
+//
+// Author: Georgey John
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+// <-- JVM NOT MANDATORY -->
+// <-- ENGLISH IMPOSED -->
+
+
+//
+// assert_close --
+//   Returns 1 if the two real matrices computed and expected are close,
+//   i.e. if the relative distance between computed and expected is lesser than epsilon.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+function flag = assert_close ( computed, expected, epsilon )
+  if expected==0.0 then
+    shift = norm(computed-expected);
+  else
+    shift = norm(computed-expected)/norm(expected);
+  end
+//  if shift < epsilon then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+    flag = assert_checktrue ( shift < epsilon );
+endfunction
+//
+// assert_equal --
+//   Returns 1 if the two real matrices computed and expected are equal.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+//function flag = assert_equal ( computed , expected )
+//  if computed==expected then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+//endfunction
+
+
+//Quadratic programming optimization
+
+H = [1 -1; -1 2];
+f = [-2; -6];
+
+A = [1 1; -1 2; 2 1];
+b = [2; 2; 3];
+lb=[0,0];
+ub=[%inf, %inf];
+intcon = [1 2];
+
+[xopt,fopt,status,output]=intqpipopt(H,f,intcon,A,b,[],[],lb,ub)
+
+
+assert_close ( xopt , [ 1 1 ]' , 0.0005 );
+assert_close ( fopt , [ -7.5 ]' , 0.0005 );
+assert_checkequal( double(exitflag) , 0 );
+printf("Test Successful");
+\ No newline at end of file
diff --git a/tests/unit_tests/intqpipopt.tst b/tests/unit_tests/intqpipopt.tst
new file mode 100644
index 0000000..fae106d
--- /dev/null
+++ b/tests/unit_tests/intqpipopt.tst
@@ -0,0 +1,72 @@
+// Copyright (C) 2017 - IIT Bombay - FOSSEE
+//
+// Author: Georgey John
+// Organization: FOSSEE, IIT Bombay
+// Email: toolbox@scilab.in
+//
+// This file must be used under the terms of the CeCILL.
+// This source file is licensed as described in the file COPYING, which
+// you should have received as part of this distribution.  The terms
+// are also available at
+// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt
+
+// <-- JVM NOT MANDATORY -->
+// <-- ENGLISH IMPOSED -->
+
+
+//
+// assert_close --
+//   Returns 1 if the two real matrices computed and expected are close,
+//   i.e. if the relative distance between computed and expected is lesser than epsilon.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+function flag = assert_close ( computed, expected, epsilon )
+  if expected==0.0 then
+    shift = norm(computed-expected);
+  else
+    shift = norm(computed-expected)/norm(expected);
+  end
+//  if shift < epsilon then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+    flag = assert_checktrue ( shift < epsilon );
+endfunction
+//
+// assert_equal --
+//   Returns 1 if the two real matrices computed and expected are equal.
+// Arguments
+//   computed, expected : the two matrices to compare
+//   epsilon : a small number
+//
+//function flag = assert_equal ( computed , expected )
+//  if computed==expected then
+//    flag = 1;
+//  else
+//    flag = 0;
+//  end
+//  if flag <> 1 then pause,end
+//endfunction
+
+
+//Quadratic programming optimization
+
+H = [1 -1; -1 2];
+f = [-2; -6];
+
+A = [1 1; -1 2; 2 1];
+b = [2; 2; 3];
+lb=[0,0];
+ub=[%inf, %inf];
+intcon = [1 2];
+
+[xopt,fopt,status,output]=intqpipopt(H,f,intcon,A,b,[],[],lb,ub)
+
+assert_close ( xopt , [ 1 1 ]' , 0.0005 );
+assert_close ( fopt , [ -7.5 ]' , 0.0005 );
+assert_checkequal( double(exitflag) , 0 );
+printf("Test Successful");
+\ No newline at end of file