integrate[e^(-b^2(y^2-a^2)^2),{y,0,inf},Assumptions->{a>0,b>0}]