diff --git a/miscprob.cc b/miscprob.cc
index c3d9b8f5ccd654ccc1cc6a8044915dedf59e2549..8d462aeadd71b37e94f55095785ca79b76b23d54 100644
--- a/miscprob.cc
+++ b/miscprob.cc
@@ -205,7 +205,7 @@ ReturnMatrix mvnrnd(const RowVector& mu, const SymmetricMatrix& covar, int nsamp
 
  float mvnpdf(const RowVector& vals, const RowVector& mu, const SymmetricMatrix& covar)
  {
-   return std::exp(-0.5*((vals-mu).t()*covar.i()*(vals-mu)).AsScalar())/(std::pow(covar.Determinant(),0.5)*std::pow(2*M_PI,vals.Ncols()/2.0));
+   return std::exp(-0.5*((vals-mu)*covar.i()*(vals-mu).t()).AsScalar())/(std::pow(covar.Determinant(),0.5)*std::pow(2*M_PI,vals.Ncols()/2.0));
  }
 
 // ReturnMatrix gammarnd(const int dim1, const int dim2,