Index: /issm/trunk-jpl/src/m/qmu/gaussian_pdf.m
===================================================================
--- /issm/trunk-jpl/src/m/qmu/gaussian_pdf.m	(revision 24315)
+++ /issm/trunk-jpl/src/m/qmu/gaussian_pdf.m	(revision 24315)
@@ -0,0 +1,3 @@
+function pdf=gaussian_pdf(x,average,sigma)
+
+pdf=1/(sqrt(2*pi*sigma^2))*exp(-(x-average).^2/(2*sigma^2));
Index: /issm/trunk-jpl/src/m/qmu/lognormal_pdf.m
===================================================================
--- /issm/trunk-jpl/src/m/qmu/lognormal_pdf.m	(revision 24315)
+++ /issm/trunk-jpl/src/m/qmu/lognormal_pdf.m	(revision 24315)
@@ -0,0 +1,4 @@
+function pdf=lognormal_pdf(x,mu,sigma)
+
+pdf=1./(x*sigma*sqrt(2*pi)) .* exp( -(log(x)-mu).^2/(2*sigma^2)  );
+