This program aims to get inner point from multiple both side linear constraints. The proof before is worked in non extreme condition, but to calculate in the numerical, there's a problem concern with error order and return norm order. The implementation now has a const. multiply order glitch. So to use this implementation, we need to adjust const. multiply before to use. Bug reports are welcomed.