  In this paper we show how to approximate the transition density of a CARMA(p,
q) model driven by means of a time changed Brownian Motion based on the
Gauss-Laguerre quadrature. We then provide an analytical formula for option
prices when the log price follows a CARMA(p, q) model. We also propose an
estimation procedure based on the approximated likelihood density.
