diff --git a/src/ExprNode.cc b/src/ExprNode.cc
index 83136997f0cd0b36e04b6257d84af2d250c08980..de3cd447bf171b7905c38605eebc15fde93e2db7 100644
--- a/src/ExprNode.cc
+++ b/src/ExprNode.cc
@@ -5939,6 +5939,17 @@ TrinaryOpNode::writeOutput(ostream &output, ExprNodeOutputType output_type,
           arg3->writeOutput(output, output_type, temporary_terms, temporary_terms_idxs, tef_terms);
           output << ")/M_SQRT2)))";
         }
+      else if (isJuliaOutput(output_type))
+	{
+	  // Julia API is normcdf(mu, sigma, x) !
+          output << "normcdf(";
+          arg2->writeOutput(output, output_type, temporary_terms, temporary_terms_idxs, tef_terms);
+          output << ",";
+          arg3->writeOutput(output, output_type, temporary_terms, temporary_terms_idxs, tef_terms);
+          output << ",";
+          arg1->writeOutput(output, output_type, temporary_terms, temporary_terms_idxs, tef_terms);
+          output << ")";
+        }
       else
         {
           output << "normcdf(";
@@ -5964,6 +5975,17 @@ TrinaryOpNode::writeOutput(ostream &output, ExprNodeOutputType output_type,
           arg3->writeOutput(output, output_type, temporary_terms, temporary_terms_idxs, tef_terms);
           output << ",2)/2)))";
         }
+      else if (isJuliaOutput(output_type))
+	{
+	  // Julia API is normpdf(mu, sigma, x) !
+          output << "normpdf(";
+          arg2->writeOutput(output, output_type, temporary_terms, temporary_terms_idxs, tef_terms);
+          output << ",";
+          arg3->writeOutput(output, output_type, temporary_terms, temporary_terms_idxs, tef_terms);
+          output << ",";
+          arg1->writeOutput(output, output_type, temporary_terms, temporary_terms_idxs, tef_terms);
+          output << ")";
+        }
       else
         {
           output << "normpdf(";