From: Michael Orlitzky <michael@orlitzky.com>
Date: Sun, 21 Jul 2013 14:14:19 +0000 (-0400)
Subject: Add another example for forward_substitute.
X-Git-Url: http://gitweb.michael.orlitzky.com/?p=numerical-analysis.git;a=commitdiff_plain;h=4a1589178a2928710f13e0cb18efb34ed8c0d883

Add another example for forward_substitute.
---

diff --git a/src/Linear/System.hs b/src/Linear/System.hs
index 58d8d78..f49c652 100644
--- a/src/Linear/System.hs
+++ b/src/Linear/System.hs
@@ -31,6 +31,11 @@ import qualified Algebra.Field as Field
 --   >>> forward_substitute m b
 --   ((1.0),(0.0))
 --
+--   >>> let m = fromList [[4,0],[0,2]] :: Mat2 Double
+--   >>> let b = vec2d (2, 1.5 :: Double)
+--   >>> forward_substitute m b
+--   ((0.5),(0.75))
+--
 forward_substitute :: forall a m. (Field.C a, Arity m)
                    => Mat m m a
                    -> Mat m N1 a