Hello,
I was trying to solve a linear programming problem with SAS and R but got different results.
Linear Programming Problem with SAS
proc optmodel; var x {1..3} >=0; max z =11*x[1]+16*x[2]+15*x[3]; con c1: x[1]+2*x[2]+1.5*x[3] <=12000; con c2: 0.66*x[1]+0.66*x[2]+x[3] <= 4600; con c3: 0.5*x[1]+0.3333*x[2]+0.5*x[3] <=2400; solve with lp;print x[1] x[2] x[3]; quit;
SAS output: Objective Value=100311.708 x1=490.31 x2=5044.7 x3=946.93
With R
> library(lpSolve)
> obj=c(11,16,15) > con=matrix(c(1,2,3/2,2/3,2/3,1,1/2,1/3,1/2),nrow=3,byrow=TRUE) > dir=c("<=","<=","<=") > rhs=c(12000,4600,2400) > lp("max",obj,con,dir,rhs) Success: the objective function is 100200 > lp("max",obj,con,dir,rhs)$solution [1] 600 5100 800
Please help me with this.
... View more