From 7ff36b680b38d73e48db3fc770cadc818ac4af4d Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?St=C3=A9phane=20Adjemian=20=28Charybdis=29?=
 <stephane.adjemian@univ-lemans.fr>
Date: Mon, 19 Dec 2011 15:54:30 +0100
Subject: [PATCH] Added unitary test.

---
 matlab/gauss_hermite_weights_and_nodes.m | 17 ++++++++++++++++-
 1 file changed, 16 insertions(+), 1 deletion(-)

diff --git a/matlab/gauss_hermite_weights_and_nodes.m b/matlab/gauss_hermite_weights_and_nodes.m
index 7db89ca26..5946ac25e 100644
--- a/matlab/gauss_hermite_weights_and_nodes.m
+++ b/matlab/gauss_hermite_weights_and_nodes.m
@@ -105,4 +105,19 @@ nodes = sqrt(2)*nodes;
 %$ t(2) = dyn_assert(ones(1,2),variance,1e-12);
 %$ t(3) = dyn_assert(zeros(1,2),expectation,1e-12);
 %$ T = all(t);
-%@eof:3
\ No newline at end of file
+%@eof:3
+
+%@test:4
+%$ n = 9; sigma = .1;
+%$ [nodes,weights] = gauss_hermite_weights_and_nodes(n);
+%$
+%$ sum_of_weights = sum(weights);
+%$ expectation = sum(weights.*nodes*.1); 
+%$ variance = sum(weights.*((nodes*.1).^2));
+%$
+%$ % Check the results.
+%$ t(1) = dyn_assert(1.0,sum_of_weights,1e-12);
+%$ t(2) = dyn_assert(.01,variance,1e-12);
+%$ t(3) = dyn_assert(0.0,expectation,1e-12);
+%$ T = all(t);
+%@eof:4