In this paper, we obtain analyticity of the inhomogeneous Navier Stokes equations. The main idea is to use the exponential operator e(phi(t)vertical bar D vertical bar), where phi(t) = delta - theta(t), delta > 0 is the analyticity radius of (rho(0) - 1, u(0)), and vertical bar D vertical bar is the differential operator whose symbol is given by parallel to xi parallel to(l1) We will show that for sufficiently small initial data, solutions are analytic globally in time in critical Besov spaces.