We present a method to compute the 3-point correlation function for a
general model of inflation driven by a single, minimally coupled
scalar field. Our approach is based on the numerical evaluation of
both the perturbation equations and the integrals which contribute to
the 3-point function. As an example, we use our method to compute the
3 point correlation function for a class of inflation models with a
feature in the potential which have been invoked to fit the "glitch"
in the power spectrum present in both the 2 and 3 year WMAP dataset.
We find this feature generically leads to large scale and shape
dependent non-gaussianities that are potentially observable with the
future Planck satellite. More generally, we use this example to
explain that while adding features to potential can improve the fit to
the 2-point function, these are generically associated with a greatly
enhanced signal at the 3-point level. Moreover, this 3-point signal
will have a very nontrivial shape and scale dependence, which is
correlated with the form of the 2-point function, and may thus lead to
a consistency check on the models of inflation with non-smooth
potentials.

The paper associated with this talk is astro-ph/0611645.