Thanks for your fast input! FWIW, I want to use bounds for my fitting params that is why I need mpfit... JF