sorry , i forget the term so the question is :

int ( q  * exp((1-b)^2* q_z^2 / lambda[1] ) * exp((1-b)^2* q^2 / lambda[2])   /  ((1+b^2*q^2+b^2*q_z^2)*(q^2+q_z^2)) , q_z= -infinity..infinity)  then w.r.t     q  from (0......infinity)

after integration we want to minimize the results , the range is as (b=0..1), ( lambda[1])  =1....100000) ,   (lambda[2])  =1....100000)