* John W. Eaton <address@hidden> [2005-06-15 11:41]:
On 15-Jun-2005, Andrey Romanenko wrote:
| Maybe we should
| consider creating an an Octave->lp_solve interface to make it automatic.
Octave (2.9.x) now includes an interface to glpk which seems to have
no touble solving the example problem that was posted. So I think
things are improving and the lp in Octave forge could be considered
obsolete soon.
Indeed, I get the following results for the problem presented by Vic
Norton running Octave 2.9.2:
octave-2.9.2:1> A6 = [
1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
0.2 0.5 0.3 0.3 0.3 0.6 0.4 0.1 0.1
0.3 0.4 0.2 0.4 0.3 0.3 0.5 0.3 0.1
0.5 0.1 0.5 0.3 0.4 0.1 0.1 0.6 0.8
];
octave-2.9.2:2> B6 = [1.0; 0.3; 0.3; 0.4];
octave-2.9.2:3> F6 = [4.1 4.3 5.8 6.0 7.6 7.5 7.3 6.9 7.3];
octave-2.9.2:4> LB6 = zeros(9, 1);
octave-2.9.2:5> UB6 = ones(9, 1)*Inf;
octave-2.9.2:6> sol6 = glpk (F6, A6, B6, LB6, UB6, 'SSSS')
sol6 =
0.50000
0.25000
0.25000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
--
Rafael