Abstract: We have developed a new Einstein-Boltzmann solver for CMB anisotropy to quantify the significance of the non-Gaussianity induced during the non-linear evolution of the primordial fluctuations after the last-scattering surface. In our numerical code, we implemented a 2nd-order line-of-sight formula, so-called “curve”-of-sight formula, which has been developed in our recent work in JCAP 1410 (2014) 051. In this talk, we show the numerical results on the non-Gaussianity of scalar perturbations at the late time. In addition, we are now extending this code to include the tensor modes, so we will also show the preliminary results on the bispectra induced by the tensor modes.