Two corrections for the same email. Must be Friday...
First, for variable bounds in an LP, the shadow price equals the reduced cost on the variable. So there's no need to add constraints as I suggested before; just retrieve the reduced cost on the variable.
Now, looking at your second email, it seems that you're looking for economic values rather than reduced costs. You can try to use the barrier algorithm (without crossover), but ultimately you should look at other methods to compute this - subgradients, perturbations, etc.