In reply to Azam Nabizadeh.

This can be a tricky problem. For linear constraints and boundaries, you can use the NLPFEA subroutine.