  We carry out a numerical calculation of the bispectrum in generalised
trajectories of canonical, single--field inflation. The trajectories are
generated in the Hamilton-Jacobi (HJ) formalism based on Hubble Slow Roll (HSR)
parameters. The calculation allows generally shape and scale dependent
bispectra, or dimensionless $f_{NL}$, in the out-of-slow-roll regime. The
distributions of $f_{NL}$ for various shapes and HSR proposals are shown as an
example of how this procedure can be used within the context of Monte Carlo
exploration of inflationary trajectories. We also show how allowing
out-of-slow-roll behaviour can lead to a bispectrum that is relatively large
for equilateral shapes.
