/*
	find one feasible solution
	coefficients are given in this program
*/

#include <stdio.h>
#include "opt.h"

/*
  find 	x0, x1, x2 >= 0
  s.t.	-x0 + 2 x1      = 3
	        x1 + x2 = 4
*/

#define	NumEqs	2
#define NumVars	3

main()
{
	int	m = NumEqs;
	int	n = NumVars;
	double	a[NumEqs][NumVars] = {{-1.00, 2.00, 0.00},
			      	      { 0.00, 1.00, 1.00}};
	double	b[NumEqs] = {3.00, 4.00};
	double	eps = 1e-4;
	
	double	ox[NumVars], opt;
	status	st;
	
	matrixfprint(stdout, m, n, a);
	matrixfprint(stdout, m, 1, b);
	
	st = LPfeasiblesolution(m, n, a, b, eps, ox);
	
	if (st == success) {
		puts("--- find one solution ---");
		matrixfprint(stdout, n, 1, ox);
	} else {
		puts("--- no feasible solution ---");
	}
}